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Abstract 

In the so-called "microscopic" models of vehicular traffic, attention is paid ex- 
plicitly to each individual vehicle each of which is represented by a "particle"; the 
nature of the "interactions" among these particles is determined by the way the 
vehicles influence each others' movement. Therefore, vehicular traffic, modeled as a 
system of interacting "particles" driven far from equilibrium, offers the possibility 
to study various fundamental aspects of truly nonequilibrium systems which are of 
current interest in statistical physics. Analytical as well as numerical techniques of 
statistical physics are being used to study these models to understand rich variety 
of physical phenomena exhibited by vehicular traffic. Some of these phenomena, 
observed in vehicular traffic under different circumstances, include transitions from 
one dynamical phase to another, criticality and self-organized criticality, metasta- 
bility and hysteresis, phase-segregation, etc. In this critical review, written from 
the perspective of statistical physics, we explain the guiding principles behind all 
the main theoretical approaches. But we present detailed discussions on the re- 
sults obtained mainly from the so-called " particle- hopping" models, particularly 
emphasizing those which have been formulated in recent years using the language 
of cellular automata. 
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1 Introduction 



The concepts and techniques of statistical physics are being used nowadays 
to study several aspects of complex systems [1] many of which, till a few 
decades ago, used to fall outside the traditional domain of physical systems 
[2]. Physical-, chemical-, earth-, biological- and social-sciences as well as tech- 
nology meet at this frontier area of inter-disciplinary research. Flow of vehic- 
ular traffic and granular matter [3,4], folding of proteins [5], formation and 
growth of bacterial colonies [6] , biological evolution of species [7] and transac- 
tions in financial markets [8] are just a few examples of exotic phenomena in 
such systems. Most of these systems are interesting not only from the point of 
view of Natural Sciences for fundamental understanding of how Nature works 
but also from the points of view of applied sciences and engineering for the 
potential practical use of the results of the investigations. Our review of the 
current status and future trends of research on the theory of vehicular traffic 
(and some related systems) will, we hope, convince you that, indeed, the re- 
sults of recent studies of complex systems have been "a conceptual revolution, 
a paradigm shift that has far reaching consequences for the very definition of 
physics" [9]. 

For almost half a century physicists have been trying to understand the fun- 
damental principles governing the flow of vehicular traffic using theoretical 
approaches based on statistical physics [3,4,10-13]. The approach of a physi- 
cist is usually quite different from that of a traffic engineer. A physicist would 
like to develop a model of traffic by incorporating only the most essential in- 
gredients which are absolutely necessary to describe the general features of 
typical real traffic. The theoretical analysis and computer simulation of these 
models not only provide deep insight into the properties of the model but also 
help in better understanding of the complex phenomena observed in real traf- 
fic. Our aim in this review is to present a critical survey of the progress made 
so far towards understanding the fundamental aspects of traffic phenomena 
from the perspective of statistical physics. 

There are two different conceptual frameworks for modeling traffic. In the 

"coarse-grained" fluid-dynamical description, the traffic is viewed as a com- 
pressible fluid formed by the vehicles but these individual vehicles do not 
appear explicitly in the theory. In contrast, in the so-called "microscopic" 
models of vehicular traffic attention is explicitly focused on individual vehi- 
cles each of which is represented by a " particle" ; the nature of the interactions 
among these particles is determined by the way the vehicles influence each oth- 
ers' movement. In other words, in the "microscopic" theories vehicular traffic 
is treated as a system of interacting "particles" driven far from equilibrium. 
Thus, vehicular traffic offers the possibility to study various fundamental as- 
pects of the dynamics of truly nonequilibrium systems which are of current 
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interest in statistical physics [14-17]. 

In order to provide a broad perspective, we describe both "macroscopic" and 
"microscopic" approaches although we put more emphasis on the latter. Some- 
times the phenomenological equations of traffic flow in the "macroscopic" 
models can be obtained from microscopic considerations in the same spirit in 
which macroscopic or phenomenological theories of matter are derived from 
their molecular-theoretic description. 

At present, even within the conceptual framework of "microscopic" approach, 
there are several different types of mathematical formulations of the dynamical 
evolution of the system. For example, the probabilistic description of vehicular 
traffic in the kinetic theory is developed by appropriately modifying the kinetic 
theory of gases. On the other hand, a deterministic description of the motion 
of individual vehicles is provided by the so-called car-following theories which 
are based on the basic principles of classical Newtonian dynamics. In contrast, 
the so-called particle-hopping models describe traffic in terms of a stochastic 
dynamics of individual vehicles. We explain the guiding principles behind all 
these formulations. But we discuss in detail the results obtained mainly from 
the investigations of the recently developed "particle-hopping" models which 
are usually formulated using the language of cellular automata (CA) [18]. 
At present, there is no traffic model yet which can account for all aspects of 
vehicular traffic. In this review we consider a wide variety of CA models which 
describe various different types of traffic phenomena. 

We map the particle- hopping models of vehicular traffic onto some other model 
systems; these mappings indicate the possibility of exploiting powerful tech- 
niques, used earlier for other systems, to study traffic models and, sometimes, 
enable us to obtain results for traffic models directly from the known results 
for models of other systems. We present pedagogical summaries of the statis- 
tical mechanical treatments of the CA models of traffic. We critically examine 
the regimes of validity of the approximation schemes of analytical calculations 
which we illustrate with explicit calculations in those limiting cases where 
these usually yield exact results. The results of the theoretical analysis of 
these models are compared with those obtained from computer simulations 
and, wherever possible, with the corresponding empirical results from real 
traffic. We also compare vehicular traffic with many other similar physical 
systems to show the ubiquity of some physical phenomena. 

Computer simulations are known to provide sufficiently accurate quantitative 
results when analytical treatments require approximations which are too crude 
to yield results of comparable accuracy. In this review we demonstrate how 
computer simulations often help in getting deep insight into various phenom- 
ena involved in traffic and in qualitative understanding of the basic principles 
governing them thereby avoiding potentially hazardous experiments with real 
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traffic. Computer sinmlations of the "microscopic" models of traffic have not 
only attracted the attention of a growing number of statistical physicists in 
the recent years, but have also been received positively by many traffic en- 
gineers. The ongoing research eflforts to utihze computer simulations of the 
microscopic models for practical applications in planning and design of trans- 
portation networks have been reviewed very recently by Nagel et al. [19]. 

Our review is complimentary to those published in the recent years by Helbing 
[13] and by Nagel et al. [19]. A large number of important papers on traffic 
published in recent years are based on the particle-hopping models. But these 
works have received very little attention in [13]. We discuss the methods and 
results for the particle-hopping models in great detail in this review after 
explaining the basic principles of all the theoretical approaches. Moreover, we 
focus almost exclusively on the fundamental principles from the point of view 
of statistical physics while Nagel et al. [19] emphasize practical applications 
which are directly relevant for traffic engineering. 

At this point a skeptic may raise a serious question: " can we ever predict traffic 
phenomena with statistical mechanical theories without taking into account 
effects arising from widely different human temperaments and driving habits 
of the individual drivers?" We admit that, unlike the particles in a gas, a driver 
is an intelligent agent who can "think", make individual decisions and "learn" 
from experience. Besides, the action of the driver may also depend on his/her 
physical and mental states (e.g., sorrow, happiness, etc.). It is also true that 
the behavior of each individual driver does not enter explicitly into the " micro- 
scopic" models of traffic. Nevertheless, as we shall show in this review, many 
general features of traffic can be explained in general terms with these mod- 
els provided the different possible behavioural effects are captured collectively 
through a probabilistic description which requires only a few phenomenologi- 
cal parameters. Similar strategies have been suggested also for capturing the 
behavioural effects of some individual traders in financial markets collectively 
through probabilistic descriptions [1]. These probabilistic descriptions make 
the dynamics of the models intrinsically stochastic. 

Throughout this article the terms vehicle (or, car) and driver are used in- 
terchangeably, although each of these terms usually refers to the composite 
unit consisting of the vehicle and the driver. The main questions addressed by 
physicists are posed as problems in Section 2. As a motivation we present some 
relevant empirically observed general features of real traffic as well as their 
plausible phenomenological explanations in Section 3. The following sections 
review the different theoretical approaches. The classification of the models 
into different classes is not unique. Mostly we choose a classification according 
to the use of discrete or continuous space, time and state variables. The con- 
ceptual basis of the older theoretical approaches, namely, the fluid-dynamical 
theories, the kinetic theories and the car-following theories, are explained, in 
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the Sections 4,5 and 6, respectively, where the corresponding recent develop- 
ments are also summarized. These model classes are continuous in space, time 
and state variables. Some coupled map lattice models of traffic are consid- 
ered in Section 7. They are discrete in time. CA models are discrete in space, 
time and state variables. The model suggested by Nagel and Schreckenberg 
(NaSch) [20] is the minimal model of traffic on highways; the theoretical results 
on various aspects of this model are discussed in Section 8 where the nature 
of the spatio-temporal organization of vehicles are also investigated and the 
fundamental question of the (im-)possibihty of any dynamical phase transi- 
tion in the NaSch model is addressed. Various generalizations and extensions 
of the NaSch model (including those for multi-lane traffic) are reviewed in 
Section 9. The occurrence of self-organized criticality in the so-called cruise- 
control limit of the NaSch model is pointed out. It is demonstrated how addi- 
tional "slow-to-start" rules of CA can give rise to metastability, hysteresis and 
phase- separation in the generalized NaSch models, in qualitative agreement 
with empirical observations. In Section 10 the formation and "coarsening" of 
platoons of vehicles are investigated in an appropriate generalization of the 
NaSch model with one type of quenched randomness; your attention is drawn 
to the formal analogy between this phenomenon and the Bose-Einstein con- 
densation. In this section the effects of other kinds of quenched disorder on the 
nature of the steady-states of the NaSch model are also considered. In Section 
11 we present some other CA models of highway traffic which are not directly 
related to the NaSch model. The Biham-Middleton-Levinc (BML) [21] model 
is the earliest CA model of traffic in idealized networks of streets in cities; it 
exhibits a first order phase transition. A critical review of this model is pre- 
sented in Section 12 together with a list of its generalizations which have been 
reported so far. Furthermore a marriage of the NaSch description of traffic and 
the BML model, which has led to the development of a novel model of city 
traffic, is explained. A brief status report of the ongoing efforts to make prac- 
tical use of the theoretical models for traffic engineering is also presented. The 
similarities between various particle-hopping models of traffic and some other 
systems far from equilibrium are pointed out in Section 13 followed by the 
concluding Section 14 where your attention is also drawn towards challenging 
open questions. Several Appendices deal mostly with more technical aspects 
of some important calculations, but are not necessary for an understanding of 
the main text. 



2 Fundamental and practical questions 

The aim of basic research in traffic science is to discover the fundamental 
laws governing traffic systems. The main aim of traffic engineering is on plan- 
ning, design and implementation of transportation network and traffic control 



9 



systems. Statistical physicists have been contributing to traffic science by de- 
veloping models of traffic and drawing general conclusions about the basic 
principles governing traffic phenomena by studying these models using the 
tools of statistical physics. Moreover, using these models, statistical physicists 
have also been calculating several quantities which may find practical applica- 
tions in traffic engineering. Furthermore, several groups of statistical physicists 
are currently also engaged in developing strategies for fast on-line simulation 
and traffic control so as to optimize traffic flow; significant contributions to 
this traditional domain of traffic engineering can reduce the financial burden 
on the governments. 



2.1 Some fundamental questions 

Because of the apparent similarities between the "microscopic" models of traf- 
fic and macroscopic samples of ionic conductors in the presence of external 
electric field, the tools of statistical mechanics seem to be the natural choice 
for studying these models. However, the actual calculation of even the steady- 
state properties of traffic from the "microscopic" models is a highly difficult 
problem because (apart from the human element involved) (a) the vehicles 
interact with each other and (b) the system is driven far from equilibrium, 
although it may attain a nonequilibrium steady-state. 

In principle, the time- independent observable properties of large pieces of mat- 
ter can be calculated within the general framework of equilibrium statistical 
mechanics, pioneered by Maxwell, Boltzmann and Gibbs, provided the sys- 
tem is in thermal equilibrium. Of course, in practice, it may not be possible 
to carry out the calculations without making approximations because of the 
interactions among the constituents of the system. Some time-dependent phe- 
nomena, e.g., fluctuation and relaxation, can also be investigated using the 
Linear Response Theory provided the system is not too far from equilibrium. 
Unfortunately, so far there is no general theoretical formalism for dealing with 
systems far from equilibrium. Moreover, the condition of detailed balance does 
not hold [22] although a condition of pairwise balance [23] holds for some spe- 
cial systems driven far from equilibrium. 

The dynamical phases of systems driven far from equilibrium are counterparts 
of the stable phases of systems in equilibrium. Some of the fundamental ques- 
tions related to the nature of these phases are as follows. 

(i) What are the various dynamical phases of traffic? Does traffic exhibit 
phase-coexistence, phase transition, criticality [24,25] or self-organized crit- 
icality [26,27] and, if so, under which circumstances? 

(ii) What is the nature of fluctuations around the steady-stores of traffic? Anal- 
ogous phenomenon of the fluctuations around stable states in equilibrium is 



10 



by now quite well understood. 

(iii) If the initial state is far from a stationary state of the driven system, how 
does it evolve with time to reach a truly steady-state? Analogous phenomena 
of equilibration of systems evolving from metastable or unstable initial states 
through nucleation (for example, in a supersaturated vapour) or spinodal de- 
composition (for example, in a binary alloy) have also been studied earlier 
extensively [28]. 

(iv) What are the effects of quenched (static or time-independent) disorder on 
the answers of the questions posed in (i)-(iii) above? 



2.2 Some practical questions 



Let is first define some characteristic quantitative features of vehicular traffic. 
The flux J, which is sometimes also called flow or current, is defined as the 
number of vehicles crossing a detector site per unit time [29]. The distance 
from a selected point on the leading vehicle to the same point on the following 
vehicle is defined as the distance-headway [29]. The time-headway is defined 
as the time interval between the departures (or arrivals) of two successive 
vehicles recorded by a detector placed at a fixed position on the highway [29] . 
The distributions of distance-headways and time-headways are regarded as 
important characteristic of traffic flow. For example, larger headways provide 
greater margins of safety whereas higher capacities of the highway require 
smaller headways. 

Let us now pose some questions which are of practical interest in traffic engi- 
neering. 

(a) What is the relation between density c and flux J? In traffic engineering, 
this relation is usually referred to as the fundam,ental diagram. 

(b) What are the distributions of the distance-headway and time-headway? 

(c) How should on- and off-ramps be designed? 

(d) Does an additional lane really lead to an improvement? 

(e) What are the effects of a new road on the performance of the road net- 
work? 

(f) What type of signaling strategy should be adopted to optimize the traffic 
flow on a given network of streets and highways? 

(g) The generahzed travelling salesman problem: Is the shortest trip also the 
fastest? 
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3 Some empirical facts and phenomenological explanations 

For several reasons, it is difficult to obtain very reliable (and reproducible) 
detailed empirical data on real traffic. First of all, unlike controlled experi- 
ments performed in the conventional fields of research in physical sciences, it 
is not possible to perform such laboratory experiments on vehicular traffic. In 
other words, empirical data are to be collected through passive observations 
rather than active experiments. Secondly, unambiguous interpretation of the 
collected data is also often a subtle exercise because traffic states depend on 
several external influences, e.g. the weather conditions. The systematic inves- 
tigation of traffic flow has a quite long history [30-32]. Although we now have 
a clear understanding of many aspects of real traffic several other controversial 
aspects still remain intellectual challenges for traffic scientists. In this section 
we give an overview of some of the well understood experimental findings, 
which arc relevant for our theoretical analysis in the following sections. More- 
over, wherever possible, we provide phenomenological explanations of these 
empirically observed traffic phenomena. Furthermore, we shall also mention 
some of the more recent empirical observations for which, at present, there 
are no generally accepted explanations. 

3.1 Acceleration noise 

In general, because of the different human temperaments and driving habits, 
different drivers react slightly differently to the same conditions on a high- 
way, even when no other vehicle influences its motion. Consequently, even 
on an empty stretch of a highway, a driver can neither maintain a constant 
desired speed nor accelerate in a smooth fashion. In addition to the type of 
the highway (i.e., the surface conditions, frequency of the curves, etc.) the 
driver-to-driver fluctuation of the acceleration also depends on the density of 
vehicles on the highway. The root-mean-square deviation of the acceleration of 
the vehicles is a measure of the so-called acceleration noise. The distributions 
of the accelerations have been measured since mid nineteen fifties and are well 
documented [33,34]. 

3.2 Formation and characterization of traffic jams 

Traffic jam is the most extensively studied traffic phenomenon. Traffic jams 

can emerge because of various different reasons. Most often traffic jams are 
observed at bottlenecks, e.g. lane-reductions or crossings of highways [35]. At 
bottlenecks the capacity of the road is locally reduced thereby leading to the 
formation of jams upstream traffic. Downstream the bottleneck, typically, a 
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free-flow region is observed. In addition, traffic accidents, which also lead to a 
local reduction of the capacity of the highway, can give rise to traffic jams. 




Fig. 1. Experimental flow speed-diagrams near on- and off ramps of a three lane 
highway. The upper part of the figure shows the empirical results. Each dot rep- 
resents a five-minute average of the local measurement. The lines serve merely as 
guide to the eyes. The lower part of the figure shows the location of the detectors 
(from [32]). 

Fig. 1 shows as an example empirical data of the velocity at a three-lane 
highway close to an on- and off-ramp [32]. The data show that downstream 
the bottleneck (at detector C) no slow vehicle has been recorded. In contrast, 
in the merging regime near detector A vehicles often have to move slowly. 
In between the on- and off-ramps the vehicles move with larger velocities 
compared to those in location A although the number of vehicles passing 
detector B is maximal. Therefore, the on-ramp causes a local reduction of the 
capacity of the highway. 

Perhaps, what makes the study of traffic jams so interesting is that jams often 
appear, as if, from nowhere (apparently without obvious reasons) suddenly on 
crowded highways; these so-called "phantom jams"[^ are formed by sponta- 
neous fluctuations in an otherwise streamlined flow. Direct empirical evidence 
for this spontaneous formation of jams was presented by Treiterer [36] by an- 
alyzing a series of aerial photographs of a multi-lane highway. In Fig. 2 the 
picture from [36] is redrawn. Each line represents the trajectory of an indi- 
vidual vehicle on one lane of the highway[^ The space-time plot (i.e., the 
trajectories x{t) of the vehicles) shows the formation and propagation of a 
traffic jam. In the beginning of the analysed time vehicles are well separated 
from each other. Then, due to fluctuations, a dense region appears which, fi- 
nally, leads to the formation of a jam. The jam remains stable for a certain 

^ "Stau aus dem Nichts" in german 

^ The discontinuous trajectories correspond to vehicles changing the lane. 
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period of time but, then, disappears again without any obvious reason. This 
figure clearly establishes not only the spontaneous formation of traffic jam but 
also shows that such jams can propagate upstream (opposite to the direction 
of flow of the vehicles). Moreover, it is possible that two or more jams coexist 
on a highway. 




Fig. 2. Trajectories of individual vehicles on a single lane of a multi-lane high- 
way. The trajectories were drawn from aerial photographs. During the analyzed 
time-interval the spontaneous formation, propagation and dissolution of a jam has 
been observed (from [36]). 

A more detailed analysis of traffic jams in absence of hindrances has been given 
by Kerner and Rehborn [37-39] who pointed out the following characteristic 
features of wide jams. They found that the upstream velocity and, therefore, 
the outflow from a jam is approximately constant. The outflow from a jam and 
the velocity of the jam fronts are now regarded as two important empirical 
parameters of highway traffic which can be used for calibrating theoretical 
models. 



3.3 Flux-density relation 

Obviously, traffic flow phenomena strongly depend on the occupancy of the 
road. What type of variation of flux and average velocity {v) with density c can 
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one expect on the basis of intuitive arguments? So long as c is sufficiently small, 
the average speed {v) is practically independent of c as the vehicles arc too far 
apart to interact mutually. Therefore, at sufficiently low density of vehicles, 
practically "free flow" takes place. However, from the practical experience 
that vehicles have to move slower with increasing density, one expects that at 
intermediate densities. 



dc 



<0, 



when the forward movement of the vehicles is strongly hindred by others 
because of the reduction in the average separation between them. A faster- 
than- linear monotonic decrease of (v) with increasing c can lead to a maximum 
[29] in the flux (J) = (cv) at c = Cm', for c < Cm, increasing c would lead 
to increasing (J) whereas for c > Cm sharp decrease of (v) with increase of 
c would lead to the overall decrease of (J). However, contrary to this naive 
expectation, in recent years some nontrivial variation of flux with density have 
been observed. The nature of the variation of the flux with the density is still 
not clearly understood [40] since the details of the complex experimental setup 
can strongly influence the empirical results. 
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Fig. 3. Empirical data for flow and occupancy. The data have been collected by 
counting loops on a Canadian highway. Both the occupancy and the flow have been 
directly measured by the detector. Each point in the diagram corresponds to an 
average over a time interval of five minutes (from [40]). 

Fig. 3 shows typical time averaged local measurements of the density and flow 
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which have been obtained from the Queen Ehzabeth Way in Ontario (Canada) 
[40] . At low densities the data indicate a hnear dependence of the flow on the 
density. In contrast strong fluctuations of the flow at large densities exist which 
prevents a direct evaluation of the functional form at high densities. 
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Fig. 4. Schematic representation of fundamental diagrams consistent with empirical 
investigations (from [40]). 

In order to use the empirical results for a theoretical analysis it is often more 
convenient to use the mean- values of the flow at a given density. Fig. 4 shows 
a collection of possible forms of averaged fundamental diagrams consistent 
with empirical data [40]. While the discontinuity of the fundamental diagram 
now seems to be well estabhshed [41] no clear answer can be given to the 
question on the form of the diagram in the free-flow or high-density regime. 
In the low density regime linear as well as non-linear functional forms of the 
fundamental diagrams have been suggested. For the high density branch no 
consistent picture for the high density branch exists. Here the results strongly 
depend on the specific road network. 

In several situations it has been observed that J does not depend uniquely 
on c in an intermediate regime of density; it indicates the existence of hys- 
teresis effects and meta-stable states. In the context of traffic flow, hysteresis 
effects have the following meaning: if a measurement starts in the free-flow 
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Fig. 5. Time-traced measurements of the flow (from [40]). Each data point cor- 
responds to a time-averaged value of the flow and density obtained from a local 
measurement. The average was performed over flve minute intervals. 

regime, an increase of the density leads to an increase of the flow. However, 
beyond a certain density, a further increase of the density leads to a discon- 
tinuous reduction of the stationary flow ("capacity drop") and jams emerge. 
The corresponding fundamental diagram (Fig. 4 A, B) then has the so-called 
'inverse- A form'. Figure 5 shows an experimental verification of a hysteresis 
loop [40] at a transition from a free flow to a congested state. 

Recent empirical observations which have been obtained near a crossing of 
highways exhibit [41] a flat plateau (i.e., a density-indepedent flux) over an 
intermediate regime of density of the vehicles. 



3.4 Microscopic states of traffic flow and phase transitions 



The results for the flux-density relation already suggest the existence of at least 
two different dynamical phases of vehicular traffic on highways, namely a free- 
flow phase and a congested phase. In the free-flow regime all vehicles can move 
with high speed close to the speed limit. The nature of the congested traffic is 
still under debate. Careful empirical observations in the recent years indicate 
the existence of two different congested phases, namely, the synchronized phase 
and the stop-and-go traffic phase [38,39,42-44]. Vehicles move rather slowly 
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in the synchronized states, as compared to the free-flow states, but the flux in 
the synchronized states can take a values close to the optimum value because 
of relatively small headways. Besides, non-trivial strong correlations between 
the density on different lanes exist in the synchronized state [38,44,45] which 
actually motivates the notation synchronized traffic. The stop-and-go traffic 
differs from the synchronized states in the sense that every vehicle inside the 
jams come to a complete halt for a certain period of time. 




Fig. 6. Fundamental diagrams corresponding to different types of synchronized flow. 
Measurements which correspond to synchronized traffic of type (i) are locaHsed in 
a narrow region of the flow density plane, while synchronized states of type (ii) are 
similar to the results of free-flow states. In contrast for synchronized traflfic of type 
(iii) the time-traced measurements show an unsystematic behaviour. 

Following Kerner three different types of synchronized traffic can be dis- 
tringuished by the time-dependent behaviour of the density and flow [38,39,44] . 
In synchronized traffic of type (i) constant values of density and flow can be 
observed during a long period of time. In synchronized traffic of type (ii) pat- 
terns of density and flow quite similar to free flow states have been observed. 
The differences between synchronized states of type (ii) and free-flow are given 
by the reduced average velocities and the alignment of the speeds on different 
lanes in synchronized traffic. Moreover irregular patterns of time-traced mea- 
surements of the flow have been found in synchronized traffic (see Fig. 6) of 
type (in). 

3.5 Time- and distance-headways 

The flux J can be written as J — N/T where T — J^^ik is the sum of 
the time-headways recorded for all the A'" vehicles. Hence, J — 1/Tav where 
Tav = (1/^) Si^i is the average TH. Therefore, the TH distribution contains 
more detailed informations on traffic flow than that available from the flux 
alone. With the variation of density c of the vehicles, T^v exhibits a minimum 
at c = where the flux is maximum [29] . 
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Fig. 7. Time-headway distribution for different density regimes. Top: In free-flow 
traffic tfie At-distribution is dominated by two peaks at 0.8 sec and 1.8 sec. Middle: 

In synchronized traffic cars with narrow time gaps are to find as well as a dominating 
peak. Bottom: In stop-and-go traffic short time headways are surpressed. The peak 
at 1.8 sec remains since vehicles are leaving the jam with a typical temporal headway 
of approximately 2 sec. 

The results discussed in the preceeding subsections are based on time averaged 
local measurements. But it is also very useful to analyze the single-vehicle 
data directly [46-48,41]. The single- vehicle data allows calculation of the time- 
headway distributions [41]. All the timc-hcadway distributions in the free-flow 
regime show a two peak structure. The first peak at At = 0.8 sec represents the 
global maximum of the distribution. On a microscopic level these short time- 
headways correspond to platoons of some vehicles traveling very fast - their 
drivers are taking the risk of driving " bumper-to- bumper" with a rather high 
speed. These platoons are the reason for the occurrence of high-flow states in 
free traffic. The corresponding states exhibit metastability, i.e. a perturbation 
of finite magnitude and duration is able to destroy such a high-fiow state [39]. 
Additionally, a second peak emerges at At — 1.8 sec which can be associated 
with a typical drivers' urge to maintain a temporal headway of ~ 2 sec (which 
is the safe distance recommended in driving schools). 

Surprisingly, the small time headways have much less weight in congested 
traffic. Only the peak at At — 1.8 sec is recovered, where the time headway 
distribution now takes the maximum value. But nevertheless, the small time 
headways (At < 1.8 sec) contribute significantly in synchronized traffic. In 
stop-and-go traffic only the 1.8 sec-peak remains and short time-headways are 
surpressed. The asymptotic behavior is rather unsystematic and refiects the 
dynamics of vehicles inside the jams. 
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Fig. 8. Speed-distance relation in free and congested traffic. It turns out that the 
typical velocities at large distance headways strongly depend on the present traffic 
state (from [41]). 

Another important result characterizing the microscopic states is the depen- 
dence of the velocity of individual vehicles on the distance headway (see Fig. 8). 
This function is also of great importance for theoretical approaches, e.g. it is 
used as input for the so-called optimal velocity model. In the free-flow regime 
it is evident that the asymptotic velocity is reached already for small distance- 
headways. The slope of the velocity function is much lower than that in the 
free-flow regime. Surprisingly, the asymptotic velocity depends strongly on 
the microscopic state as well as on the density, e.g., in dense traffic low veloc- 
ities of the vehicles are also observed even when large distance headways are 
available. 



4 Fluid-dynamical Theories of vehiculcir trcifRc 

When viewed from a long distance, say, an aircraft, flow of fairly heavy traffic 
appears like a stream of a fluid. Therefore, a "macroscopic" theory of traffic can 
be developed, in analogy with the hydrodynamic theory of fluids, by treating 
traffic as an effectively one-dimensional compressible fluid (a continuum) [49] . 
We follow the convention that traffic is flowing from left to right. 

Suppose, c{x; t) and J{x; t) are the "coarse-grained" density and flux at an ar- 
bitrary location x at an arbitrary instant of time t. The equation of continuity 
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for the fluid representing traffic is 

dc(x; t) dJ(x; t) , 

dt + = § ~ " ?i ' " 

where the first and the second terms on the right hand side take care of 

the sources and sinks, respectively, at the Jj„ on-ramps situated at Xi (i = 
1, 2, Jin) and Jout off-ramps situated at Xj (j = 1, 2, Jout)- We can write 
ai{x — xf,t) and Pj{x — Xj] t) as 

ai{x — Xi, t) — ai{t)(f)i{x — Xi) and Pj{x — Xj]t) — Pj{t)(f)j{x — Xj) (3) 

where — Xi) and — Xj) describe the spatial distribution of the in- 
coming and outgoing ffux, respectively, while a^{t) and Pj{t) account for the 
corresponding temporal variations. 

In the following subsections, for simpUcity, we shall consider a given stretch 
of highway with no entries or exits. In such special situations the equation of 
continuity reduces to the simpler form [50] 

t) dJ{x;t) ^ 
dt dx ' ^ ^ 

One cannot get two unknowns, namely, c{x; t) and J{x; t) (or, equivalently, 
v{x;t)) by solving only one equation, namely (4), unless they arc related to 
each other. In order to proceed further, one needs another independent equa- 
tion, say, for v{x;t); we shall write down such an equation later in subsection 
4.4. An alternative possibihty, which Lighthill and Whitham [50] adopted in 
their pioneering work, is to assume that J{x; t) is determined primarily by 
the local density c{x; t) so that J{x; t) can be treated as a function of only 
c(x;t). Consequently, the number of unknown variables is reduced to one as, 
according to this assumption, the two unknowns c(a;; t) and J{x] t) are not 
independent of each other. 

^.i Lighthill- Whitham theory and kinematic waves 

As a ffrst approximation, let us begin with Lighthill- Whitham assumption 
that 

Jix;t) ^j{c{x;t)) (5) 

where j(c) is a function of c. The functional relation (5) between density and 
ffux cannot be calculated within the framework of the ffuid-dynamical the- 
ory; this must be either taken as a phenomenological relation extracted from 
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empirical data or derived from more microscopic considerations. In general, 
the flux-density curve implied by equation (5) need not be identical with the 
fundamental diagram in the steady-state. 

Under the assumption (5), the x-dependence of the local fiux J{x\ t) arises only 
from the x-dependence of c{x\t). Alternatively, since J{x;t) = c{x;t)v{x;t), 
assuming v{x;t) to depend only on c{x;t) the x-dependence of v{x;t) arises 
only from the x-dependence of c{x; t). Using (5) the equation of continuity (4) 
can be expressed as 



dc{x; t) ^ dc{x; t) 
dt dx 



v[x;t) + c{x;t) — 
ac 



^ + „,^ = (6) 

dt dx 



where Vg — dJ/dc. Note that the equations (5) and (6) form the complete sys- 
tem of dynamical equations governing traffic flow in this ffist approximation. 
However, the equation (6) is non-/mear because, in general, Vg = dJ/dc = 
v{c) + c^^^ depends on c. If Vg were a constant wq, independent of c, equa- 
tion (6) would become linear and the general solution would be of the form 
c{x;t) = f{x — vot) where / is an arbitrary function of its argument. In that 
case, the solution of any particular problem would be found by merely match- 
ing the function / to the corresponding given initial and boundary conditions. 
Such a solution describes a density wave motion as an initial density profile 
would get translated by a distance vot in a time interval t without any change 
in its shape. However, the non-linearity of the equation (6) gives rise to sub- 
tleties which are essential to capture at least some aspects of real traffic. 

The solution of the nonlinear equation (6) is of the general form 

c{x;t) ^ F{x-Vgt), (7) 

where F is an arbitrary function of its arguments. If we define a wave to be 
"recognizable signal that is transferred from one part of a medium to another 
with a recognizable velocity of propagation" [51] then the solutions of the 
form (7) can be regarded as a "density wave". There arc several similarities 
between the density wave and the more commonly encountered waves like, 
for example, acoustic or elastic waves. But, the acoustic or elastic waves are 
solutions of linearized partial differential equations whereas the equation (6) 
is nonlinear and, hence, Vg is c-dependent. Besides, the waves of the type (7) 
are called kinematic waves [50-52] to emphasize their purely kinematic origin, 
in contrast to the dynamic origin of the accoustic and elastic waves. From 
the initial given density profile c{x; 0) the proffie c{x; At) at time At can be 
obtained by moving each point on the initial proffie a distance Vg{c)t to the 
right; obviously, the distance moved is different for different values of c. The 
time-evolution of the density profile can be shown graphically [50,51] on the 
space-time diagram (i.e., the x — t plane) where an arbitrary point xq on the 
t — axis moves along a straight line of slope Vg{c) if the initial density at Xq is 
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c. These straight hncs arc referred to as characteristics; different characteristics 
corresponding to different c have different slopes Vg{c). 



The speed Vg{c) of the density wave should not be confused with v{c), the 
actual speed of the continuum fluid representing traffic. In fact, at any instant 
of time v{x; t) can be obtained from the corresponding density profile c{x; t) by 
using the relation v{x\ t) = j{c(x; t))/c(x; t). Moreover, since Vg = v{c) + c^^ 
and since ^^^^^ < 0, the speed of the density wave is less than that of the 
fluid. Therefore, the density wave propagates backward relative to the traffic 
and the drivers are thereby warned of density fluctuations ahead downstream. 
Furthermore, the density wave moves forward or backward relative to the road, 
depending on whether c < Cm or c> Cm where Cm corresponds to the maximum 
in the function j{c). 

When J(c) is convex, i.e., (PJ/dc^ < 0, we have dvg/dc < 0; consequently, 
higher values of c propagate slower than lower values of c thereby distorting 
the initial density profile. On the other hand, when dVg/dc > higher values of 
c propagate faster and the distortion has the opposite tendency as compared 
to the case of dvg/dc < 0. In both the situations the distortion of the initial 
density profile is caused by the c-dependence of Vg which arises from the non- 
linearity of the equation (6). The distortion of the density profile with time 
can also be followed on the space-time diagram. If dvg/dc < 0, in regions of de- 
creasing density (i.e., c{xi) > c{x2) for xi < X2) the characteristics move away 
from each other whereas, in regions of increasing density, the characteristics 
move towards each other. 

When two characteristic lines on the space-time diagram intersect the density 
would be double-valued at the point of intersection. We can avoid this appar- 
ently impossible scenario by the following interpretation: When two character- 
istic lines intersect a shock wave is generated. By definition, a shock represents 
a mathematical discontinuity in c and, hence, also in v. The speed of a shock 
wave is given by 

^^^Jjq^ (8) 

where c"*" and c~ are, respectively, the densities immediately in front (down- 
stream) and behind (upstream) the shock while J(c+) and J(c~) represent the 
corresponding downstream and upstream fluxes, respectively. Note that the 
shock wave moves downstream (upstream) if Vg is positive (negative). Often 
the shock is weak in the sense that the relative discontinuity (c^ — c~)/c~ is 
small and in such cases the shock wave speed tends to Vg — dJ/dc. As a shock 
separates a section of high an low densities of the model, it corresponds to a 
section of a highway where a free- flow and a congested regime is present. In 
particular for large differences between and c~ the velocity of the shock 
can be interpreted as the velocity of a backwards moving jam. 
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One advantage of the kinematic approach outhned above over any dynamic 
approach is that the dynamical equation, which will be given in Sec. 4.4, is 
difficult to derive from basic first principles and usually involve quite a few 
phenomenological parameters and even a phenomenological function. On the 
other hand, the only input needed for the kinematic approach is the phe- 
nomenological function J(c) which can be obtained from empirical data. 



4-2 Diffusion term in Lighthill-Whitham theory and its effects 



Let us now make improvement over the original Lighthill-Whitham theory, 
which is based on the first approximation (5). We now assume that the local 
fiux J{x;t) is determined not only by the local density c{x; t) but also by the 
gradient of the density. In other words, we replace the assumption (5) by 

J{c)=j{c)-D^ (9) 

where D is a positive constant. Note that, for fixed c{x;t) (and, hence, fixed 
j(c)), a positive (negative) density gradient leads to a lower (higher) flux as the 
drivers are expected to reduce (increase) the speed of their vehicles depending 
on whether approaching a more (less) congested region. Using the relation (9) 
in the equation of continuity (4) we now get 

dc{x;t) , dc{x;t) _ ^ d^c{x;t) 

dt ^ ' dx ~ dx^ ^ ' 

where Vg{c) = The equation (10) reduces to the equation (6) when D — 0. 
The nonhnearity and diffusion have opposite effects: the term Vg{c)^ leads to 
"steepening" and ultimate "breaking" of the wave whereas the term Dd^c/dx^ 
smoothens out the profile. Nonvanishing D also leads to a non-zero width of 
the shock wave. 



4-3 Greenshields model and Burgers equation 



So far in the preceding subsections we have not considered any specific form 
of the function j(c) relating fiux with density. One can start with the simplest 
(differentiable) approximation capturing the basic form of the fundamental 
diagram, 

J ^VmaxC{l- C). (11) 

Note that Vmax in (H) is a phenomenological parameter and it is interpreted 
to be the maximum average speed for c — > 0. In traffic science and engineering. 
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one usually uses 1 — c/cjam instead of 1 — c in the equation (11) and the cor- 
responding form of the relation between J and c is known as the Greenshields 
[30] model. Substituting (11) into the equation (10) we get 

dc{x-t) dc{x;t) d\ 

at OX OX OX^ 

Introducing the linear transformation of variables 

~ "^maxt ) t — t (13) 

one gets 

dt' dx' ~ dx" ^ ' 

which is the (deterministic) Burgers equation [51,53]. Note that the trans- 
formation (13) takes one from the space-fixed coordinate system {x,t) to a 
coordinate system {x',t') that moves with uniform speed Vmax', so, vehicles 
moving with speed v^ax with respect to the coordinate system {x, t) do not 
move at all with respect to the coordinate system {x',t'). 

The advantage of this route to the theory of traffic flow is that the Burgers 
equation (14) can be transformed further into a diffusion equation, thereby 
getting rid of the nonlinearity, through a nonlinear transformation called the 
Cole-Hopf transformation [51]. Since it is straightforward to write down the 
formal solution to the diffusion equation, one can see clearly the role of the 
coefficient D and the nature of the solutions in the limit — > 0. 

If equation (6) is assumed to be the only equation governing traffic flow then 
an inhomogeneous initial state can lead to a shock wave but the amplitude 
of the shock wave decreases with time and eventually the shock wave fades 
out leading to a homogeneous steady state in the limit t — * oo. Leibig [54] 
has studied how a random initial distribution of steps in the density proflle 
evolves with time in this theory. No traffic jam forms spontaneously from a 
state of uniform density at this level of sophistication of the fluid-dynamical 
approach. 



4-4 Navier-Stokes-like momentum equation and consequences 



Corresponding to the assumption (9) we can write a velocity equation 

v{x;t)^v{c{x;t))--^^ (15) 

c ox 
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where v{c) — In the kinematic approach discussed so far in the preced- 

ing subsections it is imphcitly assumed that, following any change in the local 
density (and density gradient) is followed by an immediate response (without 
delay) of the velocity field. For a more realistic description, the local speed 
should be allowed to relax after a non-zero delay time td- So, it seems natural 
to treat the right hand side in (15) as a desired local velocity at x and write 
the total derivative dv/dt of v with respect to time as [55-58] 

dv dv c 
dt dx td 

where the function fsa/e(c) is identical to v{c). Note that fsa/e(c) is a mono- 
tonically decreasing function of c, i.e, dvgafe/dc < 0. Equation (16) is an 
additional dynamical equation describing the time-dependence of the velocity 
v{x; t). 

Now let us interpret the two terms on the right hand side of (16). The phe- 
nomenological function Vsafe{c) gives the safe speed, corresponding to the ve- 
hicle density c, that is achieved in time-independent and homogeneous traffic 
flow and td is the corresponding average relaxation time. Next, note that the 
term D^^^^^ takes into account the natural tendency of the drivers to acceler- 
ate (decelerate) if the density gradient is negative (positive), i.e. if the density 
in front becomes smaller (larger); therefore, it can be interpreted as propor- 
tional to the pressure gradient in the fluid describing traffic. In addition to 
these terms, another term proportional to |^ is also added to the right hand 
side of the velocity equation; this tends to reduce spatial inhomogeneities of 
the velocity field and is usually interpreted as the analogue of the viscous 
dissipation term in the Navier-Stokes equation. 

Thus, finally, in the fluid- dynamical approach, a complete mathematical de- 
scription of the vehicular traffic on highways is provided by two equations, 
namely, the equation of continuity (2) and the Navier-Stokes-like velocity 
equation [59-62] 



./e(c) - V 



-D 



dc{x; t) 
dx 



(16) 



dv 
dt 



dv 
dx 



-D 



dc 



d'^v 



dx dx 



td 



[Vsafe{c)-V] 



(17) 



where D, /i and td are phenomenological constants. 



4-5 Fluid- dynamical theories for multi-lane highways and city traffic 



One can describe the traffic on two-lane highways [63] by two equations each 
of the same form (2) and where the source term in the equation for lane 1 (lane 
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2) takes into account the vehicles which enter into it from the lane 2 (lane 1) 
while the sink term takes into account those vehicles entering the lane 2 (lane 
1) from the lane 1 (lane 2). 

A lattice hydrodynamic theory for city traffic has been formulated recently 
[64] . This fiuid-dynamical model is motivated by the CA model, developed by 
Biham ct al. [21], which will be discussed in detail later in this review. Instead 
of generalizing the Navier-Stokes equation (17) a simpler form of the velocity 
equation has been assumed. 

4-6 Some recent results of the fluid- dynamical theories and their physical 
implications 

The fiuid-dynamical model of vehicular traffic has been studied numerically 
by discretizing the partial differential equations (2) and (17) together with 
appropriate initial and boundary conditions. Both periodic boundary con- 
ditions and open boundary conditions with time-independent external flux 
af{t) — (3f{t) — 7 have been considered. 




Fig. 9. Time evolution of a density profile in the Kiihne-Kerner-Konhauser model 
[57-59] consisting of equations (2) and (17) with a special choice of fsa/e(c)- Initially 
an approximately constant density profile is assumed. The figure shows that the 
homogeneous state is not stable, already for small deviations from a flat profile 
(from [13]). 

In the fluid-dynamical theory, based on the equation of continuity (2) and the 
Navier-Stokes-like equation (17) traffic jams can appear spontaneously, even 
if the initial density profile c{x] 0) deviates very little from the homogeneous 
state (fig. 9). In order to understand the physical mechanism of the forma- 
tion of local cluster of vehicles, let us consider a local increase of the density 
Ac{x) > at some location x. Since dvgafe/dc < 0, the local increase of the 
density leads to a decrease of Vgafe- This decrease in the safe velocity forces 
drivers to reduce their average velocity v sharply if \dvsafe/dc\ is large enough. 
On the other hand, it follows from the equation of continuity that the local 
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decrease of v gives rise to further increase of c around x and, consequently, 
further subsequent decrease of v{x) in this location. This avalanche- like pro- 
cess, which tends to increase the amplitude of the local fluctuation of the 
density around the homogeneous state, competes against other processes, like 
diffusion and viscous dissipation, which tend to decrease inhomogeneities. 

To our knowledge, the first attempt to understand the physical mechanism 
of synchronized traffic within the framework of the fluid-dynamical formalism 
was made by Lee et al. [62]|3J On a finite stretch of highway, of length L, they 
installed an on-ramp and an off-ramp on the model highway with a separation 
of L/2 between them. They chose the spatial distribution of the external flux 
(f){x) in equation (3) as 

(f){x) = {2naY^/^ exp{-xy2a^) (18) 

with a = 56.7m. They also assumed the form fsa/e(c) = ^0(1 — c/c)/[l + 
E{c/c)^] for the safe velocity with adjustable parameters Vq, E, 9, c. Lee et al. 
[62] first allowed the system to reach the steady state after applying a weak 
time-independent flux = = 7 and simulating the time-evolution 

of the traffic by solving simultaneously the equation of continuity (2) and 
the Navier-Stokes-like equation (17) with a specific set of chosen values for 
the parameters td, D, fi,vo, E,9,c, etc. Since they chose 7 < 7c, an initially 
homogeneous traffic reaches a steady "free- flow" where homogeneous regions 
with different densities are separated from each other by narrow "transition 
layers" near the ramps ( no stable "free- flow" exists if 7 > 7c). Then, they 
applied a pulse of additional flux 6q at the on-ramp for a short duration 6t. 
After a transient period, which depends on the parameters of the model, the 
system was found to settle in a limit cycle in which the local density and local 
flux oscillate periodically and the oscillations are localized near the on-ramp. 
The discontinuous change of the spatio-temporally averaged velocity induced 
by the localized perturbations of finite amplitude, associated hysteresis effects 
and the stability of the limit cycle were found to be qualitatively similar to 
some of the empirically observed characteristics of synchronized flow in real 
traffic; therefore, Lee et al. [62] identified the limit cycle observed in their 
theoretical investigation as the synchronized state of vehicular traffic. They 
drew analogy between this state and a "self-excited oscillator" [62]. However, 
this mechanism of the synchronized state is not yet accepted as the true and 
only possible explanation of the phenomena associated with the synchronized 
state observed empirically. 

Meanwhile similar results have been obtained for a gas-kinetic based traffic 
model [66], also using on- and off-ramps in order to explain the transition from 
free-flow to synchronised states. This work was completed in a recent paper. 



A brief summary of that work can be found in [65] . 
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where a phase diagram was calculated, which dependends on the on-ramp ac- 
tivity and the flow on the highway [67]. Summarising the recent results of the 
macroscopic traffic models, there seems to be evidence that on-and off-ramps 
play an important role for a theoretical explanation of synchronised traffic. 
Nevertheless some experimental features are still not captured by these ap- 
proaches. E.g. the empirical results show that for synchonised traffic of type 
(iii) no correlations between density and flow exist [41], in contrast to the reg- 
ular patterns of the oscillating states found in simulations of the macroscopic 
models. 

Despite its success in capturing many aspects of traffic flow the fluid- dynamical 
approach has its limitations; for example, viscosity of traffic is not a directly 
measurable quantity. Nevertheless, the fluid dynamical approach is being pur- 
sued not only by some physicists but also by several members of the traffic 
engineering community [68,69]. 



5 Kinetic theories of vehiculeir trcifRc 



In the kinetic theory, traffic is treated as a gas of interacting particles where 
each particle represents a vehicle. The various different versions of the kinetic 
theory of vehicular traffic [12,70-76] have been developed by modifying the 
kinetic theory of gases. 

Recall that in the kinetic theory of gases [77] /(r,p; t)d^rd^p denotes the num- 
ber of molecules which, at time t, have positions lying within a volume ele- 
ment d^r about r and momenta lying within the momentum-space element 
d^p about p. The Boltzmann equation, which describes the time-evolution of 
the distribution f{x, v; t), is given by 



dt m 




/(r,p;t)=he (19) 



where the symbols and Vp denote gradient operators with respect to r and 
p, respectively, while F is the external force. The term [%] represents the 

\ "^^^ / coll 

rate of change of /, with time, which is caused by the mutual collisions of the 
molecules. 

In the first of the following two subsections we present the earliest version of 
the kinetic theory of vehicular traffic which was introduced by Prigogine and 
coworkers [70,12] by modifying some of the key concepts in the kinetic theory of 
gases and by writing down an equation analogous to the Boltzmann equation 
(19). In the subsequent subsection we discuss the kinetic theory developed 
later by Paveri-Fontana [71] to cure the defects from which the Prigogine 
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theory was found to suffer. 



5.1 Prigogine model 



Suppose f{x,v;t)dxdv denotes the number of vehicles, at time t, between x 
and X + dx, having actual velocity between v and v + dv. In addition, Pri- 
gogine and coworkers [70,12] introduced a desired distribution fdes{x,v) which 
is a mathematical idealization of the goals that the population of the drivers 
collectively strives to achieve. The actual distribution may deviate from the de- 
sired distribution because of various possible influences, e.g., road conditions, 
weather conditions or interaction with other vehicles, etc. They also argued 
that some of these influences cease after some time while the interactions with 
the other vehicles persist for ever. For example, only a short stretch of the 
road surface may be icy and strong winds or rain may stop after a short dura- 
tion; in such situations / can relax to fdes over a relaxation time Trei provided 
mutual interactions of the vehicles is negligibly small. On the basis of these 
arguments, Prigogine and coworkers [70,12] suggested that the analogue of the 
Boltzmann equation for the traffic should have the form 



dt dx 

\ / ret ^ f viii, 

where {^qi) ^ accounts for the relaxation of / towards /^es in the absence of 
mutual interactions of the vehicles while ( ^ ) accounts for the changes of 

V <9M int ^ 

f arising from mutual interactions among the vehicles. Note that the term 
(^^) . ^ on the right hand side of (20) may be interpreted as the analogue of 

the term (^^^ in the equation (19) whereas the term {j^^ ^ in equation 

(20) may be interpreted as the counterpart of the term F ■ Vpf{r,P',t) in the 
equation (19). 

Prigogine and coworkers wrote down an explicit form for the term ^ by 

generalizing that for the term (%] in the kinetic theory of gases. We shall 

V / coll 

consider this term in the next subsection. In order to write down a simple 
explicit form of the relaxation term in the equation (20) they assumed that 
(i) the collective relaxation, which would cause the actual distribution to tend 
towards the desired distribution, involves only a single relaxation time Trei so 
that 

/ fdes (21) 




and 

(ii) the desired speed distribution F^esiv) remains independent of the local 
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concentration c{x; t) so that 

fdes{x, V] t) = c{x; t)Fdes{v) (22) 

Therefore, a more exphcit form of the Boltzmann-hke equation (20) in the 
Prigogine theory is given by 

df , ^df _ /(x.r:f)-c(.r:f)i^..(tO ^ /df\ 

dt dx Trel K^Jint 

Note that, in the absence of mutual interactions of the vehicles, the distri- 
bution f{x,v;t) would relax exponentially with time. The concept of desired 
distribution fdes{x,v;t) and this scenario of collective relaxation of / towards 
fdes has subsequently come under severe criticism [71]. Analyzing a set of 
"ideal experiments" in the light of the Prigogine theory, Paveri-Fontana [71] 
showed that the results obtained from the Boltzmann-like equation (23) are 
physically unsatisfactory. 

More recently, Lehmann [72] has attempted to revive the Prigogine approach 
by reformulating it as a semi- phenomenological theory where the distribution 
f{x, v;t) is assumed to follow the simpler form 

dt dx Trel 

and the effects of the interactions are taken into account implicitly through a 
density-dependent desired distribution function fdes{v,c) which has to deter- 
mined empirically. 



5.2 Paveri-Fontana model 



In order to remove the conceptual as well as mathematical drawbacks of the 
Prigogine model of the kinetic theory of vehicular traffic, Paveri-Fontana [71] 
argued that each vehicle, in contrast to the molecules in a gas, has a desired 
velocity towards which its actual velocity tends to "relax" in the absence of 
"interaction" with other vehicles. Thus, Paveri-Fontana model is based on a 
scenario of relaxation of the velocities of the individual vehicles rather than a 
collective relaxation of the distribution of the velocities. 

In mathematical language, Paveri-Fontana introduced an additional phase- 
space coordinate, namely, the desired velocity. Suppose, g{x, v, Vdes', t)dxdvd,Vdes 
denotes the number of vehicles at time t between x and x + dx, having actual 
velocity between v and v-\-dv and desired velocity between va^s and Vdes + dvdes- 
The one- vehicle actual velocity distribution function 
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f{x,V;t) = / dVdes9{x,V,Vdes;t) 



(25) 



describes the probability of finding a vehicle between x and x + dx having 
actual velocity between v and v + dv at time t. Similarly, the one-vehicle 
desired velocity distribution function 



fo{x, Vdes] t) = j dvg{x, V, Vdes] t) (26) 

describes the probability of finding a vehicle between x and x + dx having 
desirec? velocity between Vdes and Vdes + dvdes- The local density of the vehicles 
c{x; t) at the position x at time t can be obtained from 

roo POO 

c{x;t)= dvdes dvg{x,v,Vdes;t). 
Jo Jo 

Similarly, the corresponding average actual speed {v{x;t)) and the average 
desired speed {vdesix',t)) are defined as 



{vdes{x;t)) 



c{x; t) 

dVdes /o°° dvVdes9{x, V, Vdes] t) 



c{x; t) 

Finally, the local flux J{x;t) is defined as J{x;t) = c{x;t){v{x;t)). 

Now let us assume that the desired velocity of each individual driver is in- 
dependent of time, i.e., dvdes/dt = 0. Of course, the drivers may also adapt 
to the changing traffic environment and their desired velocities may change 
accordingly. In principle, these features can be incorporated into the kinetic 
theory at the cost of increasing complexity of the formalism. 

Next, let us also assume that, in the absence of interaction with other vehicles, 
an arbitrary vehicle reaches the desired velocity exponentially with time, i.e., 
dv/dt— {vdes where r is a relaxation time. The Boltzmann-like kinetic 

equation for g{x, v, Vdes', t) can be written as 



— + v— 
dt dx 



9 + 



d_ 
dv 



Vde 




(27) 



In order to write down an explicit form of the "interaction term" we have to 
model the interactions among the vehicles. First of all, we model the vehicles 
as point-like objects. We consider the scenario where a fast vehicle, when hin- 
dered by a slow leading vehicle, either passes or slows down to the velocity of 
the lead vehicle. Let us now make some further simplifying assumptions: 
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(i) The slowing down takes place with a probability 1 — Ppass where Ppass is 
the probability of passing. 

(ii) If the fast vehicle passes the slower leading vehicle, its own velocity re- 
mains unchanged. 

(iii) The velocity of the slower leading vehicle remains unchanged, irrespective 
of whether the faster following vehicle passes or slows down. 

(iv) The slowing down process is instantaneous, i.e., the braking time is neg- 
ligibly small. 

(v) It is adequate to consider only two- vehicle interactions; there is no need 
to consider three- vehicle (or multi- vehicle) interactions. 

(vi) The postulate of "vehicular chaos", which is the analogue of the postulate 
of "molecular chaos" in the kinetic theory of gases, holds, so that the two- 
vehicle distribution function g2{x, v, Vdes, x', v', v'^^^; t) can be approximated as 
a product of two one-particle distributions g^x^v^Vdes')^) and g{x' ,v' ,v'^gg;t), 
i.e., g2{x,v,Vdes:X',v',v'a^^;t) ~ g{x,v,Vdes;t)g{x',v',v'^^/,t). Thus, the equa- 
tion (27) can be written explicitly as 



d d 
dt dx 



d 



Vdes - V 

■9 



T 

POO 



poo 

= f{x,V]t) / dv'{l - Ppass)W - '")9{x,v\Vdes]t) 
Jv 

rv 

-g(x,v,Vdes;t) dv' {1- Ppass) {v-v')f{x,v';t) (28) 

J 

where the form of the "interaction term" on the right-hand side of the equa- 
tion (28) follows from the assumptions (i)-(vi) above. The first term on the 
right hand side of (28) describes the "gain" of probability g{x,v,Vdes',i) from 
the interaction of vehicles of actual velocity v' with slower leading vehicle of 
actual velocity v while the second term describes the loss of the probability 
g{x,v, Vdes] t) arising from the interaction of vehicles of actual velocity v with 
even slower leading vehicle of actual velocity v'. 

The stationary homogeneous solution g{v,Vdes) is, by definition, independent 
of X and t. But, to our knowledge, so far it has not been possible to get even 
this solution of the Boltzmann-like integro-differential equation (28) by solving 
it analytically even for the simplest possible choice of the desired distribution 
function although numerical solutions [74,76] provide some insights into the 
regimes of validity of the equation (28) and gives indications as to the di- 
rections of further improvements of the Paveri-Fontana model. For example, 
the finite sizes of the vehicles must be taken into account at high densities. 
Besides, the assumption (iv) of instantaneous relaxation has also been relaxed 
in a more recent extension [74] . 

Normally passing would require more than one lane on the highway. There- 
fore, the models discussed so far in the context of the kinetic theory may be 
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regarded, more appropriately, as quasi- one- dimensional. These neither deal ex- 
plicitly with gi{x, V, Vdes] t) ^or the individual lanes (labeled by i) nor take into 
account the process of lane- changing. Besides, all the vehicles were assumed 
to be of the same type. Now, in principle, we can generalize the formalism of 
the kinetic theory of traffic to deal with different types of vehicles on multi- 
lane highways. Suppose gf{x,v,Vdes',t) is the distribution for vehicles of type 
a on the z-th lane of the highway. Obviously, the Boltzmann-like equations for 
the different lanes are coupled to each other. However, one needs additional 
postulates to model the lane-changing rules [78,79]. 

Very little work has been done so far on developing kinetic theories of two- 
dimensional traffic flow which would represent, for example, traffic in cities. 
Suppose, for simplicity, that the network of the streets consists of north-south 
and east-west streets and that east-west streets allow only east-bound traffic 
while only north-bound traffic flow takes place along the north-south streets. 
Let Pxx (Pxy) denote the probabilities of an east-bound vehicle passing another 
east-bound (north-bound) vehicle. Similarly, suppose, Py^ {Pyy) denote the 
probabilities of a north-bound vehicle passing an east-bound (north-bound) 
vehicle. Let gx{x,y,u,Ufies]t) and gy{x,y,v,Vfies]t) represent the distributions 
for the east-bound and north-bound vehicles, respectively, where u and v refer 
to the actual their actual velocities whereas Udes and Vdes refer to the corre- 
sponding desired velocities. The Boltzmann-like equations governing the time 
evolutions of these distributions are given by [80] 



d d d 
dt dx dy 
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du 
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dt 
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where 



'dgx 
dt 



x.coll 



roo 

-■fx{x,y,u;t) / du'{l - Pxx){u' - u)gx{x,y,u',Udes;t) 

Ju 

ru 

■gx(x,y,u,Udes]t) du'{l- Pxx){u-u')fa:{x,y,u';t) (30) 
Jo 



and 



'dgx 
dt 



y,coll 



roo roo 

■ K,o dv'fy{x, y, v'; t) du'{l - Pxy)u'gx{x, y, u', Udes] t) 

roo 

gx{x,y,u,Udes]t) dv\l- Pxy)ufy{x,y,v']t) (31) 

^0 



The first term on the right hand side of the equation (30) describes gain of 
population of east-bound vehicles with velocity u because of interaction with 
other east-bound vehicles with velocity u' >u while the second term describes 
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the loss of population of east-bound vehicles because of interaction with east- 
bound vehicles of velocity u' < u. The right hand side of the equation (31) 
is based on the assumption that when an east-bound vehicle interacts with a 
north-bound vehicle at a crossing, it either passes or stops. 

5.3 Derivation of the phenomenological equations of the macroscopic fluid- 
dynamical theories from the microscopic gas-kinetic models 

In this section we discuss the results of the attempts to derive the phenomeno- 
logical equations of traffic flow in the macroscopic fluid-dynamical theories 
from the microscopic gas-kinetic models. Several attempts have been made so 
far to derive the equation of continuity and the Navier-Stokes-like equation 
for traffic from the corresponding Boltzmann-like equation in the same spirit 
in which the derivations of the equation of continuity and Navier-Stokes equa- 
tion for viscous fluids from the Boltzmann-equation have been carried out. 
However, because of the postulate of "vehicular chaos", the equation (28) is 
expected to be valid only at very low densities where the correlations between 
the vehicles is negligibly small whereas traffic is better approximated as a 
continuum fluid at higher densities! 

Let us define the moments 

mk,e{x;t) = j dv j dvdesv''v^d^^g{x,v,Vdes;t) (32) 

Note that c = mo,o, (v) = mi^Q. Integrating the Boltzmann-like equation (28) 
over the actual velocities we get 

d d _ 

g^foi^^ '^des; t) + '^des] t)fo{x, Vdes] t)] = (33) 

where ^(a;, v^es; ^) is defined as 

-I ,^ fdvvg{x,v,Vdes;t) 

v{x,Vdes;t) = — (34) 

j0{X,Vdes;t) 

The equation (33) is an equation of continuity for each desired speed Vdes 
separately; it is consequence of the assumption that dvdes/dt — 0, i.e., no 
driver changes the desired speed. 

Using the Boltzmann-hke equation (28) and the definition (32) we can get 
separate partial differential equations for the moments of v, moments of Vdes 
and the mixed moments of v and Vdes- Unfortunately, these lead to a hierarchy 
of moment equations where each evolution equation for moment of a given 
order involves also moments of the next higher order. In order to close this 
system of equations, one needs to make appropriate justifiable assumptions. 
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6 Car-following theories of vehicular traffic 

In the car-following theories [10,11,81] one writes, for each individual vehicle, 
an equation of motion which is the analogue of the Newton's equation for each 
individual particle in a system of interacting classical particles. In Newtonian 
mechanics, the acceleration may be regarded as the response of the particle to 
the stimulus it receives in the form of force which includes both the external 
force as well as those arising from its interaction with all the other particles 
in the system. Therefore, the basic philosophy of the car-following theories 
[10,11,81] can be summarized by the equation 

[Response]n oc [Stimulus]n (35) 

for the n-th vehicle (n = 1, 2, ...). Each driver can respond to the surrounding 
traffic conditions only by accelerating or decelerating the vehicle. Different 
forms of the equations of motion of the vehicles in the different versions of the 
car-following models arise from the differences in their postulates regarding the 
nature of the stimulus (i.e., "behavioural force" or a "generalized force" [82]). 
The stimulus may be composed of the speed of the vehicle, the difference in 
the speeds of the vehicle under consideration and its lead vehicle, the distance- 
headway, etc., and, therefore, in general, 

Xn = fsti{vn, Aa;„, Avn) (36) 

where the function fsu represents the stimulus received by the n-th vehicle. 
Different versions of the car-following models model the function fsu differ- 
ently. In the next two subsections we discuss two different conceptual frame- 
works for modeling fsu- 

6.1 Follow-the-leader model 

In the earliest car-following models [83,84] the difference in the velocities of 
the n-th and (n + l)-th vehicles was assumed to be the stimulus for the n-th 
vehicleQ. In other words, it was assumed that every driver tends to move with 
the same speed as that of the corresponding leading vehicle so that 

Xnit) = - [Xn+lit) - Xn{t)] (37) 

r 

where r is a parameter that sets the time scale of the model. Note that 1/r in 
the equation (37) can be interpreted as a measure of the sensitivity coefficient 

^ In the following we label the vehicles in driving direction such that the (n + l)-th 
vehicle is in front of the n-th vehicle. 
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S of the driver; it indicates how strongly the driver responds responds to unit 
stimuhis. According to such models (and their generalizations proposed in the 
fifties and sixties) the driving strategy is to follow the leader and, therefore, 
such car-following models are collectively referred to as the follow-the-leader 
model. 

Pipes [84] derived the equation (37) by differentiating, with respect to time, 
both sides of the equation 

AXn{t) = Xn+l{t) - Xn{t) = {/^x)safe + rXn{t) (38) 

which encapsulates his basic assumption that (a) the higher is the speed of the 
vehicle the larger should be the distance-headway, and, (b) in order to avoid 
collision with the leading vehicle, each driver must maintain a safe distance 
{Ax)safe from the leading vehicle. 

It has been argued [85] that, for a more realistic description, the strength of 
the response of a driver at time t should depend on the stimulus received from 
the other vehicles at time t — T where T is a response time lag. Therefore, 
generahzing the equation (37) one would get [85] 

Xn{t + T)=S[xn+i{t)-Xn{t)] (39) 

where the sensitivity coefficient <S is a constant independent of n. 

According to the equations (37) and (39), a vehicle would accelerate or de- 
celerate to acquire the same speed as that of its leading vehicle. This implies 
that, as if, slower following vehicle are dragged by their faster leading vehi- 
cle. In these linear dynamical models the acceleration response of a driver is 
completely independent of the distance-headway. Therefore, this oversimpli- 
fied equation, fails to account for the clustering of the vehicles observed in 
real traflic. Moreover, since there is no density-dependence in this dynamical 
equation, the fundamental relation cannot be derived from this dynamics. In 
order to make the model more realistic, we now assume [86] that the closer is 
the n-th vehicle to the {n + l)-th the higher is the sensitivity of the driver of 
the n-th car. In this case, the dynamical equation (39) is further generalized 
to 

where k is a constant. An even further generalization of the the model can be 
achieved [87,88] by expressing the sensitivity factor for the n-th driver as 

4v^(t + r)r 
[x„+,it) - x„it)Y ^ ^ 
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where I and m arc phenomenological parameters to be fixed by comparison 
with empirical data. These generahzed follow-the-leader models lead to cou- 
pled non-/mear differential equations for a;„. Thus, in this "microscopic" the- 
oretical approach, the problem of traffic flow reduces to problems of nonlinear 
dynamics. 

So far as the stability analysis is concerned, there are two types of analysis 
that are usually carried out. The local stability analysis gives information on 
the nature of the response offered by the following vehicle to a fluctuation in 
the motion of its leading vehicle. On the other hand, the manner in which a 
fluctuation in the motion of any vehicle is propagated over a long distance 
through a sequence of vehicles can be obtained from an asymptotic stability 
analysis. 

From experience with real traffic we know that drivers often observe not only 
the leading vehicle but also a few other vehicles ahead of the leading vehicle. 
For example, the effect of the leading vehicle of the leading vehicle can be 
incorporated in the same spirit as the effect of "next-nearest-neighbour" in 
various lattice models in statistical mechanics. A linear dynamical equation, 
which takes into account this "next-nearest-neighbour" within the framework 
of the follow-the-leader model, can be written as [81] 

Xn{t + T) = - Xr,{t)\ + - Xn{t)\ (42) 

where S^^^ and S^^^ are two phenomenological response coefficients. 

The weakest point of these theories is that these involve several phenomeno- 
logical parameters which are determined through "calibration", i.e., by fitting 
some predictions of the theory with the corresponding empirical data. Besides, 
extension of these models to multi-lane traffic is difficult since every driver is 
satisfied if he/she can attain the desired speed! 



6.2 Optimal velocity models 



We can express the driving strategy of the driver in the car-following models 
in terms of mathematical symbols by writing 

Xnit)^-\v^''''^'{t)-vM (43) 

where l^'^''*'"^'^(t) is the desired speed of the n-th driver at time t. In all follow- 
the-leader models mentioned above the driver maintains a safe distance from 
the leading vehicle by choosing the speed of the leading vehicle as his/her own 
desired speed, i.e., v^^^'^^^'^{t) = Xn+i- 
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An alternative possibility has been explored in recent works based on the car- 
following approach [89-104]. This formulation is based on the assumption that 
ydestred (^^pQ^ids ou the dzstancc- headway of the n-th vehicle, i.e., v^e««»"e'^(^^) — 
V°P\Axn{t)) so that 

Xn{t) = - \V"P'{/^Xr,{t)) - vM (44) 

where the so-called optimal velocity function V"^^{Axn) depends on the corre- 
sponding instantaneous distance- headway Axn{t) = Xn^i{t) — Xn{t). In other 
words, according to this alternative driving strategy, the n-th vehicle tends 
to maintain a safe speed that depends on the relative position, rather than 
relative velocity, of the n-th vehicle. In general, V°p^{Ax) — > as Ax — > and 
must be bounded for Ax oo. For explicit calculations, one has to postulate 
a specific functional form of V°^'^{Ax). Car- following models along this line of 
approach have been introduced by Bando et al. [89,90]. For obvious reasons, 
these models are usually referred to as optimal velocity (OV) models. 

Since the equations of motion in the foUow-the-leader models involve only the 
velocities, and not positions, of the vehicles these can be formulated as essen- 
tially first order differential equations (for velocities) with respect to time. In 
contrast, since the equations of motion in the OV model involve the positions 
of the vehicles explicitly, the theoretical problems of this model are formu- 
lated mathematically in terms of second order differential equations (for the 
positions of the vehicles) with respect to time [89,90]. 

The simplest choice for V^^Ax) is [91,92] 

V^'iAx) ^VmaM^x-d), (45) 

where d is a. constant and is the Heavyside step function. According this 
form of y°^*(Aa;), a vehicle should stop if the corresponding distance-headway 
is less than d; otherwise, it can accelerate so as to reach the maximum allowed 
velocity Vmax- A somewhat more realistic choice [92,96] is 

for Ax < Ax A, 
y°P*(Ax) = { fAx for Axa < Ax < Axb, (46) 
for Axb < Ax. 

The main advantage of the forms (45) and (46) of the OV function is that exact 
analytical calculations, e.g. in the jammed region, are possible [92]. Although 
(45) and (46) may not appear very realistic, they capture several key features 
of more realistic forms of OV function [89], e.g., 

yP^Ax) = tanh[Ax - Ax^] + tanh[Axc] . (47) 
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for which analytical calculations are very difficult. For the convenience of nu- 
merical investigation, the dynamical equation (44) for the vehicles in the OV 
model has been discretized and, then, rewritten as a difference equation [104]. 
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Fig. 10. Fundamental diagram of the OV model. The solid line shows the OV func- 
tion and the dots simulation data. One can distinguish five different density regimes 
with respect to the stable stationary state (from [92]). 

The main question addressed by the OV model is the following: what is the 
condition for the stability of the homogeneous solution 



where b = {Ax)av = L/N is the constant average spacing between the vehicles 
and c is the constant velocity. It is not difficult to argue that, in general, in 



One can distinguish five different density regimes with respect to the stability 
of microscopic states (see Fig. 10). At low and at high densities the homoge- 
neous states are stable. For intermediate densities one can distinguish three 
regimes where jammed states exist. In region III the jammed state is stable 
whereas in regions II and IV both homogeneous and jammed states form stable 
structures. Beyond the formation of jams also hysteresis effects have been ob- 
served. Thus the OV model is able to reproduce many aspects of experimental 
findings. 

Modified Korteweg-de Vries (KdV) equation has been derived from the equa- 
tion (44) in a special regime of the parameters [93] and the relations between 
its kink solutions and traffic congestion have been elucidated [99]. A gener- 
alization to two-lane traffic can be found in [105]. In order to account for 



x„ = bn + at 



(48) 




dAx Aa;=6 
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traffic consisting of two different types of vehicles, say, cars and trucks, Mason 
et al. [97] generalized the formulation of Bando et al. [89] by replacing the 
constant r by t„ so that 



(49) 



where now depends on whether the n-th vehicle is a car or a truck. Since 
a truck is expected to take longer to respond than a car we should assign 
larger r to trucks and smaller r to cars. Some other mathematically motivated 
generalizations of the OV model have also been considered [95,98,100]. 

As mentioned earlier, drivers often receive stimulus not only from the leading 
vehicle but also a few other vehicles ahead of the leading vehicle. One possible 
way to generalize the OV models for taking into account such multi-vehicle 
interactions [102] is to write the dynamical equations as 



where Sj are sensitivity coefficients; one of the commonly used explicit forms 
of the optimal velocity function, for example (47) can be chosen for that of 
the function V°p^ in (50). 

Before concluding our discussions on the OV models let us mention a very 
simple model [94] where one assumes that 



with a > 0. Obviously, Axc may be interpreted as the safety distance. More- 
over, a restriction Vmin ^ Vn < Vmax is imposed on the allowed velocities of the 
vehicles by introducing the allowed minimum and maximum velocities Vmin 
and Vmax, respectively. Note that in this oversimplified model ,t„ depends on 
the corresponding distance-headway and, therefore, has some apparent simi- 
larities with the OV models. But, unlike the more general OV models, in this 
model Xn does not depend on its instantaneous velocity. The main reason for 
considering such a oversimplified scenario is that velocity of the propagation 
of jams can be calculated analytically. 

Some ideas of the OV model has been utilized by Mahnke and Pieret [106] in 
their master equation approach to the study of jam dynamics. They assumed 

that, at a time, only one vehicle can go into or come out of a jam; this, 
naturally, does not take into the merging or splitting of jams. Under this 
assumption, the master equation (see e.g. App. F, equation (F.l)) governing 
the probability distribution P(n; t) of the jam sizes n is given by 




(50) 




for Axn > 

for Axn < Axc- 



(51) 
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^^^^ = Vr+(n - l)P{n -l;t) + W_{n + l)P(n + 1; t) 

- [W+{n) + Pr_(n)]P(n;t) (52) 

where W+ and are the "growth" and "decay" transition rates, respec- 
tively. It has been argued [106,107] that W-{n) = 1/r = constant since a 
vehicle would require a constant average time r to come out of a jam. How- 
ever, W+{n) would depend on V^*(Aa;) since the time taken by a free-flowing 
vehicle immediately behind a jam to get into the jam would depend on Ax as 
well as on V°^*{Ax) although the actual expression of W+{n) may be compli- 
cated in a reasonable ansatz [107]. Moreover, several assumptions of the model 
will have to be relaxed before the results of this approach can be compared 
with those from real traffic. 

Before concluding this section we would like to emphasize that while formu- 
lating the dynamical equations for updating the velocities and positions of the 
vehicles in any "microscopic" theory the following points should be considered: 

(i) in the absence of any disturbance from the road conditions and interac- 
tions with other vehicles, a driver tends to drive with a desired velocity Vdes', 
if the actual current velocity of the vehicle v{t) is smaller (larger) than v^es, 
the vehicle accelerates (decelerates) so as to approach v^es- 

(ii) In a freely-flowing traffic, even when a driver succeeds in attaining the 
desired velocity Vdes, the velocity of the vehicle fluctuates around Vdes rather 
than remaining constant in time. 

(iii) The interactions between a pair of successive vehicles in a lane cannot 
be neglected if the gap between them is shorter than Vdes', in such situations 
the following vehicle must decelerate so as to avoid collision with the leading 
vehicle. 

Clearly, the reliability of the predictions of the OV model depends on the 
appropriate choice of the optimal velocity function. 



7 Coupled-map lattice models of vehicular traffic 



Recall that, in the car-following models, space is assumed to be a contin- 
uum and time is represented by a continuous variable t. Besides, velocity 
and acceleration of the individual vehicles are also real variables. However, 
most often, for numerical manipulations of the differential equations of the 
car-following models, one needs to discretize the continuous variables with 
appropriately chosen grids. In contrast, in the coupled-map lattice approach 
[108], one starts with a discrete time variable and, the dynamical equations 
for the individual vehicles are formulated as discrete dynamical maps that 
relate the state variables at time t with those at time t + 1, although position. 
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velocity and acceleration are not restricted to discrete integer values. The unit 
of time in this scheme (i.e., one time step) may be interpreted as the reaction 
time of the individual drivers as the velocity of of a vehicle at the time step t 
depends on the traffic conditions at the preceding time step t — 1. 

Keeping in mind the general points raised at the end of the preceding section 
regarding the formulation of the dynamical equations for updating the veloci- 
ties and positions of the vehicles, the general forms of the dynamical maps in 
the coupled-map lattice models can be expressed as: 

Vn{t + 1) ^ Mapn[Vn{t),Vaes: ^Xn{t)], (53) 
Xnit+l)^Vnit)+Xnit) (54) 

where Vdes is a desired velocity. In general, the dynamical map Map[vn{t),Vdes, Aa;„ 
takes into account the velocity Vn{t) and the distance- headway Axn{t) of the 
n-th vehicle at time t for deciding the velocity Vn{t + 1) at time t + 1. The ef- 
fects of the interactions among the vehicles enter into the dynamical updating 
rules (53), (54) only through the distance-headway Axn- 



7.1 The model of Yukawa and Kikuchi 

Yukawa and Kikuchi [109-111] have studied coupled- map models based on the 
map 

v{t + 1) = F{v{t)) := -fv{t) + /3tanh ~^^^^ ) + ^ 

for the uninfluenced motion of a single vehicle, is the preferred velocity of 
the vehicle and 7, 6 and e are parameters. For 7 close to 1 the map becomes 
chaotic, but acceleration and deceleration are approximately constant far from 
v^. Their magnitude is determined by the parameter [3. e controls the differ- 
ence of the acceleration and deceleration capabilities. Although the model is 
deterministic, fluctuations in the velocity are introduced through determin- 
istic chaos. These fluctuations around are determined by the parameter 
5. 

If there is more than one vehicle on the road one needs an additional decel- 
eration mechanism to avoid collisions. This can be achieved by introducing 
a deceleration map. Assuming that deceleration is dominated by the head- 
way, two models have been studied in [109]. In model A, the deceleration map 
describes a sudden braking process. If the front-bumper to front-bumper dis- 
tance Axn to the next vehicle ahead is less than the current velocity Vn{t) of 
the following vehicle, then the velocity is reduced to Ax„ — / where / is the 
length of the vehicles. The corresponding map is B{Axn{t)) = Axn{t) — I. 
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Model B has a more complex deceleration map: 



Vn{t+l) = G{AXn{t),Vn{t)) 



[AXnit) - I - Vn{t)] 



{a - l)Vn{t) 
(for Vnit) < Axn{t) - I < avn{t)). (56) 



The parameter a determines the range within which the deceleration map 
G{Ax, v) is used. For headways less than aVn{t) the map G is used instead of 
F{v). Note that G{Ax,v = Ax - I) = Ax - I and G{Ax,v = {Ax - I) /a) = 
F[v = {Ax — I) /a), i.e. G interpolates between the free-motion map F and 
the sudden braking map of model A. The full velocity map of Model B is thus 
given by 



Mapr,{Vn{t),AXn{t)) 



'F{vn{t),v^) for avn{t) < Axn{t), 

G{AXn{t),Vn{t)) for Vn{t) < AXn{t) < aVn{t)), 
B{^^n{t)) for AXn{t) < Vn{t)). 



(57) 



Local measurements of the flow for a system of vehicles with different preferred 
velocities produce a fundamental diagram of inverse-A shape (see Fig. 4A) 
[111]. Here the non-uniqueness of the flow has a simple explanation. Due to 
the different v^, platoons form behind the slowest vehicles. Whenever such a 
platoon passes the measurement region, a flow value on the lower branch is 
recorded. Otherwise, the flow corresponds to the upper branch. 

Measurements of the power spectral density of temporal density fluctuations, 
i.e. the Fourier transform of the time-series of local densities, show a 1//"- 
behaviour with a ~ 1.8 in the free-flow regime. Due to the deterministic 
dynamics, the system evolves into a state with power-law fluctuations. In 
[110,111] it has been suggested that the origin of the l//"-fluctuations is 
the power-law distribution (cx l/{Ax)^'^) of the headways Ax, since these are 
related to density waves. The occurence of jams destroys long-time correlations 
since vehicles loose their memory of current fluctuations when they are forced 
to stop in a jam [110]. Therefore no 1/ /"-behaviour can be observed in the 
jammed regime. 

In [112] a coupled- map model based on optimal- velocity functions has been 
introduced by discretizing the time variable of the OV model (see Sec. 6.2). 
This allows to study systems with open boundaries and multilane systems. 
Furthermore a multiplicative random noise can be imposed in the velocity 
update so that the velocity map is given by 

V{t + 1) = [v{t) + a {Vopt{Ax) - V{t))] (1 + fnoiseO (58) 
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where ^ e [—1/2, 1/2] is a uniform random variable and f noise the noise level. 
7.2 The model of Nagel and Herrmann 

Nagel and Herrmann (NH) [113] have introduced a coupled-map model which 
is related to the continuum limit of the Nagel-Schreckenberg cellular automata 
model (see Sec. 8). A generahzation of the NH model has later been presented 
in [114]. Vehicles are characterized by a maximal velocity Vmax ^-nd a safety 
distance a. The velocity map for the NH model is given by 



In the velocity update step, vehicles which have a headway Ax smaller than 
the safety distance a decelerate. The headway distance after deceleration is 
determined by the parameter S. Vehicles which have a large enough headway, 
on the other hand, accelerate. The acceleration coefficient a is determined by 
(l,Aa;,(t)/7). 

Since the dynamics of the model is deterministic, the behaviour depends 

strongly on fluctuations of the initial state [114]. For equidistant starting po- 
sitions of the vehicles the fundamental diagram consists of two linear branches 
with maximum flow f{cm) = VmaxCm at density Cm = ^/{vmax + P)- For homo- 
geneous starting positions the system is free-flowing up to a critical density 
Ccrit- Beyond this density free-flowing and congested areas coexist. 

7.3 The model of Krauss, Wagner and Gawron 

Krauss ct al. [115,116] introduced a whole class of stochastic models by consid- 
ering necessery conditions for the collision-free motion of vehicles. The models 
are continuous in space and discrete in time. 

The vehicles are characterized by a maximum velocity v,nax, their acceleration 
and deceleration capabilities a{v) and b{v), respectively, and their length / 
which will be taken to be / = 1 in the following. Then the update rules for the 
velocity v and the space coordinate x of each vehicle are as follows: 

Step 1: Determine desired velocity. 




max(Ax„(i) — 5, 0) for Vn{t) > Ax„(i) — a, 

TSm\{Vn{t) a, Vmax) for Vn{t) < AXn{t) - (3. 



(59) 



Vdes = min [Vmax, V + aiv),Vsafe\ 
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Step 2: Randomization. 

V = max [0, rand{vdes - a, Vdes)] 



Step 3: Vehicle movement. 



X 



X + V 



Here rand{vi,V2) denotes a random number uniformly distributed in the in- 
terval [vi,V2) and Vgafe is a velocity which guarantees collision motion of the 
vehicles. It is given explicitly by 



^safe 



V + b(v) 



(60) 



where Vp is the velocity of the preceding vehicle located at Xp and g — Xp—x — 1 
is the headway, i.e. the distance to the preceding vehicle. 

In the simplest case the acceleration and deceleration capabilities do not de- 
pend on the velocity, i.e. a{v) = a — const and b{v) — b — const. The 
behaviour of the model can be classified in three different families (see Fig. 
11). 




Fig. 11. Schematic sketch of the three different classes of models in the KWG model. 

The three families of models sketched schematically in Fig. 11 can be charac- 
terized as follows: 



Class I: High acceleration 

Here no spontaneous jamming exists. For a — > Vmax and 6^1 the behaviour 
is similar to that of a cellular automata model without velocity memory in- 
troduced in [117] which is closely related to the Kasteleyn model [118] of 
statistical physics. It can also be interpreted as 5- vertex model [119]. Another 
model belonging to this class has been introduced by Fukui and Ishibashi [120] 
(see Sec. 11.1). 

Class II: High acceleration - low deceleration 

The outflow from a jam is identical to the maximal possible flow. The jamming 
transition is not a true phase transition, but rather a crossover. The limit 
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6 — > oo, a = 1 corresponds to a continuum version of the Nagel-Schreckenberg 
cellular automata model [20] which will be introduced in the next section. 

Class III: Low acceleration - low deceleration 

These models exhibit phase separation and met ast ability. The jamming tran- 
sition is of first order. The outflow from a jam is not maximal. For a <^ Vmax 
and h <^ Vmax the model is closely related to the Gipps model [88] discussed 
in Sec. 6.1. Other models belonging to this class are the Kerner-Konhauser 
model (see Sec. 4.4), the optimal- velocity model (Sec. 6.2) and the models 
with slow-to-start rules which will be introduced in Sec. 9.1.2. 

On a macroscopic level, classes I, II and III can be distinguished by the or- 
dering of the densities c/ and Cc, where c/ is the density of the outflow from a 
jam and Cc is the density where homogeneous flow becomes unstable [116]. For 
Cc > Of the outflow from a jam is stable and the system phase-separates into 
free- flow and jammed regions. Furthermore, metastable states can be found. 
This is the type of the behaviour found in class III. For Cc < Cf, on the other 
hand, the outflow from a jam is unstable and no metastable states or phase 
separation can be found. This is the typical behaviour of classes I and II. 
These classes can further be distinguished since in class I one does not find 
any structure formation, like spontaneous jamming, in contrast to class II. 

A related model has been studied before by Migowsky et al. [121]. In this model 
vehicles are also characterized by a maximum velocity Vmax and a bounded 
acceleration capability {—bmax < x < amax) which determines the safety dis- 
tance dg necessary to avoid accidents. The investigations in [121] focussed on 
the effect of so-called driving strategies. These strategies are characterized by 
a vector (/^,, /q, Z^), where and fs are the fraction of the vehicle's max- 

imal velocity, acceleration and safety distance actually used[^, respectively. 
This can lead to the possibility of accidents and allows to study the number 
of crashes as a function of the driving strategies. 

In [121] different strategies have been compared. Furthermore, dynamical 
changes of strategies can be introduced which allow the drivers to adapt to 
the local traffic conditions. In general this leads to a decrease in the number of 
accidents, jams and fuel consumption, but at high densities the flow is reduced 
compared to the case of fixed strategies. 



^ In other words, the maximum velocity is given by Vmax = fv Vmax etc. 
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8 Nagel-Schreckenberg cellular automata model of vehiculcir traffic 
on highways 

In general, CA are idealization of physical systems in which both space and 
time are assumed to be discrete and each of the interacting units can have only 
a finite number of discrete states. Note that for a discretization of differential 
equations, e.g. those of the hydrodynamic approach, space and time variables 
are discrete, but the state variable still is continuous. The concept of CA was 
introduced in nineteen fifties by von Neumann while formulating an abstract 
theory of self-replicating computing machines [122]. However, it received the 
attention of a wider audience in the nineteen seventies through Conway's game 
of life [123]. The family of one- dimensional CA was studied systematically, 
in the nineteen eighties, from the point of view of dynamical systems and 
popularized by Wolfram [18,124]. Since then the concept of CA has been 
applied to model a wide variety of systems [125-127]. To our knowledge, the 
first CA model for vehicular traffic was introduced by Cremer and Ludwig 
[128]. 

In the CA models of traffic the position, speed, acceleration as well as time 
arc treated as discrete variables. In this approach, a lane is represented by a 
one-dimensional lattice. Each of the lattice sites represents a "cell" which can 
be either empty or occupied by at most one "vehicle" at a given instant of time 
(see Fig. 12). At each discrete time step t ^ t + 1, the state of the system is 
updated following a well defined prescription (a summary of various possible 
different schemes of updating is given in Appendix A). The computational 
efficiency of the discrete CA models is the main advantage of this approach 
over the car-following and coupled-map lattice approaches. 
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Fig. 12. A typical configuration in the NaSch model. The number in the upper 

right corner is the speed of the vehicle. 

In the NaSch model, the speed v of each vehicle can take one of the Vmax + 1 
allowed m^eg'er values v = 0,1, ...,Vmax- Suppose, a;„ and v„ denote the position 
and speed, respectively, of the n-th vehicle. Then, dn = Xn+i — Xn, is the gap 
in between the n-th vehicle and the vehicle in front of it at time t. At each 
time step t — t + 1, the arrangement of the N vehicles on a finite lattice of 
length L is updated in parallel according to the following "rules": 

Step 1: Acceleration. If Vn < Vmax, the speed of the n-th vehicle is increased 
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by one, but Vn remains unaltered if 

Vn min(v„ + l,v 

max ) (ui). 

Step 2: Deceleration (due to other vehicles). If dn < Vn, the speed of the n-th 
vehicle is reduced to d„ — 1, i.e., 

Vn min(t'„, dn - 1) (U2). 

Step 3: Randomization. If f n > 0, the speed of the n-th vehicle is decreased 
randomly by unity with probability p but Vn does not change if Vn — 0, i.e., 

Vn max(i>„ — 1,0) with probability p (U3). 

Step 4-' Vehicle movement. Each vehicle is moved forward according to its new 
velocity determined in Steps 1-3, i.e. 

Xn^Xn + Vn (U4). 



The NaSch model is a minimal model in the sense that all the four steps are 
necessary to reproduce the basic features of real traffic; however, additional 
rules need to be formulated to capture more complex situations. The step 
1 reflects the general tendency of the drivers to drive as fast as possible, if 
allowed to do so, without crossing the maximum speed limit. The step 2 is 
intended to avoid collision between the vehicles. The randomization in step 3 
takes into account the different behavioural patterns of the individual drivers, 
especially, nondeterministic acceleration as well as overreaction while slowing 
down; this is crucially important for the spontaneous formation of traffic jams. 
Even changing the precise order of the steps of the update rules stated above 
would change the properties of the model. E.g. after changing the order of 
steps 2 and 3 there will be no overreactions at braking and thus no spontaneous 
formation of jams. The NaSch model may be regarded as stochastic CA [18]. 
In the special case v^ax = 1 the deterministic limit of the NaSch model is 
equivalent to the CA rule 184 in Wolfram's notation [18] and some abstract 
extensions of this CA-184 rules [129] have been studied in the more general 
context of complex dynamics and particle flow. 

Why should the updating be done in parallel, rather than in random sequential 
manner, in traffic models like the NaSch model? In contrast to a random se- 
quential update, parallel update can lead to a chain of overreactions. Suppose, 
a vehicle slows down due the randomization step. If the density of vehicles is 
large enough this might force the following vehicle also to brake in the decel- 
eration step. In addition, if p is larger than zero, it might brake even further 
in Step 3. Eventually this can lead to the stopping of a vehicle, thus creating 
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a jam. This mechanism of spontaneous jam formation is rather reahstic and 
cannot be modeled by the random sequential update. 

The update scheme of the NaSch model is illustrated with a simple example 
in Fig. 13. 
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Fig. 13. Step-by-step example for the application of the update rules. Wc have 
assumed Vmax = 2 and p = 1/3. Therefore on average one third of the cars qualifying 
will slow down in the randomization step. 

Space-time diagrams showing the time evolutions of the NaSch model demon- 
strate that no jam is present at sufficiently low densities, but spontaneous 
fluctuations give rise to traffic jams at higher densities (Fig. 14(a)). From the 
Fig. 14(b) it should be obvious that the intrinsic stochasticity of the dynamics 
[20], arising from non-zero p, is essential for triggering the jams [20,113]. For 
a realistic description of highway traffic [20], the typical length of each cell 
should be about 7.5m which is the space occupied by a vehicle in a dense jam. 
When Vmax — 5 each time step should correspond to approximately 1 sec of 
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Fig. 14. Typical space-time diagrams of the NaSch model with Vmax = 5 and (a) 
p = 0.25, c = 0.20, (b) p = 0, c = 0.5. Each horizontal row of dots represents the 
instantaneous positions of the vehicles moving towards right while the successive 
rows of dots represent the positions of the same vehicles at the successive time steps. 

real time which is of the order of the shortest relevant timescale in real traffic, 
namely the reaction time of the drivers. 

Almost all the models of traffic considered in this review, including the NaSch 
model, have been formulated in such away that no accident between successive 
vehicles is possible. However, accident of the vehicles is possible if the condition 
for safe driving is relaxed. For example, Boccara et al. [130] replaced the 
update rule of the NaSch model by the rule 

if Vn+i{t) > then Xn{t + 1) = Xn{t) + Vn{t + 1) + Av (61) 

where Av is a Bernoulli random variable which takes the value 1 with proba- 
bility Pcareiess and zero with the probability 1 —Pcareiess- The probability Pac of 
accident per vehicle per time step is a non-monotonic function of the vehicle 
density c [130,131]. 

8.1 Relation with other models 

8.1.1 Relation with totally asymmetric simple exclusion process 

Now we point out the similarities and differences between the Vmax = 1 
limit of the NaSch model and the totally asymmetric simple exclusion pro- 
cess (TASEP) which is the simplest prototype model of interacting systems 
driven far from equilibrium [16,14,15]. In the TASEP (Fig. 15) a randomly cho- 
sen particle can move forward, by one lattice spacing, with probability q if the 
lattice site immediately in front of it is empty. It corresponds to a Kawasaki dy- 
namics [132] for exchange of a charged particle and hole on nearest-neighbour 
lattice sites at infinite temperature and in the presence of an infinite electric 
field [133] (see Appendix B for some technical aspects of TASEP). Several 
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different generalized variants of tlie TASEP liave been considered. For exam- 
ple, in the fc-hop model [134] a particle can exchange its position with the 
nearest-hole on its right with probability q, provided the separation of the two 
sites under consideration is not more than k lattice spacings. The k-hop model 
reduces to the TASEP in the special case k — 1. 



C5 o 



Fig. 15. ASEP with periodic boundary conditions. 

Note that in the NaSch model with Vmax = 1 every vehicle moves forward with 
probability q = 1 — p in the time step t + 1 if the site immediately in front 
of it were empty at the time step t; this, is similar to TASEP. But, updating 
is done in parallel in the NaSch model whereas that in the TASEP is done 
in a random sequential manner. Nevertheless, the special case of Vmax — 1 
for the NaSch model achieves special importance from the fact that so far it 
has been possible to derive exact analytical results for the NaSch model only 
in the special limits (a) p = and arbitrary Vmax (which we have already 
considered), and (b) Vmax = 1 and arbitrary p. 



8.1.2 Relation with surface growth models and the phenomenological fluid- 
dynamical theories of traffic 

The NaSch model with Vmax — 1 can be mapped onto stochastic growth mod- 
els of one-dimensional surfaces in a two-dimensional medium, the single-step 
model [135]. Corresponding to each configuration {crj} of the NaSch model 
in the site-oriented description, one can obtain a unique surface profile {Hj} 
through the relation Hj — | I]j<fc(l — 2(7^) [135,136]. Pictorially one can inter- 
pret this mapping as shown in Fig. 16. Particle (hole) movement to the right 
(left) correspond to local forward growth of the surface via particle deposition. 
In this scenario a particle evaporation would correspond to a particle (hole) 
movement to the left (right) which is not allowed in the NaSch model. It is 
worth pointing out that any quenched disorder in the rate of hopping between 
two adjacent sites would correspond to columnar quenched disorder in the 
growth rate for the surface [137]. 

The surface growth model described above is known to be a discrete counter- 
part of continuum models of growing surfaces whose dynamics are governed 
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Fig. 16. The schematic pictorial representation of the mapping of the NaSch model 
with Vmax = 1 onto a stochastic model of surface growth. 

by the so-called Kardar-Parisi-Zhang equation [135,136]. Since the Kardar- 
Parisi-Zhang equation can be mapped onto the Burgers equation [53] using 
the Cole-Hopf transformation [135], it is not surprising that several features 
of vehicular traffic are described by the NaSch model at the microscopic level 
and by the noisy Burgers equation for the coarse-grained continuum of the 
fluid-dynamical theory [138] . 



8.2 Limiting cases of the NaSch model 



In spite of the fact that the deterministic limits p = and p = 1 of the NaSch 
model do not capture some of the most essential features of vehicular traffic 
it may be instructive to examine these limits to gain insight into the features 
of this simpler scenario. Another limiting case which exhibits a surprisingly 
complex behaviour is the case Vmax — oo. 



8.2.1 NaSch model in the deterministic limit p — 

The NaSch model, a stochastic CA, becomes a deterministic CA in the limit 
p = 0. In this special case, the deterministic dynamical update rules of the 
model can be written as 



Vnit + 1) = min[f„a^, Vnit) + 1, - 1], 

Xn{t^l)=Xn{t)+Vn{t + l) (62) 

which can lead to two types of steady states depending on the density c [113]. 
At low densities, the system can self-organize so that dn > v^ax foi' all n and, 
therefore, every vehicle can move with v^ax-i i-e., Vn{t) = Vmax, giving rise to 
the corresponding flux cumax- This steady-state is, however, possible only if 
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enough empty cells are available in front of every vehicle, i.e., for c < cj^* = 
{vmax + and the corresponding maximum flux is J^ax = '^max/ (vmax + !)• 
On the other hand, for c > c^*, — 1 < min[t>„(t) + 1, Vmax] and, therefore, 
the relevant steady-states are characterized by Vn(t) = dn(t) — l, i.e., the flow is 
limited by the density of holes. Since the average distance-headway is 1/c — 1, 
the fundamental diagram in the deterministic limit p = of the NaSch model 
(for any arbitrary Vmax) is given by the exact expression 

J = mm{cvmax, 1-c). (63) 

Note that the result = 1/c — 1 is identical with Greenshields ansatz f „ = 
1/c - l/cjam if we identify Cjam = 1- 

8.2.2 NaSch model in the deterministic limit p = 1 

Aren't the properties of the NaSch model with maximum allowed speed Vmax, 
in the deterministic limit p = 1, exactly identical to those of the same model 
with maximum allowed speed Vmax — 1? Although this expectation may seem 
to be consistent with the observation that J = for all c in the special case 
Vmax = 1 = the answer to the question posed above is: NO. To understand 
the subtle features of the deterministic limit p = 1 one has to consider Vmax > 
1. You can easily convince yourself that if, for example, Vmax = 2, then, for 
c > 1/3, all stationary states correspond to J = because at least one vehicle 
will have only one empty cell in front (i.e. dn = 2) and it will never succeed in 
moving forward. For Vmax = 2 and p = 1, although there are stationary states 
corresponding to J 7^ for all c < 1/3, such states are metastable in the sense 
that any local external perturbation leads to complete breakdown of the flow. 
If the initial state is random, such metastable states cannot lead to non-zero 
J because they have a vanishing weight in the thermodynamic limit. Hence, 
if p = 1, all random initial states lead to J = in the stationary state of the 
NaSch model irrespective of Vmax and c! 

8.2.3 NaSch model in the limit Vmax = 00 

The limit Vmax = 00 has been introduced in [139]. One has to be aware that 
there are several possible ways of performing this limit since only finite systems 
of length L can be treated in computer simulations. In [139] the case Vmax = L 
has been investigated!^ but other limiting procedures are also possible, e.g. 
Vmax oc with a > or even Vmax = oo independent of the system size. 
In principle, these different limiting procedures could lead to different results, 
but up to now no systematic study has been performed. We therefore restrict 
ourselves to the case Vmax = L studied in [139]. 

See also [113], where the case p = was studied. 
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Surprisingly one finds that the fundamental diagram has a form quite different 
from that of the case of a finite v^ax- The fiow does not vanish in the limit c — > 
since already one single car produces a finite value of the flow, J(c — > 0) = 1. 
Due to the hindrance effect of other cars, J(c) is a monotonically decreasing 
function of the density c (see Fig. 17). Another characteristic feature of the 
fundamental diagram is the cxistancc of a plateau at flow Jp where the value 
Jp depends on the randomization p, but not on the system size L. The length 
of the plateau, on the other hand, increases with L. 
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Fig. 17. Fundamental diagram of the NaSch model with Vmax = L for a system size 
of L = 1000 (from [139]). 

What is the microscopic structure of the stationary state leading to such a 
fundamental diagram ? At low densities, where flow J is larger than the plateau 
value Jp, the cars tend to be uniformly distributed just as in the deterministic 
case p = (see Sec. 8.2.1). For densities in the plateau regime, however, one 
jam exists in the system, whereas for higher densities there is more than one 
jam. In the thermodynamic limit, one expects a phase transition at c = 
between a jamless phase with J = 1 and a phase with one jam and flow 
Jp [139]. Increasing the density further, more jams develop and the plateau 
ceases. Note that this behaviour is completely different from the prediction of 
mean-field theory in that limit [140] (see Sec. 8.3.1) showing the importance 
of correlations. 

8.3 Analytical theories of the NaSch model with periodic boundary conditions 

CA are, by design, ideal for large scale computer simulations. However, proper 
interpretations of the numerical data obtained from computer simulations are 
not always quite straightforward because of the finite-size effects and "nu- 
merical noise" . One cannot deny the importance of exact analytical results in 
providing a testing ground for the computer codes. On the other hand, the 
parallel updating makes exact analytical solution of CA models very difficult. 
Nevertheless, even in those situations where exact solutions are not possible. 
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a combination of approximate analytical treatments and computer simulation 
often turns out out to be very powerful method of analysis of a problem. This 
approach has been quite successful in recent years in the studies of the NaSch 
model and its generahzations. 

Before we proceed with the analytical theories in the non-deterministic NaSch 
model, we would like to point out that the fundamental diagram J(c) is known 
exactly for arbitrary v^ax and p in the two limits c ^ and c ^ 1. In the 
former case, J ~ cvp where vp = Vmax — p is the free-flow velocity. On the 
other hand, in the latter case, J ~ (1 — — c) as flow is determined by 
holes moving backwards at a speed 1 — p. 



8.3.1 Site- oriented naive mean-field theory for the NaSch model 

In the "site-oriented" theories one describes the state of the finite system 
of length L by completely specifying the state of each site, i.e., by the set 
((Ti, (72, (Jl) where aj {j — 1,2, L) can, in principle, take Vmax + 2 values 
one of which represents an empty site while the remaining Vmax + 1 correspond 
to the Vmax + 1 possible values of the speed of the vehicle occupying the site j. 
In some of the analytical calculations of steady-state properties of the NaSch 
model one follows, for convenience, the sequence 2 — 3 — 4 — 1, instead of 
1 — 2 — 3—4 of the stages of updating [140] as this merely shifts the starting step 
and, therefore, does not influence the steady-state properties of the model. The 
advantage of this new sequence is that, in a site-oriented theory, the variable 
aj can now take Vmax + 1 values as none of the vehicles can have a speed v — 
at the end of the acceleration stage of the updating. 

Let us introduce the lattice gas variables n{i;t) through the following defini- 
tion: n{i; t) = if the site labeled by i is empty and n{i; t) = 1 if it is occupied 
by a vehicle (irrespective of the speed). Obviously, the space-average of n{i;t) 
is the density of the vehicles, i.e., J2in{i;t)/ L — c. Suppose, Cv{i;t) is the 
probability that there is a vehicle with speed v {v — 0,1,2, Vmax) at the site 
i at the time step t. Obviously, c{i;t) = J2l='o' Cv{i',t) is the probability that 
the site i is occupied by a vehicle at the time step t and d{r, t) = 1 — c{i; t) is 
the corresponding probability that the site i is empty at the time step t. 

In the naive site-oriented mean-field (SOMF) approximation for the NaSch 
model one writes down the equations relating Cy{i;t-\-l) {v — 1, v^ax) with 
the corresponding probabilities at time t and, then, solves the equations in the 
steady-state (see the Appendix C for a detailed derivation of these equations 
for arbitrary Vmax)- In the simplest case of v^ax = 1 and periodic boundary 
conditions one gets [140] 
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co{i; t + l) = c{i; t)c{i + 1; t) + pc(i; t)d{i + 1; t), 
Ciii] if: + 1) = qc{i — 1; t)d{i] t). 



(64) 
(65) 



The equation (64) expresses the simple fact that at the time step t + 1 the 
speed of the vehicle at the i-th site can be zero either because the site i + 1 
was occupied at time t or because of random deceleration (if the site i + 1 was 
empty at time t). Similarly, the equation (65) implies that the speed of the 
vehicle at the site i can be 1 at time t + 1 if at the time step t the site i was 
empty while the site i — 1 was occupied by a vehicle which did not decelerate 
during the random deceleration stage of updating. 

In the steady state, Cy{i,t) are independent of t. Besides, if periodic bound- 
ary conditions are imposed, the i-dependence of c„(i) also drops out in the 
translational-invariant steady-state. Therefore, in the steady-state 

J = ci = gc(l — c) (66) 

It turns out [140] that the naive SOMF underestimates the flux for all Vmax- 
Curiously, if instead of parallel updating one uses the random sequential up- 
dating, the NaSch model with Vmax = 1 reduces to the TASEP for which the 
equation (66) is known to be the exact expression for the corresponding flux 
(see, e.g., [20])! 

8.3.2 Paradisical mean- field theory of the NaSch model 

What are the reasons for these differences arising from parallel updating and 
random sequential updating? There are "Garden of Eden" (GoE) [141] states 
(dynamically forbidden states) [142] of the NaSch model which cannot be 
reached by the parallel updating whereas no state is dynamically forbidden 
if the updating is done in a random sequential manner. For example, the 
configuration shown in Fig. 18 is a GoE stateP^ because it could occur at 
time t only if the two vehicles occupied the same cell simultaneously at time 
t - 1. 





1 


2 















Fig. 18. A GoE state for the NaSch model with Vmax ^ 2. 

The naive SOMF theory, discussed in the preceding subsection, does not ex- 
clude the GoE states. On the other hand, results of the paradisical mean-field 

The configuration shown in Fig. 12 is also a GoE state! 
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(pMF) theory are derived by repeating the calculations of the naive SOMF 
theory excluding all the GoE states from consideration. The exact expression, 
given in the next subsection, for the flux in the steady-state of the NaSch 
model with v^ax = 1 is obtained in the pMF theory (see Appendix D for 
detailed calculations), thereby indicating that the only source of correlation 
in this case is the parallel updating [142]. But, for Vmax > 1, there are other 
sources of correlation because of which exclusion of the GoE states merely 
improves the naive SOMF estimate of the flux (Fig. 19) but does not yield 
exact results [142,143]. 




Fig. 19. Fundamental diagram for Vmax = 2 and p = 0.5 (left) and p = 0.1 (right). 
Comparison of paradisical MFT (full line) with results from computer simulations 
(•) and the naive MFT (dotted line). 

8.3.3 Site- oriented cluster-theoretic approach to the NaSch Model 

The site-oriented cluster theoretic approach leads to a systematic improvement 
of the naive SOMF theory of the NaSch model. We define a n-cluster to be 
a collection of n successive sites and denote the probability of finding an n- 
cluster in the state (cti, a"2, a"„) in the steady-state of the system by the 
symbol P„((Ti, (T2, o"„). In the general n-cluster approximation [140], one 
divides the lattice into "clusters" of length n such that two neighbouring 
clusters have n — 1 sites in common (see Fig. 20); an n-cluster is treated 
exactly and the cluster is coupled to the rest of the system in a self-consistent 
way, as we shall show in this subsection. Even without any calculation, one 
would expect that, for a given Vmaxi the n-cluster approximation should yield 
more accurate results with increasing n and should give exact results in the 
limit n ^ 00. Fortunately, in the special case Vmax = 1 exact results are 
obtained already for n as small as 2, i.e., the results of 2-cluster calculations 
are exact for Vmax — 1 [140] . 

Let us first explain the key concepts involved in the cluster theory. It is 
straightforward to verify, for example, in the special case of Vmax = 1; that the 
state of the 2-cluster Uj, CTj+i at time t-\-l depends on the state of the 4-cluster 
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Fig. 20. Decomposition of a lattice into (a) 1-clusters (b) 2-clusters and (c) 3-clusters 
in the SOMF theory. 



(tj-i, Tj, Tj+i, Ti+2) at time t. In general, in the n-cluster approximation for an 
arbitrary v^ax one has to take into account the vehicles that can enter an 
n-clustcr from one of the Vmax cells to its left and can leave it to occupy one 
of the Vmax cells to its right. Therefore, in general, the state cxj, cxj+i, ...cxj+n-i 
of an n-cluster at time t -\- 1 depends on the state of a n + 2vmax cluster 
^j-vma.^ ^3+1^ -.^j+n-i+vma. ^t time t. Thcrefore, in the special 

case of Vmax — 1) the master equations 



P2{cri, (Jj+i; i + 1) = ^ W{(Ti, (7i+i|Tj_i, Tj, Tj+i, ri+2)-P4(Ti_l, Tj, Tj+i, ri+2; t) 

(67) 



governing the time evolution of the 2-cluster probabilities P2(ci, '^i+i) involve 
the 4-cluster probabilities for all those configurations (rj_i, Tj, Tj+i, rj_|_2; t) which 
can lead to the 2-cluster configuration (cij, (jj+i; t -|- 1) under consideration as 
well as the corresponding transition probabilities W{ai, (Ti+i|Tj_i, r^, r^+i, Ti+2). 
Similarly, the master equation governing the time evolution of the 4-cluster 
probabilities i-*4(Tj_i, Tj, Tj+i, Tj+2) involve 6-cluster probabilities, and so on. In 
order to obtain a closed set of equations one has to truncate this hierarchy 
in an appropriate manner; in the n-cluster approximation one expresses the 
(n-|-2t'maa;)-cluster probabilities in terms of products of n-cluster probabilities. 

The n-cluster approximation represented geometrically in Fig. 20 for n = 1 
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can be expressed mathematically as 

i=Tj-2 

Thus, 1-cluster approximation is equivalent to the naive SOMF approxima- 
tion. The 2-cluster approximation represented geometrically in Fig. 20 can be 
expressed mathematically as [140] 



Tj-l, Tj, Tj+i, Tj+2, Tj+s) OC -P2(Tj-2, Tj_i)P2(rj-l, Tj)P2{Tj, Tj+i) 

•^2(Tj+1, T;-+2)P2(Tj+2, Tj+3) (69) 

or, more precisely, 

P{^j-2, Tj-l, Tj, Tj+i, Tj+2, Tj+s) = P2{Tj_2\Tj_i)P2{Tj-i \Tj)P2{Tj, Tj+i) 

■P2{Tj+^\Tj+2)P2{Tj+2\Tj+3) (70) 

where 

P2{Tj-l\Tj) = — — (71) 

Z^Tj-i ^2[Tj-l,Tj) 

are 2-cluster conditional probabilities. Similarly, the 3-cluster approximation 
consists of the approximate factorization 

Tj, Tj+i, Tj+2, Tj+3, Tj+4) = P3{Tj_2\ rj-l, Tj ) P3{Tj _i\ Tj , Tj+i ) 
■Psi^j, ^j+l, ^j+2)P3{ rj+l, rj+2 ki+3)^3( ^j+2: ~^-+3 kj+4)- (72) 

Analoguous factorizations hold for an arbitrary number of sites on the left- 
hand-side of (68), (70) and (72). 

Let us now illustrate the scheme of the cluster calculations for the NaSch 
model by carrying out the calculation for the simplest case, namely, the 2- 
cluster calculations for Vmax = 1- For convenience, one follows the sequence 
2—3—4—1, instead of 1—2— 3— 4 of the stages of updating so that, for f^aa; = 1, 
a two-state variable a is adequate to describe the state of a lattice site; u = 0, 1 
correspond, respectively, to an empty site and a site occupied by a vehicle with 
speed 1. Thus, in the special case of Vmax = 1 one would need only four 2- 
cluster probabilities, namely, ^2(0, 0), ^2(1, 0), ^2(0, 1), -P2(l, !)• Interestingly, 
the constraints 

P2(l,0) + P2(l,l) = c (73) 
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and 



P2(0,0)+P2(0,l) = 1-C 



(74) 



together with the particle-hole symmetry 



P2(1,0)=P2(0,1) 



(75) 



leave only one of the four 2-cluster probabilities, say, P2(l,0), as an inde- 
pendent variable which one needs to calculate by solving the corresponding 
master equation. For general Vmax, the n-cluster approximation on the right 
hand side of the master equation leads to {vmax + 1)" nonhnear equations; the 
number of independent equations gets reduced by the so-called Kolmogorov 
consistency conditions [144]. 

Using (70) one factorizes the 4-cluster probabihties on the right-hand-side of 
(67) for P2(l, 0) in terms of 2-cluster conditional probabilities. In the first col- 
umn of the table in Fig. 21 we list all those configurations (Ti_i, Tj, Tj+i, Tj+2; t) 
which can lead to the configurations, shown in the second column, which is 
the exhaustive list of the 4-cluster configurations each having ai — 1, (Tj+i = 0; 
the corresponding transition probabilities VF(1, 0|Ti_i, Tj, Tj+i, Tj+2) are given 
in the third column. 



Fig. 21. The four circles in the first two columns represent the states of the sites of 
a 4-cluster; the open and filled circles correspond, respectively, to empty sites and 
sites occupied by vehicles with speed v = 1. 

Using the configurations at t and t + 1 as well as the corresponding transition 
probabilities given in table in Fig. 21 the master equation (67) for P2(l,0) 
reduces to the quadratic algebraic equation 



where we have used the shorthand notation y — P2(l, 0). Solving this quadratic 
equation we get [140] (see also [145]) 
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(76) 




(77) 
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and, hence, P2(l, 1), -P2(0, 0), P2(0, 1) from the equations (73), (74) and (75). 
Moreover, the expression (77) estabhshcs that P2(l, 0) > Pi(l)Pi(0) = c(l — c), 
which indicates an effective particlc-hole attraction (particle-particle repul- 
sion) in the NaSch model with Vmax = 1- Furthermore, from equation (77) one 
gets the expression 



J(c,p)=gP2(l,0) = - 



1 - ^Jl-Aqc{l-c) 



(78) 



which can be proved [140] to be the exact expression for the corresponding 
flux. It is not difficult to carry out 2-cluster calculations for higher values of 
Vmax, but one gets only approximate results for Vmax > 1 [140]. 
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Fig. 22. The fundamental diagram in the NaSch model for (a) v^ax = 1 and (b) 
Vmax > 1, both for p = 0.25. The data for all Vmax > 1 have been obtained through 
computer simulations. 

An interesting feature of the expression (78) is that the flux is invariant under 
charge conjugation, i.e., under the operation c — (1 — c) which interchanges 
particles and holes. Therefore, the fundamental diagram is symmetric about 
c— 1/2 when Vmax = 1 (see Fig. 22(a)). Although this symmetry breaks down 
for all Vmax > 1 (see Fig. 22(b)), the corresponding fundamental diagrams 
appear more realistic. Moreover, for given p, the magnitude of Cm decreases 
with increasing Vmax as the higher is the Vmax the longer is the effective range 
of interaction of the vehicles (sec Fig. 22). Furthermore, for Vmax = 1, flux 
merely decreases with increasing p (see eqn. (78)), but remains symmetric 
about c — 1/2 — Cm- On the other hand, for all Vmax > 1, increasing p not 
only leads to smaller flux but also lowers Cm (Fig. 23). 

For Vmax > 1; o^e uccds to carry out higher order cluster calculations [140,143] 
to get more accurate results than those obtained in the 2-cluster approxima- 
tion. For Vmax — 2, the fundamental diagrams obtained from the n-cluster 
approximation {n = 1,2, ..,5) are compared in Fig. 24 with the Monte Carlo 
data. This comparison clearly establishes a rapid convergence with increasing 
n and already for n = 4 the difference between the cluster calculation and 
MC data is extremely small. In [143] the cluster probabilities for Vmax — 2 
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Fig. 23. The fundamental diagrams in the NaSch model with Vmax = 5 are plotted 
for a few values of p. 

have been obtained from computer simulations. The results suggest that the 
n-cluster approximation for n > 3 becomes asymptotically exact in the limit 




density c 

Fig. 24. The fundamental diagrams in the NaSch model with Vmax = 2 in the 
n-cluster approximation (n = 1, 2, .., 5). 

8.3.4 Car-oriented mean-field theory of the NaSch model 

In the "car-oriented" theories the state of the traffic system is described by 
specifying the positions and speeds of all the vehicles in the system [146]. 
Suppose, Pn{t) is the probability to find at time t exactly n empty sites 
immediately in front of a vehicle. Another auxiliary quantity, which turns 
out to be very convenient to use in several different calculations, is g{t), the 
probability at time t that a vehicle moves in the next time step. These two 
sets of quantities, namely, Pn{t) and g{t) are related. For example, in the 
NaSch model with v^ax = 1, a vehicle will move in the next time step if 
there is at least one empty cell in front of it (probability J2n>iPn{t)) and 
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if it does not decelerate in the randomization step (probability q); therefore, 
g(t) = q[J2n>i Pn(t)] = q[l — Po(t)]. Hence the flux J{c,p) can be obtained 
from J{c,p) = eg = cq[l — Pq]. 

The essence of the car-oriented mean-field (COMF) approximation [146] is to 
neglect the correlations between the gaps in front of the successive car^^. 
The equations describing the time evolution of the probabilities -P„(t), under 
this approximation (see Appendix E for the derivation of these equations), 
can be solved in the steady-state using a generating function technique [146]. 
Following this approach, in the special case Vmax = 1, one recovers the exact 
expression (78) for J{c,p) = cq[l — Pq]. 

For Vmax = 2 one has to distinguish between P„(f = 1) and Pn{v = 2). 
Moreover, one has to generalize the quantity g to Qa, the probability that the 
vehicle moves a cells [a = 1, 2) in the next time step. Applying the same 
generating function techniques as for Vmax = 1, one can also solve the coupled 
sets of steady-state equations for Pn{v = 1) and P„(f = 2) for Vmax = 2 but 
gets only approximate results [146]. 

Interestingly, finite size of the system affects the equations for Vmax = 2 in 
a much more dramatic way [143] than those for Vmax = 1 thereby revealing 
the intrinsic qualitative differences in the nature of correlations in the NaSch 
model for Vmax = 1 and Vmax > 1- 

Comparisons with Monte Carlo simulations show that in contrast to the 3- 
cluster approximation for Vmax = 2 COMF does not become asymptotically 
exact in the limit p — 0. This implies that even in this limit correlations 
between the headways are not negligible. It is interesting, however, that for 
the fundamental diagram one finds an excellent agreement between MC sim- 
ulations and the predictions of COMF [146] for p 0. The reason is that 
in the deterministic limit many configurations exist which produce the same 
fiow. COMF is not able to identify the dominating structures correctly, but 
nevertheless can predict the correct current. 



8. 3. 5 Microscopic structure of the stationary state 

As we have seen MFT underestimates the fiow in the stationary state con- 
siderably. Deviations become larger for higher velocities Vmax- This shows the 
importance of correlations. As described above a particle- hole attraction exist. 
Using the 2-cluster probabilities for Vmax = 1 this attraction can be expressed 
in mathematical form as P2(l; 0) > Pi(l)Pi(0) = c(l — c). 



A similar approach, the so-called interparticle distribution functions technique, is 
used for studying reaction-diffusion systems [147]. 
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For Vmax = 1 all improvements of MFT (2-Cluster, COMF and pMFT) are 
exact. Here only correlations between neighbouring cells are important. The 
fact that pMFT is exact shows that no 'true' correlations exist. All correlations 
have their origin in the existence of GoE states. This also helps to understand 
why for random- sequential dynamics already MFT is exact for Vmax = 1 and 
the stationary state is uncorrelated. The reason is simply that for random- 
sequential dynamics no GoE states exist! 

The situation changes for higher velocities Vmax > 1- Here pMFT is no longer 
exactp^ Therefore 'true' correlations exist. This corresponds to the obser- 
vation made in [140] that the NaSch model shows a qualitatively different 
behaviour for Vmax = 1 and Vmax > 1- Furthermore it explains why so far the 
exact determination of the stationary state for Vmax > 1 has not been possible. 

It is interesting to investigate how the microscopic structure of the stationary 
state depends on the randomization p. For p = we have seen in Sec. 8.2.1 
that for densities c < l/{vmax + 1) the vehicles arrange themselves in such a 
way that all headways are at least Vmax- This is no longer possible for larger 
densities, but still the vehicles have the tendency to maximize their headway. 
Furthermore, for p = no spontaneous formation of jams exists since over- 
reactions are not possible. The behaviour in this limit can be interpreted as 
coming from a kind of "repulsive interaction" between the vehicles. 

The behaviour for p = 1 is a little bit different. Here we have seen in Sec. 8.2.2 
that metastable states with finite flow exist for c < 1/3 and Vmax > 1- 

For < p < 1 the microscopic structure interpolates between these two lim- 
iting cases. This can be seen by analysing the 3-cluster probabilities obtained 
from Monte Carlo simulations [143] . For small p the microscopic structure of 
the stationary state is determined by the 'repulsive interactions' between ve- 
hicles. With increasing p one finds a tendency towards phase separation into 
jammed and free-flow regions. A standing vehicle is able to induce a jam even 
at low densities since the restart probability is small. The jams formed are 
typically not compact, but of the form '.0.0.0. ' since a vehicle approaching the 
jam slows down in the randomization step with a rather high probability. 

Concluding one might say that the microscopic structure for < p < 1 is 
determined by the competition of the two "fixed points" p = and p = 1. 



Note that for random-sequential dynamics also MFT is no longer exact! 
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8-4 Spatio-temporal organization of vehicles; is there a phase transition? 

The density Cm corresponding to maximum flux is an obvious first candidate 
for a critical density separating the regimes of free-flow and congested flow 
in the NaSch model. We shall show in this subsection that this, indeed, is a 
critical point provided p = 0. However, in spite of strong indications that, 
probably, a noise-induced smearing of the transition takes place when p 0, 
rigorous proofs are still lacking. 

• Order parameter 

For a proper description of a phase transition one should introduce an appro- 
priate order parameter which can distinguish the two phases because of its 
different qualitative behavior within the two phases [24,25]. 

A first candidate [148] for the NaSch model would be the average fraction of 
vehicles at rest, i.e., with instantaneous speed v — O.ln the deterministic limit 
p — this, indeed, serves the purpose of the order parameter for the sharp 
transition at cf^f from the free-flowing dynamical phase to the congested dy- 
namical phase. But, in the general case of non-zero p, there is a non- vanishing 
probability that a vehicle comes to an instantaneous rest merely because of 
random braking even at extremely low density c; this probability is p for 
Vmax — 1 a-nd decreases with increasing c. 

The next obvious choice would be [149] the density of nearest-neighbor pairs 
in the stationary state 

II T L 

^ ^ t=lj=l 

where, as defined earlier, rij — for an empty cell and rij = 1 for a cell 
occupied by a vehicle (irrespective of its velocity). Because of the step 2 of 

the updating rule (deceleration due to other vehicles) m gives the spacc-time- 
averaged density of those vehicles with velocity which had to brake due to 
the next vehicle ahead. 

The Fig. 25(a) shows that m vanishes at cflf if p = 0. On the other hand if 
p ^ 0, m docs not vanish even for c < Cm although m becomes rather small 
at small densities (see Fig. 25(b)). 

We now present a heuristic argument to point out why any quantity related 

to the fraction of jammed vehicles is non-zero at any density c > and, hence, 
inadequate to serve as an order parameter [150]. To slow down to Vmax — 2 a 
vehicle must be hindered by one randomly braking vehicle in front. Similarly, 
to reach a speed Vmax — 3 a vehicle must find two vehicles within the range of 
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Fig. 25. The density-dependence of the order parameter m in (a) the deterministic 
Umit p = of the NaSch model {vmax = 1)2) and (b) the NaSch model with 
Vmax =^,P = 0.5. 

interaction, and so on. The probability for n vehicles to be found in the close 
vicinity of a given vehicle is proportional to c". Therefore, the probability 
Pv{c) of finding a vehicle with speed v < Vmax — 1 is proportional to c^^''=^~^~'" 
and, hence, even for v — 0, Pv{c) is, in general, non-zero for all c —> 0. 

• Spatial correlations 

A striking feature of second-order phase transitions is the occurrence of a 
diverging length scale at criticality and a corresponding algebraic decay of 
the correlation function [24,25]. Using lattice gas variables rij, the equal-time 
density-density correlation function is defined by 

CW^^yEEw^-c'- (80) 
^ ^ t=ij=i 

which measures the correlations in the density fluctuations that occur at the 
same time at two different points in space separated by a distance r. 

Again it is very instructive to consider first the deterministic case p — 
(Fig. 26(a)). Since, as argued before, there are exactly v„iax empty sites in 
front of each vehicle at c = c^* the correlation function at c = c^* is given by 



Qr^\^ - C^') iOT r = Omod {Vmax + l), ,g^x 

"(Cm*)^ else. 



For all Vmax the correlation function for small non-zero p (Fig. 26(b)) has es- 
sentially the same structure as that for p = (Fig. 26(a)) but the amplitude 
decays exponentially [149] for all c. In the general case of non- vanishing p, the 
asymptotic behaviour (r oo) of the correlation length C, can be obtained 
analytically [143] only for Vmax = 1- H turns out that, for given p, ^ is max- 
imum at c = 1/2 = Cm and that ^(c = 1/2) oc Thus, for v^ax = 1; C 
diverges only for p = but remains finite for all non-zero p. For Vmax > 1 the 
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Fig. 26. The density-density correlation function G{r) (a) the deterministic hmit 
p = of the NaSch model {vmax = 2) and (b) the NaSch model with Vmax = 2, 
j9= 1/128. 

trend of variation of ^ with c (Fig. 27(a)) in the vicinity of is the same 
as that for v^ax — 1 [149]. Moreover, for v^ax > 1; the maximum value of 
the correlation length, ^max plotted against p (Fig. 27(b)), is also consistent 
with the corresponding trend of variation for Vmax = 1- Thus, the correla- 
tion function G{r) gives a strong indication that the NaSch model exhibits a 
second order phase transition, at c = c^*, only for p = but this transition 
is smeared out if p 7^ 0. This noise-induced smearing of the phase transition 
in the NaSch model is very similar to the smearing of critical phenomena by 
finite-size effects. 
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Fig. 27. The dependence of (a) ^ (c, p) on c for four different values of p {vmax = 2) 
and (b) ^max on p {Vmax = 2, 3). 



• Distribution of lifetimes of jams 

Another quantity which should be able to give information about the nature 
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of the transition from free-flow to the jammed regime is the distribution of 
hfetimes of jams. Following Nagel [151] each vehicle which has a velocity less 
than Vmax before the randomization step will be considered jammed. This def- 
inition is motivated by the cruise-control limit (see Sec. 9.1.1) where it is more 
natural than in the NaSch model. One expects, however, that the long-time 
behaviour of the lifetime distribution is independent of the exact definition of 
a jam. The short-time behaviour, on the other hand, might differ strongly, e.g. 
for "compact jams" where a jam is defined as a series of consecutive standing 
vehicles without any empty cells in between. 
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Fig. 28. Lifetime distribution in the NaSch model for Vmax = 5 and p = 0.5 and 
various densities below and above Cm ~ 0.085. 

Fig. 28 shows the results of Monte Carlo simulations for the lifetime distri- 
bution in the NaSch model for different densities near the transition region, 
c Cm = 0.085 ± 0.005 (for Vmax = 5, p = 0.5), where Cm is the density where 
the flow is maximal. The most interesting feature of the lifetime distribution 
is the existence of a cutoff near = 10000. It has been shown [151] that this 
cutoff is neither a finite-size nor a finite-time effect. For times smaller than Tc 
a scaling regime exists where the distribution decays algebraically. 

• Dissolution of a megajam 

Gerwinski and Krug [150] tried to find an intuitive criterion which allows the 
distinction of free-flow and jammed phases. It is based on the investigation of 
jam dissolution times. Starting from a megajam configuration, i.e. a block of 
consecutive cells occupied by vehicles with the remaining L — N cells being 
empty, they determined the time until the jamP^ has dissolved completely. 

A simple estimate gives the density at which the lifetime is expected to become 
infinite. Suppose that the jam dissolves with velocity vj. Since the first vehicle 

^"^In [150] the same definition of a jam as in the previous point "Distribution of 
lifetimes of jams" (see [151]) has been used. 
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move freely with average velocity vp = Vmax —p it will reach the end of the jam 
at the same time as the dissolution wave if the condition (L — N)/vf — N/vj 
is satisfied. The corresponding density is then given by 



Vj + Vf 



Vj + Vmax - P 



(82) 



For Vjnax = 1 vehicles accelerate immediately to Vmax- In this case one has 
Vj = q = 1 — p. For higher velocities, q = 1 — p is only an upper bound for vj. 
Inserting vj — 1 — p into (82) one therefore obtains an upper bound for the 
density c*. Taking into account interactions between vehicles in the outflow 
region of the jam, one can derive an effective acceleration rate qeff, and thus 
the jam dissolution velocity vj — qeff, as a function of p [150]. 

Computer simulations show a sharp increase of the lifetime near the density c*. 
It becomes "infinite", i.e. the jam does not dissolve within the measurement 
time, only at a higher density c* which is considerably larger than the density 
Cm of maximum flow. At intermediate densities c* < c < cl the jam docs not 
dissolve during the flrst lap, but later due to fluctuations of the two ends of 
the jam. During this time other jams have usually formed. All results found 
in [150] are consistent with the measurements of the lifetime distributions 
presented in the previous point. 

• Relaxation time 

A characteristic feature of a second order phase transition is the divergence of 
the relaxation time at the transition point. For p = this has been studied flrst 
by Nagel and Herrmann [113]. They found a maximum Tmax oi the relaxation 
time at the density c = i / (vmax+i) which diverges in the thermodynamic limit 
L ^ oo as Tmax oc L. FoY fluite p the behaviour of the relaxation time is more 
complicated. For technical reasons Csanyi and Kertesz [152] made no direct 
measurements of the relaxation time, but used the following approach: Starting 
from a random conflguration of cars with velocity vj = the average velocity 
v{t) is measured at each time step t. For i — > oo the system reaches a stationary 
state with average velocity (^^oo)- The relaxation time is characterised by the 
parameter [152] 



v*{t) denotes the average velocity in the acceleration phase t — >■ for low ve- 
hicle density p ^ 0. In this regime, due to the absence of interactions between 
the vehicles, one has v*{t) = (1 — p)t. Thus the relaxation time is obtained 
by summing up the deviations of the average velocity {v{t)) from the values 
of a system with one single vehicle which can move without interactions with 



oo 




(83) 
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other cars (p — > 0). Note that for a purely exponentiaUy decaying quantity 
v{t) = + Cexp(— t/r') the definition (83) is proportional to r', i.e. the 
standard definition of the relaxation time. One finds a maximum of the relax- 
ation parameter near, but below, the density of maximum flow for p — 0.25 
and Vmax = 5 (see [152]). 
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Fig. 29. Relaxation parameter near the transition density for v^qx — 

5 and p = 0.25. 

The relaxation time (83) shows interesting behavior [152,139,149] which is 
difficult to interpret in terms of critical point phenomena. One finds a max- 
imum of the relaxation parameter near, but below, the density of maximal 
flow. This maximum value increases with system size, but the width of the 
transition region does not seem to shrink [139]. Two scenarios are possible: 
1) The relaxation time converges to a large but flnite value for large system 
sizes beyond the present computer power; 2) The relaxation time diverges for 
L — > oo. Scenario 1) appears to be more plausible in view of the flnite lifetimes 
of jams discussed above. Complicated interactions between jams could in prin- 
ciple lead to a divergence. Keeping in mind the unusual scaling behaviour of 
the width, this should occur in a flnite interval, not at a special (critical) point. 

The interpretation of the parameter r as a relaxation time can be problematic. 
This can be seen clearly for c > CtransiUon, where r can become negative 
[149]. Here it is possible that during relaxation the system can temporarily 
reach states with a higher average velocity than in the stationary state. This 
overreaction can be divided into two phases for p > 0. Within the flrst few 
time steps small clusters which occur in the initial conflguration vanish. The 
second phase is characterized by the growth of surviving jams. More and more 
cars get trapped into large jams and therefore the average flow decreases to its 
stationary value. This decrease causes negative values of r at large densities. 

• Distribution of distance-headways 

In order to get information on the spatial organization of the vehicles, one can 
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calculate the distance- headway distribution Vdh{^x) by following either a site- 
oriented approach [153] or a car-oriented approach [146] if Axj 



X 



i.e., if the number Axj — 1 of empty lattice sites in front of the j-th vehicle is 
identified as the corresponding distance-headway. 

Stated precisely, Vdh{k) is the conditional probability of finding a string of k 
empty sites in front of a site which is given to be occupied by a vehicle. A 
comparison between the naive mean-field expression 



(84) 



for the distance-headway distribution in the NaSch model with Vjnax — 1 and 
the corresponding Monte Carlo data [153] reveals the inadequacy of equation 
(84) at very short distances which indicates the existence of strong short-range 
correlations in the NaSch model that are neglected by the mean-field treat- 
ment. This is consistent with our earlier observation that there are particle- 
hole effective short-range attraction in the NaSch model with Vmax = 1- Again, 
this correlation disappears when a random sequential updating is carried out! 
The exact distance-headway distribution in the NaSch model with Vmax — 1 
is found to be [153,146] 



-p2c 

' dh 



y 



c(l - c) [ (1 - c) 



(J = l,2,...) 



(85) 



where y = ^2(1, 0) is given by the equation (76). 
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Fig. 30. Distributions of distance-headways in the NaSch model for p = 0.5 and 
different densities for (a) Vmax — 

1 and (b) 

Vmax — 5. n denotes the mimber of empty 
cells in front of a vehicle and is related to the distance- headway by n = Ax — 1. 

For all Vmax > 1, at moderately high densities, Vdhi^x) exhibits two peaks, in 
contrast to a single peak in the distance-headway distributions for Vm,ax — 1 
at all densities (Fig. 30); the peak at Ax = 1 is caused by the jammed vehi- 
cles while that at a larger Ax corresponds to the most probable distance- 
headway in the free-fiowing regions. At first sight, the simultaneous exis- 
tence of free- flowing and jammed regions may appear analogous to the co- 
existence of gaseous and liquid phases of matter in equilibrium. In fact, when 
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the two-peaked structure of the distance-headway distribution was first ob- 
served [153,115], it was erroneously interpreted as a manifestation of the co- 
existence of two dynamical phases, namely, the free-flowing phase and the 
jammed phase. But, later works [154] estabhshed that the analogy between 
the coexistence of free-flowing and jammed regions in the NaSch model and 
the coexistence of the gas and liquid phases of matter cannot be pushed too 
far because the analogue of the gas-liquid interfacial tension is zero in the 
NaSch model. Thus, one can not conclude that the NaSch model exhibits a 
flrst order dynamical phase transition. 

• Distributions of jam sizes and gaps between jams 

One can identify a string of k successive stopped vehicles as a jam of length 
A; (by deflnition, such jams are compact). Similarly, when there are k lattice 
sites between two successive jams, each occupied by a moving vehicle or is 
vacant then we say that there is a gap of length k between the two successive 
jams. Analytical expressions for the distributions of the jam sizes as well as of 
the gaps between jams can be calculated for the NaSch model (and some of its 
extensions) using the 2-cluster approximation or or COMF [143,153,155,156]. 
The expressions are exact in the case v^ax = 1 with periodic boundary condi- 
tions. For higher velocities the results are only approximative. In COMF the 
probability to flnd a jam of length k is given by 

whereas in the 2-cluster approach one flnds 



■J^J v=l 

For the n— cluster approximation similar expressions can be derived. 

Both distribitions (86) and (87) decay exponentially for large jam sizes. COMF 
always predicts a monotonous distribution with (j^^^^^^ > Cj^^^^^ . In con- 
trast, the jam size distribution in the n— cluster approximation can in principle 
exhibit a maximum at small jam sizes 1 < A; < n. 

• Distribution of time-headways 

Since flux is equal to the inverse of the average time-headway, much more 
detailed information is contained in the full distribution of the time-headway 
than in the fundamental diagram. The time-headway distribution contains 
information on the temporal organization of the vehicles. 
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Suppose, Vmiti) is the probability that the following vehicle takes time ti 
to reach the detector, moving from its initial position where it was located 
when the leading vehicle just left the detector site. Suppose, after reaching 
the detector site, the following vehicle waits there for t — ti time steps, either 
because of the presence of another vehicle in front of it or because of its own 
random braking; the probability for this event is denoted by Q{t — ti|ti). The 
distribution Vth{T), of the time-headway r, can be obtained from [154,157] 

T-l 

'Pthir) = E 'Pm{ti)Q{T - ti\h) (88) 
ti=i 



Substituting the expressions for Vm{ti) and Q{t — ti\ti) for Vmax = 1 in (88) 
we, finally, get [154,157] 



Qy 

c-y 



c 



t-l 



Qy ^ Qy 



c-y d-y 



+ 
P 



Qy 



oy 
d 



d-y 



(89) 



where, q = l— p^d = l — c and, for the given c and p, y can be obtained from 
equation (77). The expression is plotted in Fig. 31(a) for a few typical values of 
c for a given p. A few typical time-headway distributions in the NaSch model 
for Vmax > 1) obtained through computer simulation, are shown in Fig. 31(b). 
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Fig. 31. Time-headway distribution in the NaSch model with (a) Vmax = 1 and (b) 
• Temporal correlations 

In order to probe the spatio-temporal correlations in the fluctuations of the 
occupation of the cells, one can study the space-time correlation function 



1 1 



L T 



C{r,T) = ^-^EX^ni(t)nj+^(t + r) -c^ 



(90) 



j=i t=i 
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which, by definition, vanishes in the absence of any correlation. In [158] three 

different regimes have been distinguished. 

Free-flow (0 < p < pi): "Minijams" occur which resolve immediately. The 
correlation function shows anticorrelations around propagating peak. 
Jamming {pi < p < P2)' Free-flow and jamming coexist, i.e. jams with a flnite 
lifetime and vehicles moving with Vmax occur. This coexistence is reflected in 
the behaviour of the correlation function which exhibits a double-peak struc- 
ture. 

Superjamming (p2 < P < 1): The whole system is congested. Jamming waves 
are connected and form an inflnite wave. As a consequence, the propagating 
peak in the correlation function has disappeared. 



Neubert et al. [159] have introduced a special autocorrelation function of the 
density in order to study the velocity of jams. They have determined jam 
velocities for several variants of the NaSch model which will be introducted in 
later sections. 

• Structure factor 

Structure factors are known to give valuable information about driven systems 
[14]. For the NaSch model the static structure factor 



Hie 



ikj 



(91) 



has been investigated in [160]. Again nj denotes the occupation number of cell 
j. Note that S{k) is related to the Fourier transform of the density-density 
correlation function (80). 

For all densities, S{k) exhibits a maximum at ko ~ 0.72 which corresponds 
to the characteristic wavelength Aq = 27r/A;o of the density fluctuations in the 
free- flow regime. For Vmax > 1 one flnds ko{vmax + 1) = const. 

In [161] these investigations have been extended to the dynamical structure 
factor in velocity-space. 



Sy{k,Lj) = 



NT 



^i(kn—u>t) 



n.t 



(92) 



with k = iT^rakjN ^ lo = iT^m^jT ^ where N is the number of vehicles and 
and are integers, v^if) is the velocity of the n-th vehicle at time t. 

Compared to the dynamical structure factor in real space, (92) has the advan- 
tage that the free-flow regime only contributes white noise, S'^(A;, i^)\free-fimii = 
const. Therefore it is easier to study jamming properties. It is found in [161] 
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that Sjj{k,u!) exhibits one ridge with negative slope, corresponding to back- 
ward moving jams. One finds that the velocity of the jams is a function of 
the randomization parameter p only. It is independent of the density c and 
the maximal velocity Vmax [161]. This is consistent with results from a direct 
study of the autocorrelation function [159][^ Above a transition density, an 
algebraic behavior S^ik, uj) ~ k~'^ of the structure factor is found. This 

ui/k=Vj 

has been interpreted as an indication of critical behavior in [161]. However, 
due to the difficulties involved in the calculation of (92) only relatively short 
times T < 2048 have been considered in [161]. This is much smaller than the 
cutoff found in the lifetime of jams (see the discussion above) and lies well 
in the region where an algebraic decay is found. Therefore the results for the 
dynamical structure factor (92) and the lifetime measurements are consistent, 
but the algebraic decay is not to be interpreted as an evidence for the exis- 
tence of a critical point in the NaSch model. In order to see the cutoff, times 
T > 10'' would have to be considered. 



8.5 Exact solution of the NaSch model with v^ax = 1 one? open boundary 
conditions 



The analytical methods presented in Section 8.3 are well suited for the in- 
vestigation of translationally-invariant stationary states which are achieved 
by imposing periodic boundary conditions. For both practical and theoreti- 
cal reasons sometimes different boundary conditions are preferable. Imagine 
a situation where a multilane road is reduced to one lane, e.g. due to road 
construction. Such a bottleneck can be modeled by using a NaSch model with 
open boundaries. The multilane part of the road acts as a particle reservoir. If 
the first cell of the one-lane part is empty a car is inserted here with probability 
a. At the other end a car is removed from the last cell with probability (3 (see 
Fig. 32). These boundary conditions break the translational invariance and in 
general one can expect stationary states with a non-trivial density profile (tj). 
From a more theoretical point of view such models have been studied inten- 
sively as prototypes of systems exhibiting so-called boundary- induced phase 
transitions [162,163]. In contrast to what one expects from experience with 
equilibrium systems, one-dimensional driven nonequilibrium systems can ex- 
hibit phase transitions, even when the interactions are short-ranged, just by 
'slightly' changing the boundary conditions. 

The probability distribution characterizing the steady state of the TASEP 
with parallel dynamics (i.e., NaSch model with Vmax = 1) and open boundary 
conditions has been obtained recently in [164] and [165] using generalizations 

^^Measurements of the jam dissolution speed in [150], however, show a decrease 
with increasing Vmax and saturation for large Vmax- 
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Fig. 32. Definition of the ASEP with open boundary conditions. 



of techniques based on the Matrix Product Ansatz (MPA) (see Appendix F 
for a more technical introduction). By varying the boundary rates a and (3 
one obtains a surprisingly rich phase diagram (see Fig. 33) which is qualita- 
tively the same for all types of dynamics. Three phases can be distinguished 
by the functional dependence of the current through the system on the system 
parameters. In the low-density phase A (a < /5, ctdp)) the current is indepen- 
dent of (3. Here the current is limited by the rate a which then dominates 
the behaviour of the system. In the high-density phase B (/3 < a,/3c(p)) the 
behaviour is dominated by the output rate (3 and the current is independent 
of a. In the maximum current phase C (a > a^ip) and (3 > f3c{p)) the limiting 
factor for the current is the bulk ratef^ q = 1 — p. Here the current becomes 
independent of both a and /?. 

High- and low-density phase can be subdivided into two phases AI, AH and BI 
and BII, respectively. These subphases can be distinguished by the asymptotic 
behaviour of the density profiles at the boundaries. 

The transitions between the phases can be characterized by the behaviour of 
two correlation lengths and which only depend on p and a or /?. These 
lengths can be obtained explicitly from the exact solution. Apart from C,a and 
^13 also a third length = |,^~^ — ^^^| plays an important role. 

The transition from AH (BII) to C is continuous with diverging correlation 
length i^is)- The transition from the high- to the low-density phase is of 
first order. Here both C,a and C,i3 are finite, but ^ diverges. On the transition 
line one finds a linear density profile created by the diffusion of a domain wall 
between a low-density region at the left end of the chain and a high-density 
region at the right end. 

For the case of parallel dynamics, i.e. the NaSch model with Vmax = 1, the 



" Note that conventionally the hopping rate in the ASEP is denoted as p. Since in 
the NaSch model p is the braking probabihty the hopping rate in the ASEP (for 
Vmax = 1) becomes q = 1 - p. 
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Fig. 33. Phase diagram of the ASEP for parallel dynamics. The inserts show typical 
density profiles. 

currents in the three phases are given by [164-166] 



a{q-a) (3{q - (3) x 

Ja = Jb = Jc = 77 (1 - Vp) • (93) 

q — a'^ q — p 2 



The corresponding bulk densitiesp^ are 

a{l-a) q-p 1 

Ca = —, Cb = ^, cc = -. (94) 

g — q — p 2 

The phase boundaries are determined by the critical rates 

adp) = f3c{p) = 1 - v^. (95) 

In Fig. 33 also the special line (1 — a;)(l — /5) = p is indicated. Here the density 
profile is fiat (i.e. constant). On this line the exact solution can be obtained 
by the 2-cluster approach of Sec. 8.3.3 [166]. Since it goes through all three 
phase these results are sufficient to obtain exact analytic expressions, e.g. for 
the currents, once the structure of the phase diagram is established (e.g. by 
Monte Carlo simulations). 

The stationary state of the ASEP can also be obtained for other types of 
updates (see App. A), e.g. random-sequential [167,168], ordered-sequential 
[166,169,170] and sublattice-parallel update [171-173,166]. One finds that the 
phase diagram has the same basic structure for all updates [166]. The func- 
tional dependence of the currents, density profiles etc. on the model param- 
eters differs, however. For the important case of random- sequential updating 



In the maximum current phase no real bulk density can be defined due to the 
algebraic behaviour of the density profile, cc is therefore just 
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(93)-(95) have to be replaced by [167,168] 



Ja 



qa{l 



a) 



Pb = 1- p, P^^\' 



(96) 



PA 



(97) 



ac(p) 



1 
2 



(98) 



Results for other updates can be found in [166]. For a discussion of the calcu- 
lation of difiision constants and shock profiles we refer to the reviews [22,174] 
and references therein. 

In [162] the behaviour of the ASEP for (3 = 1 was explained by postulating 
a maximal- current principle. According to this principle, independent of the 
details of the dynamics, the system tries to maximize the stationary current 



Here J(c) is the fundamental diagram (for periodic boundary conditions) and 
c_ is the density at the left (input) boundary, i.e. c_ = a in the case described 
above. 

In [175] a nice physical picture has been developed which explains the structure 
of the phase diagram not only qualitatively, but also (at least partially) quan- 
titatively. It is determined by the dynamics of a domain wallsP^. In nonequi- 
librium systems, a domain wall is an object connecting two possible stationary 
states. The notion of domain walls in the ASEP can be illustrated in the limit 
aL <^ 1 and (3L ^ 1 of small boundary rates. At late times there will be a 
low-density region at the left end of the chain and a high-density region at 
the right end, with a domain wall in between. Schematically this state can be 
depicted as 000011111. For general values of the rates the wall not be sharp 
in general, but spread over a few lattice sites. 

For late times the dynamics of the system can then be interpreted in terms of 
the motion of the domain wall. A particle entering the system leads moves the 
wall one cell to left, and a particle leaving the system moves it one cell to the 
right. Therefore the domain wall performs a biased random walk with drift 
velocity vn = [3 — a and diffusion coefficient D = {a + (3)/2. For a < (3 the 
domain wall moves to the right until it reaches the end of the system which 
is thereafter in the low-density stationary state. For a > (3 the wall moves to 
the left until it reaches the left end and the systems goes into the high-density 

A somewhat related approach has been used to obtain an approximate solution 
for the special case of parallel dynamics with deterministic bulk dynamics {p = 0) 



J: 



J = max Jic 

ce[o,c_] 



(99) 



[176]. 
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stationary state. In the case a = (3 there is no net drift in the position of the 
walL It fluctuates with its rms displacement increasing with time as {DtY^'^, 
i.e. it can be anywhere in the system resulting in a linear density profile. 

In order to understand the case of general a and (3 one has to introduce a 
second kind of domain wall separating a maximum current phase from the 
high-density phase. Since the maximal possible flow for periodic boundary 
conditions is Jmax = 1/4 (for p = and random-sequential update) the dy- 
namics for a = 1/2 is dominated by the overfeeding at the left boundary. The 
injection rate could support a current larger than 1/4, but in the bulk it can 
not exceed this value. Therefore at the left boundary a maximum current state 
is formed. If the particles are not extracted fast enough at the right boundary 
a high-density region will develop there. These two regions are separated by 
a new kind of domain wall, the maximum current /high-density domain wall. 
Schematically it can be represented as mmmmllll. Again this domain wall 
performs a biased random walk. 

In order to obtain more quantitative predictions one goes to a coarse-grained 
picture. Then it is useful to replace the boundary rates a and P by particle 
reservoirs with densities c_ and c+. The continuity equation dc/dt+dJ/dx = 
in the continuum limit has traveling wave solutions of the form c{x — Vot) with 
the domain wall velocity 

VD = ^-^—^ (100) 

c+ — c_ 

which can be obtained by integration over the chain. 

For the low- density /high- density domain wall one hasP^ c_|_ = 1 — /3, J_|_ = 
J(c+) = (3{1 — (3) and c_ = a, J__ = J{cJ) = a{l — a). This gives indeed vn = 
(3 — a which should be valid for a,(3 < 1/2. For the maximum current/high- 
density domain wall c_ takes the value c_ = 1/2 so that J_ = 1/4 and thus 
VD = (3- 1/2. 

The arguments described above can be generalized to any process where the 
fundamental diagram J(c) of the periodic system has only one maximum at a 
density c*. For all currents J < J{c*) there exist two corresponding densities 
Ci and C2 with J(ci) = J = J {02)- For fundamental diagrams with more than 
one maximum, more than two densities might exist for a given current J. 
This implies the existence of a larger number of domain wall types. The phase 
diagram of the open system than exhibits a larger number of phases [177]. 
The maximal-current principle (99) for the TASEP with /? = 1 is generalized 



We assume p = and random-sequential dynamics. 
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to the extremal- current principle [177,15] 



maxce[c+,c_] J{.c) for c_ > c+, 
niince[c_,c+] J{c) for c_ < c+. 



(101) 



Since the above phenomenological picture does not depend on the microscopic 
details of the dynamics, it plausible that the phase diagrams for different 
updates are qualitatively the same. Boundary-induced phase transitions have 
recently been observed [178] in measurements on a German motorway [41]. 
One finds a first-order noncquilibrium phase transition between a free-flow and 
a congested phase. This transition is induced by the interplay between density 
waves induced by an on-ramp and a shock wave moving on the motorway 



9 Generalizations and extensions of the NaSch model 

As stated earlier, the NaSch model is a minimal model. The first obvious 
possible generalization would be to replace the acceleration stage of updating 
rule (Ul) to 



where, a„, acceleration assigned to the n-th vehicle, need not be unity and, 
in general, may depend on n. In the following subsections we consider more 
non-trivial generalizations and extensions of the NaSch model. 

9.1 Single-lane highways 

In the next few subsubscctions we shall demonstrate the rich variety of traffic 
phenomena that can be observed by appropriate modifications of the ran- 
dom braking. We have earlier mentioned in the context of empirical results 
that traffic flow exhibits metastability and the related hysteresis effects. Such 
phenomena have been observed in continuum formulations of "microscopic" 
models, i.e., in coupled-map lattice models [115]. However, the NaSch model 
is too simple to account for these phenomena. We now briefly describe a few 
generalizations of the NaSch model, each of which is based on modifications 
of the braking rules of the original NaSch model; one common feature of all 
of these generalized models is that they show metastability and hysteresis. 

Before we begin our discussions on specific generalized versions of the NaSch 
model, which exhibit metastability, we make some general remarks. In the 
schematic stationary fundamental diagram of Fig. 34, the low density branch 



[178]. 
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Fig. 34. Metastability and hysteresis. 



corresponds to homogeneous free-flow states, while the high density branch 
corresponds to configurations, where jammed states are present. Obviously, at 
densities Ci < c < C2, the flow depends non- uniquely on the global density. 

In order to establish the existence of meta-stable states one can follow two 
basic strategies. In the first method, the density of vehicles is changed adia- 
batically by adding or removing vehicles from the stationary state at a certain 
density. Starting in the jamming phase (large densities) and removing vehicles, 
one obtains the lower branch of the hysteresis curve. On the other hand, by 
adding vehicles to a free flowing state (low densities), one obtains the upper 
branch. This method is closely related to the experimental situation, where 
the occupancy of the road varies continuously. 

The second method does not require changing the density. Instead one starts 
from two different initial conditions, the m,ega-jam and the homogeneous state. 
The mega-jam consists of one large, compact cluster of standing vehicles. In 
the homogeneous state, vehicles are distributed periodically with same con- 
stant gap between successive vehicles (with one larger gap for incommensurate 
densities). Then, for c > Ci the homogeneous initialization leads to a free- 
flow state, while the mega-jam initialization leads to the jammed high-density 
states. 



9.1.1 Cruise- control limit and self-organized criticality 

In the cruise-control limit of the NaSch model [179] vehicles moving with their 
desired velocity v^ax are not subject to noise. This is exactly the effect of a 
cruise-control which automatically keeps the velocity constant at a desired 
value. In this model the acceleration, deceleration (due to other vehicles) and 
movement stages of updating are identical to those in the general case of 
the NaSch model; however, the randomization step is applied only to vehicles 
which have a velocity v < v^ax after step 2 of the update rule. We can express 
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this more formally by recasting the randomization stage of the update rules 
in the NaSch model as follows: 



max(0, V — 1) 



with probability 



P 



\P' 
\P 



ii V — 

if V < Vr, 



(102) 



where v is the velocity of the vehicle at the end of the step 2 of the update 
rule, i.e., after deceleration due to blocking by other vehicles. In the original 
formulation of the NaSch model Pvmax = P- ^^e other hand, the cruise- 
control limit corresponds to Pvmax 0,p ^ 0. 

For Pvmax <S 1, at sufficiently low densities, all the vehicles move deterministi- 
cally with the velocity Vmax] this deterministic motion is, however, interrupted 
by small perturbations at a vanishingly small rate. Consequently, the system 
gets enough time to relax back to the state corresponding to the deterministic 
algorithm before it is perturbed again. This effectively separates completely 
the time scales for perturbing the system and the response of the system. 




Fig. 35. Fundamental diagram in the cruise-control limit of the NaSch model for 
Vmax = 5 and = 0, P = 0.5. 

First, let us consider the periodic boundary conditions which is easier to treat 
than than the open boundary conditions. In this model, a sharp transition 
from the "free-flowing" dynamical phase to the "congested" phase takes place 
at a critical concentration c^{p, Vmax) which depends on p as well as Vmax and, 
for all p 7^ 0, c^{p,Vmax) is smaller than c^* — ^/{vmax + !)• ^ot a given 
Vmax, c*(p) increases with decreasing p and, in the deterministic limit p — > 0, 

C*(0, fmaa;) ^ ' 



In this model a jam is defined to consist of vehicles all of which have their 
instantaneous velocities smaller than Vmax- For all c < c*, jams present in 
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the initial configuration eventually disappear and in the jam-free stationary 
state every vehicle moves with the velocity v„mx- Therefore, in the density 
regime c < c* the flux increases linearly with density following J = cVmax, 
just hke that in the deterministic limit p = of the NaSch model (Fig. 35). 
But, unlike the deterministic limit p — 0, the cruise-control limit of the NaSch 
model exhibits metastability in the interval c* < c < c^*. In this context, 
the metastability means that, in the interval < c < cf^f, on appropriate 
initialization, the system can reach apparently steady states where no jam 
appears and where the fluxes are higher than J(c^.); but, perturbations of 
such a "metastable" state creates long-lived jams thereby reducing the flux 
to a level consistent with the stable branch of the fundamental diagram. At 
all c > c* jams present in the initial configuration never disappear completely 
and, in this density regime, the stable steady-state is a mixture of laminar 
flow regions and jams. The long-lived jams lower the flux beyond and the 
flux decreases linearly with density (Fig. 35). 

Let us assume that at densities slightly above c*, only one jam of length Ljam 
containing Njam vehicles exists in the system. Then, because of the periodic 
boundary conditions, the total number of vehicles N is conserved and, hence, 

N — CjamL jam + Cout{L — Ljam) (103) 

where Cjam = Njam/ Ljam and Cout = {N - Njam)/ {L - Ljam) are the densities 
of the vehicles in the jam and in the outflow region, respectively. Dividing 
both sides of (103) by L we get 

Ljam , [ -, Ljam\ f-\r\A\ 
C ^ Cjam^^ + Cout ] (104) 

Since in the cruise-control limit of the NaSch model Ljam must vanish as 
c — > c*, we conclude [179] that we must have Cout = c*, i.e., the average density 
in the outflow region of a jam is equal to the critical density c*. 

In order to study the traffic at the critical point of the cruise-control limit, 
Nagel and Paczuski [179] used a a special boundary condition which enables 
the system to select automatically the state of maximum throughput, i.e, 
the system exhibits self- organized criticality. This special boundary condition 
consists of an infinite jam from — oo to (i.e., at the left boundary) while the 
right boundary is open. Vehicles emerge from the infinite jam in a jerky way, 
before attaining the velocity Vmax- Far away from the infinite jam all vehicles 
move with the same velocity Vmax- In order to show that the state selected this 
way is "critical" [179] we perturb a vehicle, far downstream from the infinite 
jam, slightly by reducing its velocity from Vmax to Vmax — 1- This particular 
vehicle initiates a chain reaction and gives rise to a jam if the following vehicle 
is sufficiently close to it although it itself accelerates and, eventually, attains 
Vmax- This phantom jam has a time-dependent size n{t), measured by the 
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number of vehicles n in this jam at time t and it has a hfetime rufe- The 
statistics of this features of the phantom jam can be obtained by repeating 
the computer experiment sufficiently large number of times; sometimes the 
phantom jam created is small and has a short lifetime and sometimes it is large 
and has quite long lifetime. Interestingly, the characteristic quantities like, for 
example, the distributions of the sizes of the jams, lifetimes of the jams, etc. do, 
indeed, exhibit power-laws which are hall mark of the self-organized criticality 
[26,27]. E.g. the branching behaviour of the jams gives rise to intermittent 
dynamics with a 1// power law spectrum [179]. 1// noise in real traffic has 
been discovered by Musha and Higuchi [180]. They recorded transit times 
of vehicles passing underneath a bridge. The corresponding power spectral 
density of the flow fluctuations shows 1// behaviour at low frequencies. 

The exponents associated with the various power laws in the cruise-control 
limit of the NaSch model can be calculated analytically, at least for v^ax = 1? 
by utilizing a formal relation with one-dimensional unbiased random walk. If 
Vmax = Ij all the vehicles in the jams have velocity v = 0. Moreover, the jams 
are compact so that the number of vehicles in a jam is identical to its spatial 
extent. The probability distribution P(n; t) for the number of vehicles n in 
such a jam, at time t, is determined by the following equation: 



P{n; t+l) = {l-rin- rout)P{n; t) + rinP{n - l;t) + routP{n + 1; t) (105) 



where the phenomenological parameters and Tout are the rates of incoming 
and outgoing vehicles. Of course, rj„ depends on the density of the vehicles 
behind the jam. For large n and t, taking the continuum limit of the equation 
(105) and expanding we get 



If Tin > Tout the jams would grow for ever. On the other hand, the jams would 
shrink, and eventually disappear, if rj„ < Tout- If f^in = f^out, the first term on 
the right hand side of the equation (106) vanishes and the resulting equation 
governing the time evolution of P(n; t) is identical to that of the probabihty of 
finding, at time t, an unbiased one-dimensional random walker at a distance 
n which was initially at the origin. Thus, when rj„ = Tout, the jams exhibit 
large ("critical") fiuctuations which can be characterized by critical exponents. 
Using this formal mapping onto unbiased random walk, we find (a) that the 
mean size of jam at time t corresponds to the mean displacement of the random 
walker from the origin after time interval t and (b) that the lifetime of a jam 
corresponds to the time taken by the random walker to return to the origin for 
the first time. Hence, using the well known results from the theory of random 
walks [181,182], we get 




dP 



dP , Tin + Tout d'^P 



(106) 



n(t) oc t^^^, and P(riif(.) oc 



3/2 



(107) 
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It turns out that the power-law exhibited by the size of the jams, the dis- 
tributions of the hfetimes, etc. are not restricted merely to the special case 
Vmax = 1 of the cruise-control limit but is also shown by the corresponding 
computer simulation data also for arbitrary Vmax- The power- law distributions 
of P{Tiife) in the cruise-control limit of the NaSch model is in sharp contrast 
with the exponential distribution observed in the NaSch model [151]. Thus, 
in the cruise-control limit of the NaSch model, the large jams are fractal [183] 
in the sense that there are smaller sub-jams inside larger jams, ad infinitum. 
In other words, in between sub-jams, there are holes of all sizes. 

9.1.2 Slow-to-start rules, metastability and hysteresis 

The slow-to-start rules can lead not only to metastability and, consequently, 
hysteresis, but also to phase separated states at high densities, as we now 
show. 

• Takayasu-Takayasu slow-to-start rule 

Takayasu and Takayasu (TT) [184] were the first to suggest a CA model with 
a slow-to-start rule. Here, we investigate the generalization, as suggested in 
[185], of the original slow-to-start rule. According to this generahzed version, 
a standing vehicle (i.e., a vehicle with the instantaneous velocity v — 0) with 
exactly one empty cell in front accelerates with probability qt = 1 — pt, while 
all other vehicles accelerate deterministically. The other steps of the update 
rule (U2-U4) of the NaSch model are kept unchanged. 

As in the case of the NaSch model, it is instructive to consider first the de- 
terministic limits of the TT model [184,186]. The TT model reduces to the 
NaSch model in the limit pt = 0. What happens in the other determinis- 
tic limit, namely, pt = 1? In the latter deterministic limit, a stopped vehicle 
can move only if there are at least two empty cells in front [184]. Obviously, 
completely blocked states exist for densities c > 0.5, where the number of 
empty cells in front of each vehicle is smaller than two. However, in the region 
0.5 < c < 0.66 the number of blocked configurations is very small compared to 
the total number of configurations and states with a finite fiow exist. Precisely 
at c = 0.5, there are only two blocked states and the time to reach these states 
diverges exponentially with the system size. 

The fundamental diagram for the TT model with Vmax — 1 has been derived 
analytically by carrying out (approximate) 2-cluster calculations in the site- 
oriented approach [185]. But, the fundamental diagrams of the TT model 
for all Vmax > 1 have been obtained so far only numerically by carrying out 
computer simulations (see Fig. 36). 

Comparing these fundamental diagrams with the corresponding ones for the 
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Fig. 36. The fundamental diagram in the TT model {vmax = 5, 
p = 0.01, pt = 0.75, L = 10000) obtained using two different initial conditions, 
namely, a completely jammed state (jam) and a homogeneous state (hom) and 
averaging over 100000 sweeps through the lattice. 

NaSch model {pt = 0), we find the following effects of the TT slow-to-start 
rule: (i) for a given density c, the flux J(c) is smaller in the TT model as 
compared to that in the NaSch model; (ii) the particle-hole symmetry is not 
exhibited by the TT model for any Vmax (not even for Vmax = 1) and (iii) 
the TT model exhibits metastability and hysteresis which are absent in the 
NaSch model. Note that the mechanism for meta-stability in the case pt = 1 
is different from that for the metastability observed for < pt < 1. 

Because of the slow-to-start rules, the separations between the vehicles coming 
out of a jam are larger than those between the vehicles coming out a jam in 
the NaSch model. Since the density far downstream is smaller than the density 
of maximum flow, the vehicles can propagate freely in the low density regions 
of the lattice where spontaneous formation of jams is highly unlikely, if the 
parameter p is sufficiently small. Therefore, the phase-separated steady-states 
at high global densities consist of a macroscopic jam and a macroscopic free- 
flow regime both of which coexist simultaneously (Fig. 37). 

• The BJH model of slow-to-start rule 

Benjamin, Johnson and Hui (BJH) [187] modifled the updating rules of the 
NaSch model by introducing an extra step where their "slow-to-start" rule is 
implemented; this slow-to-start rule is different from that introduced by TT 
[184]. According the BJH "slow-to-start" rule , the vehicles which had to brake 
due to the next vehicle ahead will move on the next opportunity only with 
probability 1 — Ps- The steps of the update rules can be stated as follows: 

Step 1: Acceleration. Vn — > min(v„ -|- l^Vmax)- 

Step 2: Slow-to-start rule: If flag — 1, then — > with probability Ps- 
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Fig. 37. Typical space-time diagram of the TT model with Vmax — 5 and c = 0.20, 
p = 0.01 and pt = 0.75. Each horizontal row of dots represents the instantaneous 
positions of the vehicles moving towards right while the successive rows of dots 
represent the positions of the same vehicles at the successive time steps. 

Step 3: Blockage (due to other vehicles). Vn — > mm{vn, dn — 1) and, then, 
flag = 1 if = 0, else flag — 0. 

Step 4- Randomization. Vn — >■ max(t'„ — 1,0) with probability p. 



Step 5: Vehicle movement. x„ — > x„ + v,, 



Here flag is a label distinguishing vehicles which have to obey the slow-to- 
start rule {flag — 1) from those which do not have to {flag — 0). 

Obviously, for = the above rules reduce to those of the NaSch model. 
The slow-to-start rule of the TT model is a 'spatial' rule. In contrast, the 
BJH slow-to-start rule requires 'memory', i.e. it is a 'temporal' rule depending 
on the number of trials and not on the free space available in front of the 
vehicle. The fundamental diagram of the BJH model (Fig. 38) clearly shows 
the existence of mctastablc states and, therefore, expected to exhibit hysteresis 
effects [155]. But, in the special case of Vmax = 1, for which approximate 
analytical calculations can be carried out [185], no meta-stable states exist. 

Since for all v^ax > 1 in the BJH model, just as in the TT model, the outflow 
from a jam is smaller than the maximal flow, the phase-separated steady- 
states at high global densities consist of a macroscopic jam and a macroscopic 
free- flow regime both of which coexist simultaneously (Fig. 39) [155]. 
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Fig. 38. The fundamental diagram in the BJH model [vmax = 5, p = 0.01, = 0.75) 
obtained using two different initial conditions, namely, a completely jammed state 
(jam) and a homogeneous state (hom). 
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Fig. 39. Typical space-time diagram of the BJH model with Vmax = 5 and c = 0.20, 
p = 0.01 and ps = 0.75. Each horizontal row of dots represents the instantaneous 
positions of the vehicles moving towards right while the successive rows of dots 
represent the positions of the same vehicles at the successive time steps. 

However, the macroscopic jam is not compact. The typical size of the macro- 
scopic free-flow regime can be estimated by measuring the distribution of the 
gaps between the successive jams [155]. A peak occurs in this distribution for 
headways of the order of the system size (see the inset of the right part of 
Fig. 40). The position of the peak indicates the typical size of the macroscopic 
free-flow regime. 

• The NaSch model with a velocity-dependent slow-to-start rule 

Although the NaSch model does not exhibit metastable states and hysteresis, 
a simple generalization exists which is able to reproduce these effects. It is 
the so-called Velocity-Dependent-Randomization (VDR) model [188]. Here, in 
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Fig. 40. Distribution of gaps between jams in the BJH and the NaSch model for 
Vmax = 5, p = 0.05, ps = 0.5 and density c = 0.5. 

contrast to the original NaSch model, the randomization parameter depends 
on the velocity of the vehicle, p = p{v) . The rules (see Sec. 8) are supplemented 
by a new rule. 

Step 0: Determination of the randomization parameter. The randomization 
parameter used in step 3 for the n-th vehicle is given hy p = p{vn{t)). 

This new step has to be carried out before the acceleration step 1. The ran- 
domization parameter used in step 3 depends on the velocity Vn{t) of the n— th 
vehicle after the previous timestep. In order to implement a simple slow-to- 
start rule one chooses [188] 



with po > P- This means that vehicles which have been standing in the previous 
timestep have a higher probability po of braking in the randomization step than 
moving vehicles. 

The rules of the VDR model can be recast in a form similar to those of the 
BJH model. We define a label flag which distinguishes between vehicles which 

have to obey the slow-to-start rule {flag — 1) from those which do not have to 
{flag = 0). flag = 1 if f „ = at the beginning of a time step, else flag = 0. 
Explicitly, the update rules are as follows: 

Step 1: Acceleration. Vn — > min(v„ + l^Vmax)- 

Step 2: Blockage (due to other vehicles). Vn — > min(^;„, d„ — 1), 




for V = 0, 
for V > 0, 



(108) 
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Step 3: Randomization. t'„ — max(f„ — 1, 0) with probability po if flag — 1. 
else, Vn — > max(i>„ — 1,0) with probability p. 

Step 4-' Vehicle movement. x„ — > x„ + 

Let us compare this VDR model with the cruise-control limit of the NaSch 
model. The vehicles with velocity v = v^ax (at the end of the step 2) are 
treated deterministically in the cruise-control limit whereas in the VDR model 
velocities of all those with the velocity w > (just before the step 3) are 
updated stochastically, but using different values of the braking parameter. 

Typical fundamental diagrams look like the one shown in Fig. 41 where, over 
a certain interval of c, J(c) can take one of the two values depending on the 
initial state and, therefore, exhibit metastability. 
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Fig. 41. The fundamental diagram in the NaSch model with a velocity-dependent 
slow-to-start rule {vmax = 5, = 0.75, p = 1/64) obtained using two different 
initial conditions, namely, a completely jammed state (jam) and a homogeneous 
state (horn). 

Moreover, typical space-time diagrams of the VDR model (see Fig. 42) clearly 
demonstrate that mctastablc homogeneous states have a lifetime after which 
their decay leads to a phase separated steady state. The microscopic structure 
of these phase-separated high-density states is quahtative similar to those 
observed in the high-density regimes of the TT and BJH models but differs 
drastically from those found in the NaSch model. 

It is instructive to compare the fundamental diagram of the VDR model with 
those of the corresponding NaSch models. We now present a simple derivation 
of the fundamental diagram of the VDR model on the basis of heuristic argu- 
ments utilizing the observed structures of the steady- states. For small densities 
c <^ 1 there are no slow vehicles in the VDR model since interactions between 
vehicles are extremely rare. In this regime every vehicle can move with the 
free-flow velocity Vf ^ (I- p)vmax + p{vmax - 1) = Vmax - P and, therefore. 
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Fig. 42. Typical space-time diag ram of the VDR model with Vjnax — 

5 and c = 0.20, 

p = 0.01 and pQ = 0.75. Each horizontal row of dots represents the instantaneous 
positions of the vehicles moving towards right while the successive rows of dots 
represent the positions of the same vehicles at the successive time steps. 

the flux is given by 

Jhom{c) = c{Vmax " P) (109) 

which is identical to the NaSch model with randomization p. On the other 
hand, for densities close to c = 1, the vehicles are likely to have velocities v = 
or V = 1 only and, therefore, the random braking is dominated by Pq, rather 
than p, while the flow is determined by the movement of the holes. Hence, for 
large densities, i.e., l — c<^ 1, the flow is given by J(c) ^ (1 —po){l — c) which 
corresponds to the NaSch model with randomization pq. This expression for 
flux in the high-density regime can also be derived as follows. In the phase- 
separated state the vehicles are expected to move with the velocity Vf = 
Vmax — p in the free-flow region. Neglecting interactions between vehicles in 
the free-flowing region (which is justifled because of the corresponding low 
density) , the average distance of two consecutive vehicles in the free- flow region 
is given by Ax = cj^ = TyjVf + l where the average waiting time of the first 
vehicle at the head of the megajam is given by = j^^- In other words, 
the density in the free-flow regime c/ is determined by the average waiting 
time Tyj and Vf. Now suppose that Nj and Np are the number of vehicles in 
the megajam and free-flowing regions, respectively. Using the normalization 
L = Nj + NpAx we find that for the density c = ^f+Nj ^ ^-^^ Jsep{c) is 
given by Jsep{c) = ^{vmax -p) and, hence, 

Jser,{c) = {l-po){l-c). (110) 

Obviously, c/ is precisely the lower branching density Ci, because for densities 
below c/ the jam-length is zero. It should be noted that the heuristic arguments 
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presented above remain valid for po ^ P and Vmax > 1- The condition p <^ 1 
guarantees that the jams are compact in that hmit. In the case v^ax = 1, 
vehicles can stop spontaneously, even in the free-flow regime and these vehicles 
might initiate a jam. This is the basic reason why hysteresis is usually not 
observed for Vmnr — 1- 
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Fig. 43. Peaks of the jam-gap distribution at macroscopic distances. 



Analogous to the BJH model phase separation can be directly identified using 
the results of the jam-gap distribution. Fig. 43 shows that the size of the 
free-flow regime is proportional to the system size. 

The results for the slow-to-start models discussed above have been obtained 
by computer simulations of periodic systems of finite length. It was shown that 
the fundamental diagram which is sketched in Fig. 34 is generic for all models 
under consideration. Now it is self-evident to ask what kind of stationary 
states are realised in the thermodynamic limit L — > oo. The simulation results 
indicate that Ac = Ci — C2 decreases with larger system sizes and is expected 
to vanish for L 00, i.e. the jammed branch is stable in that limit. This is 
readily understood if one analyses the typical configurations which lead to an 
emerging jam or, vice versa, the mechanism of the dissolution of a jam. Jams 
emerge if overreactions of drivers lead to a chain reaction. This is possible in 
dense regions of the free-fiow state where the gap between the vehicles is not 
larger than Vmax- Obviously the probability to find large platoons of vehicles 
driving with small spatial headways increases with the system size (for fixed 
density). In addition to that the jammed states are phase separated, i.e. the 
size of the jam is of the order of the system size. During a simulation run the 
size of the jam fiuctuates due to the stochastic movement and acceleration of 
the vehicles. Jams can dissolve if the amplitude of these fiuctuations are of 
the order of the length of the jam, which is impossible in the thermodynamic 
limit. 
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Therefore the non-unique behaviour of the fundamental diagram is only ob- 
servable if finite system sizes are considered or if the vehicles move determin- 
istically in the free flow regime. Nevertheless the results discussed above are 
highly relevant for practical purposes, because the hysteresis effects have been 
observed at reahstic system sizes (e.g. L — 10000 corresponds to a highway of 
length 75 km). 

9.1.3 Flow- optimization and meta-stable states 

Hysteresis effects and meta-stable states are not only of theoretical interest, 
but also have interesting applications. From the previous discussion of the 
slow-to-start models it is evident, that one can optimize the maximum flow, 
if the homogeneous state is stabilized by controlling the density so that it 
never exceeds C2. This strategy was followed in minimizing frequent jams in 
the Lincoln- and the Holland-Tunnels in New York. Before the traffic lights 
were installed at the entrance of the tunnels jams used to form spontaneously 
within the tunnel because (a) the vehicle density used to be sufficiently high 
and (b) the drivers used to drive more carefully inside the tunnel thereby 
giving rise to stronger fluctuations which caused the jams. But, the traffic 
hghts installed at the entrance of the tunnels do not allow the density to 
exceed C2 and, consequently, jams are not formed spontaneously by the decay 
of any metastable high-density state. 



Normal part of the road Bottleneck 




traffic light 



Fig. 44. A schematic representation of a tunnel on a highway. Larger values of the 
braking parameters are used for the right part, which mimics the tunnel. 

One can mimic the situation of Lincoln and Holland-tunnels within the frame- 
work of the CA models in the following way [188-190] . The tunnel is considered 
as part of the road, where the braking probabilities p*,Po higher compared 
to the remaining part of the lattice {p,po, see Fig. 44). Therefore, if one allows 
for an uncontrolled inflow of the vehicles, jams typically appear inside the 
"tunnel" and the system capacity is governed by Pq. 

The situation differs drastically if traffic lights are implemented [188-190]. 
As shown in Fig. 45, a considerable increase of the maximum capacity can 
be achieved for an optimal combination of the red-/green-signal periods. The 
gain of capacity obtained for the optimal intervals of the signal is of the same 
order as for the reahstic examples [191]. 
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Fig. 45. Average flow for different combinations of red-/green-signal periods for a 
system of lengtli L = 1000 and Vmax = 5. Inside the tunnel = 0.15, Pq = 0.60 and 
outside p = 0.01, pq = 0.20 has been used. 

9.2 Multi-lane highways 

For a realistic description of traffic on highways the idealized single-lane mod- 
els must be generalized to develop CA models of multi-lane traffic; the main 
ingredient required for this generalization being the lane-changing rules. Sev- 
eral attempts have been made so far in this direction [192-199]. The lane 
changing rules for two-lane traffic can be symmetric or asymmetric with re- 
spect to the lanes. Similarly, if there are two (or more) different types of 
vehicles (say, cars and trucks) with two different fmax, the lane-changing rule 
can be symmetric or asymmetric with respect to the vehicles. 

In general, the update in the two-lane models is divided into two sub-steps: 
in one sub-step, the vehicles may change lanes in parallel following the lane- 
changing rules and in the other sub-step each vehicle may move forward effec- 
tively as in the single-lane NaSch model. Drivers must find some incentive in 
changing the lane. Two obvious incentives are (a) the situation on the other 
lane is more convenient for driving, and (b) the need to make a turn in near 
future. Two general prerequisites have to be fulffiled in order to initiate a lane 
change: ffist, there must be an incentive and second, the safety rules must 
be fulffiled [200]. Lane changing rules according to this scheme have been in- 
troduced by Rickert et al. [193]. They suggested that vehicles are allowed to 
change the lane if the following four criteria are satisfied: 

(CI) gap{i) < /, 
(C2) gapo{i) > L, 
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(C3) gaPo,back{'i) > lo,back, 

(C4) randO < Pc- 

Here gap{i) and gapo{i) are the gaps in front of vehicle i on the own lane and 
the other lanep'^ respectively, gapo^hackij) is the gap on the other lane to the 
next vehicle behind. /, /o, lo,back and pc are parameters specifying the rule and 
rand{) is a random number in the interval [0, 1]. 

The first rule CI represents the incentive criterion, i.e. if the gap gap{i) in 
front of the vehicle is not sufficiently large vehicles want to change the lane. 
Typical choices of the parameter I are given by / = min(f +1, Vmax)- This choice 
of the minimal headway ensures that vehicles driving in a slow platoon try to 
change the lane if possible. In the next rule C2 it is checked if the situation 
on the other lane is indeed more convenient. This motivates the choice I = lo- 
The third rule C3 avoids too small distances to following vehicles on the other 
lane. Rickert and coworkers suggested lo,hack = Vmax- It is also important to 
perform lane changing stochastically. Even if the incentive and safety criteria 
are fulfilled a lane change is performed only with probability pc (C4). This 
avoids, at least partially, so called ping-pong lane changes, i.e. multiple lane- 
changes of vehicles in consecutive timestepsp^. Already implementations of 
the NaSch model using the basic lane-changing rules revealed quite realistic 
results. Nevertheless several variants of the basic rules have been developed in 
order to improve the realism. A large number of lane changing rules considered 
in the literature have been tabulated and compared by Nagel et al. [196] (see 
Fig. 46). 

The lane changing rules for two-lane traffic can be symmetric or asymmetric 
with respect to the lanes [193] . If symmetric lane changing rules are applied the 
rules do not depend on the direction of the lane changing maneuver. In contrast 
also asymmetric lane changing rules have been considered. Lane changing rules 
can be asymmetric in two ways. First it is possible that it is preferred to drive 
on the right lane at low densities. This behavior can be implemented simply 
by leaving out the first rule for a change from the left to the right lane. Second 
it is also possible that it is even forbidden to overtake a vehicle on the right 
lane, e.g. on german highways. Then the single lane dynamics on the right 
lane depends on the configuration on the left lane. These examples show the 
fiexibility of the CA approaches. Moreover the simulations also show that the 
details of the lane changing rules may lead to considerable changes of the 
model results [194,196]. 

In multi-lane traffic it is of particular interest to investigate systems with differ- 

"^^ The gap on the other lane is defined in the same way as the gap on the own lane 
by imagining that the vehicle occupies the site parallel to its current position. 

This artifact of the parallel update was already pointed out by Nagatani [192], 
who simulated a two-lane system with Vmax = 1- 
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Fig. 46. The table gives an overview over different choices of the lane changing rules 
discussed in the literature. The numbers of the references correspond to those in 
[196]. The first column shows different incentive criteria and the second column 
denotes the corresponding minimal gap on the other lane. In general the incentive 
criteria can be different for a change from left to right {L ^ R) and {R L) (from 
[196]). 

ent types of vehicles. For CA models this has been done first by Chowdhury et 
al. [195] who simulated a periodic two-lane system with slow and fast vehicles, 
i.e. vehicles with different Vmax- The simulation results have been shown that 
already for small densities the fast vehicles take on the average the free-flow 
velocity of the slow vehicles, even if only a small fraction of slow vehicles have 
been considered. Analogous results have been obtained by Helbing and Huber- 
man [212] who used a different CA model for the in-lane update (see Sec. 11.3 
for the deflnition of the model). In addition to that Nagel et al. [196] have been 
shown that for a suitable choice of the lane changing rules and different types 
of vehicles even the phenomenon of "lane inversion" which has been observed 
at german highways can be reproducedP^. The results discussed so far show 

On german highways the left lane is considered for overtaking vehicles only. There- 
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the strong influence of slow vehicles in multi-lane systems. They fit fairly well 
the empirical results, which show an alignment of the speeds on different lanes 
and and of different types of vehicles. Nevertheless recent simulation results 
of Knospe et al. [197] indicate that the influence of slow vehicles seems to be 
overestimated by the multi-lane variants of the NaSch model. In particular 
for symmetric lane changing rules even a single slow vehicle can dominate the 
dynamics close to the optimal value of the flow. In order to weaken the effect 
of slow vehicles they suggested to consider anticipation effects, i.e. the driver 
estimates the velocity of the vehicle in the next timestep [197]. 

Another interesting quantity to look at is the frequency of lane changes at 
different densities. Here the simulation results show that close to the density 
of maximal flow the number of lane changing maneuvers drastically decreases 
if the small values of the braking noise are considered in CA models where 
the velocities of vehicles are solely determined by the distance to the vehicle 
ahead. This is due to the fact that for homogeneous states at high densities no 
sufficiently large gaps exist. For larger values of the braking noise large density 
fluctuations are observable. Therefore the local minimum of the lane-changing 
frequency is not found for larger values of p. 

In general the simulation results show that some generic multi-lane effects can 
be pointed out. First of all the maximal performance of multi-lane systems is 
slightly increased compared to corresponding single-lane network. In addition, 
slow vehicles lead to an ahgnment of velocities of different type of vehicles 
already at low densities which is confirmed by empirical observation. This 
effect is quite robust for different choices of the CA model as well as for 
different lane changing rules. It can be weakened most efficiently if anticipation 
effects are applied. The details of the lane changing rules, however, may have 
strong influence on the lane usage characteristics. 



9.3 Bidirectional traffic 

Simon and Gutowitz [201] have introduced a two-lane CA model where the 
vehicles move in opposite directions. Passing may be allowed on one or on 

both lanes. It is only attempted if there is a chance to complete the pass. 
Therefore drivers measure the local density, i.e. the density of vehicles in front 
that have to be passed. If it is sufficiently low, a pass will be attempted. This 
means that at high global densities the lanes are effectively decoupled since 
only very few passes will occur. 



fore, at low densities, the right lane is used more often. Surprisingly, at higher den- 
sities not simply a balancing of the lane usage has been observed, but for densities 
close to the optimum flow the left lane is even higher frequented. 
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In principle, three types of jams can occur on a bidirectional road: 1) Sponta- 
neous jamming and start-and-stop waves on one of the lanes; 2) jams caused 
by drivers who try to pass but can not return to their home lane since there is 
not enough space and 3) "super jams" when an adjacent pair of drivers tries 
to pass simultaneously. These super jams halt traffic on both lanes and can 
be prevented by breaking the symmetry between the lanes. 

The precise rules of the CA are in the same spirit as the rules for multilane 
traffic described in the previous subsection 9.2. First, the situation on the own 
lane is examined. If the motion is hindered by another vehicle (moving in the 
same direction), a pass is attempted. This will only be initiated if the safety 
criteria are satisfied: 1) The gap on the other lane has to be sufficiently large, 
and 2) the number of vehicles to be passed has to be small. Even if these 
criteria are satisfied a lane change occurs only with probability p change- After 
this lane changing step the vehicles move forward similar to the dynamics of 
the NaSch model. There are, however, important differences. Passing vehicles 
never decelerate randomly. In order to break the symmetry between the two 
lanes moving vehicles which are on their home lane and see oncoming traffic 
decelerate deterministically by one unit. This rule prevents the occurance of 
a super jam. 

The results of [201] show the expected behaviour, namely that passing makes 
traffic more fluid. Start-stop waves are surpressed if the density is not too 
large. The improvement of the flow on one lane compared to the the one-lane 
model depends on the density of vehicles on the other lane. It is maximal for 
very small densities (c 0) on the passing lane. If the density on the other 
lane is small (c < 0.25) the flow may be lower than in the one-lane model since 
passing oncoming vehicles create an additional hindrance. For large densities 
on at least on of the lanes there is little difference between the one- and two- 
lane models. 

Lee et al. [202] have proposed a toy model for bidirectional traffic based on a 
multispecies generalization of the ASEP. Here no passing is allowed. Instead 
oncoming traffic on the opposite lane reduces the hopping rates of the vehi- 
cles. The dynamics on each lane is given by that of the ASEP with random- 
sequential update and v^nax — but the hopping rate from an occupied cell 
j to an empty cell j + 1 on lane 1 depends on the occupancy of cell j + 1 on 
the opposite lane (lane 2). When this cell is empty, vehicles hop with rate 1, 
otherwise with rate 1/13. On lane 2 vehicles move in the opposite direction and 
the hopping rate from cell j -|- 1 to cell j depends on the occupancy of cell j 
on lane 1. It is given by 7 when this cell is empty and by 7//3 if it is occupied. 

For 7 < 1 the uninffuenced hopping rate on lane 2 is smaller than that of lane 1. 
The vehicles on lane 2 might therefore be interpreted as trucks. The interlane 
interaction parameter (5 can be interpreted as a kind of road narrowness. For 
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(3 = 1 vehicles arc not slowed down by oncoming traffic. This corresponds to 
a highway with divider. The case [3^0 corresponds to a narrow road being 
completely blocked by the oncoming traffic. 

The behaviour of the model with only one truck is rather similar to that of 

the NaSch model with quenched disorder (see Sec. 10). For (3 > (3c the system 
segregates into two phases, a high-density phase in front of the truck and a 
low-density phase behind it. 

By forbidding trucks and cars to occupy parallel cell j simultaneously the 
model can be mapped onto an exactly solvable 2-species variant of the ASEP. 
Using the matrix-product Ansatz (see App. F) many steady-state properties 
for the single-truck case can be obtained exactly. Two phase can be distin- 
guished: A low-density phase for c(3 < 1 and a jammed phase for c(3 > 1 
where c is the density of vehicles on lane 1. In contrast to the case of a fixed 
defect site (see Sec. 10.3) only one critical density Carit = 1/(3 exists since the 
particle-hole symmetry is broken. 

Generalizations of this model to other updates and higher velocities can be 
found in [203]. 



10 Effects of quenched disorder on traffic 

10.1 Randomness in the braking probability of drivers and Bose-Einstein-like 
condensation 

We have seen how modifications of the random braking probability or the 
rule(s) for random braking in the NaSch model can give rise to a rich variety 
of physical phenomena, e.g., self-organized criticality, metastability and hys- 
teresis, etc. Now we consider the effects of quenched randomness in the random 
braking probability p, i.e., we study the effects of assigning randomly different 
time-independent braking probabilities pi to different drivers i in the NaSch 
model. Such "quenched" (i.e., time- independent) randomness in the random 
braking in the NaSch model can lead to exotic phenomena [204,197] which 
are reminiscent of "Bosc-Einstein-like condensation" in the TASEP where 
particle-hopping rates are quenched random variables [205,206]. Various as- 
pects of these phenomena have been thoroughly reviewed by Krug [207] and, 
therefore, we'll restrict our discussion to only the essential points. 

Let us first consider the special case of the NaSch model with Vmax = 1. As 
explained earlier, this model reduces to the TASEP if the parallel updating 
is replaced by random sequential updating scheme. If the same hopping prob- 
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ability q is assigned to every particle except one for which the the hopping 
probability is q' < q, then the single " impurity" particle is the slowest moving 
one. The faster particles can be allowed to overtake the slow one at a non-zero 
rate [174,208]; however, if this rate of overtaking vanishes the slow particle will 
give rise to a platoon of particles behind it. This phenomenon is very similar 
to the formation of platoons of vehicles in a traffic behind the slow vehicles 
(e.g., trucks). 

Here we are interested in a more general situation of quenched "disorder" in 
the form of a distribution of intrinsic hopping probabilities of the vehicles in 
the system rather than that of the single "defect" particle. In such situations 
random initial conditions can lead to the formation of platoons if (a) slow 
particles are sufficiently rare and (b) if the density of vehicles is sufficiently 
low. Following their formation, starting from a random initial condition, the 
platoons grow through coalescence. The coarsening of the platoons has been 
investigated in the same spirit in which coarsening of domains (the so-called 
Oswald ripening) is monitored while studying spinodal decomposition in, for 
example, binary alloys [28]. Suppose, ^{t) is the typical platoon size at time t. 
Starting from a homogeneous spatial distribution of the vehicles, ^{t) can be 
monitored as a function of time t to find out the law of "growth" of the size 
of the platoons. 

Before describing the effects of the quenched randomness in the hopping prob- 
abihties on the steady-states of the TASEP and the NaSch model, we consider 
an even simpler model of platoon formation [209,210] which was developed us- 
ing the language of aggregation phenomena. In this model an initial velocity 
Vj is assigned to each vehicle j, drawn randomly from a continuous probability 
density f{v). The particles then move ballistically along a line and coalesce 
whenever a faster vehicle catches up with a slower one in front. It has been 
found that ^{t) increases indefinitely according to the power law 



where the exponent n characterizes the behaviour of f{v) in the vicinity of 
the minimal velocity Vmin, i-e., f{v) ~ A{v — VminY as v ^ Vmin with some 
positive constant A. An attempt has been made to develop a coarse-grained 
description of this phenomenon [211]. 

It has been shown [205,206] that if quenched random hopping probabilities are 
assigned to each particle in the TASEP, there are small gaps between particles 
in the high-density congested phase but in the inhomogeneous low-density 
phase there is a macroscopically large empty region in front of the slowest 
particle (i.e., the particle with smallest hopping probability) behind which a 
platoon is formed. The phase transition from the low-density inhomogeneous 
phase (which consists of a macroscopic free region and a platoon) to the high- 
density congested phase is, in many respects, analogous to the Bose- Einstein 



i{t) ~ i("+i)/('»+2) 
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transition. 



In order to see this analogy, let us imagine that the empty sites are bosons 
and the state of a boson is determined by which particle it is immediately 
in front of. In the language of the ideal Bose gas, in the high-density phase 
the bosons are thinly spread over all the states. On the other hand, in the 
low-density phase there is a finite fraction of the empty sites are condensed 
in front of the slowest particle in such " Bose- Einstein- like condensed" state. 
The steady-state velocity of the particles is the analogue of the fugacity of the 
ideal Bose gas. What makes the system interesting is the fact that the platoon 
appears at low-density rather than at high density of the vehicles. 

The Bose-Einstein-like-condensation in the TASEP with quenched random 
hopping probabilities of the individual particles survives when the random 
sequential updating is replaced by parallel updating [206]. Finally, it is worth 
emphasizing that, the analogy with the ideal Bose gas is only formal as the 
empty sites in the TASEP are not non-interacting quantum particles. 

The qualitative features of the dynamical phases and phase transitions ob- 
served in the NaSch model with random braking probabilities, for v„iax = 1 
as well as for larger Vjnax-, are very similar to those described above for the 
TASEP with random hopping probabilities [204]. Typical snapshots of the 
system at three different stages of evolution from a random initial state are 
shown in Fig. 47. 




Fig. 47. Typical space-time diagram of the NaSch model with Vmax = 2 and c = 0.1 
and random braking probabilities. Each horizontal row of dots represents the in- 
stantaneous positions of the vehicles moving towards right while the successive rows 
of dots represent the positions of the same vehicles at the successive timesteps. 



102 



The typical size of the platoons ^(t) can be computed directly [204] by com- 
puting the correlation function (80) and identifying the separation r = Rq of 
the first zero-crossing of this correlation as ^{t). Following this procedure, it 
has been observed that ^{t) follows the power law (111) when the distribution 
of the random braking probabilities is given by P{p) — 2"{n + — p)"- 

10.2 Random v^ax 

The two important parameters of the NaSch model are p and Vmax- In the 

preceding subsection we have seen the effects of randomizing p assigning the 
same Vmax to all the vehicles. In this subsection, on the other hand, we inves- 
tigate the effects of randomizing Vmax, assigning a non-random constant p to 
every driver. 

The simplest possible model to investigate the effects of quenched randomness 
in Vmax is that considered by Ben-Naim et al. [209] which was discussed to 
motivate the phenomenon of platoon formation in the preceding subsection. 

In order to model traffic consisting of two different types of vehicles, say, for 
example, cars and trucks, of which a fraction ffast are intrinsically fast (say, 
cars) while the remaining fraction 1 — ffast ai'c intrinsically slow (say, trucks), 
Chowdhury et al. [195] assigned a higher Vmax (e-g-, Vmax = 5) to a fraction 
ffast of vehicles chosen randomly while the remaining fraction 1 — ffast were 
assigned a lower Vmax (e.g., v^ax = 3). As the density of the vehicles increases, 
the vehicles with higher Vmax find it more difficult to change lane in order 
pass a vehicle with lower Vmax ahead of it in the same lane. This leads to 
the formation of "coherent moving blocks" of vehicles each of which is led 
by a vehicle of lower Vmax [212]. Two main causes of traffic accidents, namely, 
differences in vehicles speeds and lane changes, are reduced considerably in this 
state thereby making this state of traffic much safer. It is worth mentioning 
that even a small number of slow vehicles in 2-lane models, where overtaking 
is possible, can have a drastic effect. For details we refer to [197] and the 
discussion in Sec. 9.2. 



10.3 Randomly placed bottlenecks on the roads and the maximum flux prin- 
ciple 

So far we have investigated the effects of two different types of quenched 
randomness both of which were associated with the vehicles (i.e., particles). 
We now consider the effects of yet another type of quenched randomness which 
is associated with the road (i.e., lattice). 
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In order to anticipate the effects of such randomness associated with the high- 
way, let us begin with the simplest possible caricature of traffic with a "point 
defect" [213]: a single "impurity" (or, "defect") site in the deterministic limit 
p = of the NaSch model with Vmax = 1- In this model, vehicles move forward, 
in parallel, by one lattice spacing if the corresponding site in front is empty; 
each vehicle takes Ti^p (> 1) timesteps to cross the "impurity" site but only 
one time step to cross a normal site when the next site is empty. The impurity 
sites acts like a blockage for all T^^p > 1. As explained in section 8.2.1, in the 
absence of the impurity, J = c for < c < 1/2 and J = 1 — c for 1/2 < c < 1. 
Note that, if the impurity is present, l/T^^p vehicle passes through the im- 
purity site per unit time. Therefore, the bottleneck created by the impurity 
introduces an upper cut-off of the flux, viz.. 1 /Timp- Obviously, J = c < 1 /Timp 
so long as c < Ci = l/Tj^p. Similarly, J = 1 — c < l/T^j^p for c > C2 where 
C2 = 1 — Ci- In the density interval Ci < c < 1 — Ci, the bottleneck at the 
impurity is the flow-limiting factor and, hence, in this regime, J = l/T^^p is 
independent of c. Thus, in the simple caricature of traffic under consideration 
one would expect the flux to vary with density following the relation 

c if < c < Ci, 

J = I 1/Timp if Ci < C < C2, (112) 

1 — C if C2 < c < 1, 

where 

Ci = 7T-^ and C2 = ^ i^j^m, — (113) 

l + (Ai),„p l+{M)imp ^ ' 

and Timp = 1 + (Ai)j^p such that {/S.t)im,p — for the normal sites but 
{/S.t)imp > for the impurity site. 

The fundamental diagram, obtained numerically through computer simulation 
of the NaSch model with a single defect and non-zero p (Fig. 48) is in quahta- 
tive agreement with those of the fundamental diagram (112). The qualitative 
features of the fundamental diagram in Fig. 48 are also similar to those of the 
TASEP with a single defect [214] where the hopping probability q is smaller 
than that at all the normal sites. Equation (112) also indicates that the larger 
is {At)imp the lower is the maximum flux 1/[1 -|- {At)imp\ and the wider is the 
interval Ci < c < C2 over which the flux remains constant. 

What makes the problem of a single "point defect" nontrivial is that, over the 
interval Ci < c < C2 of the density of the vehicles, where the J is maximum 
and independent of c, the localized blockage has global effects whereby the traf- 
fic exhibits macroscopic phase segregation into high-density and low-density 
regions. Evidence for such macroscopic phase segregation can be obtained di- 
rectly from the density profiles (see Fig. 49). Fig. 49 implies that so long 
as c < Ci the particles will not pile up but a local increase of density will 
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Fig. 48. Fundamental diagram of the NaSch model with Vmax = 1 and a blockage 
site located at the site 1. The hopping probability in the bulk is given hy q = 0.75 
and at the defect by q = 0.25. 
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Fig. 49. Density profiles for different values of c in the NaSch model with Vmax = 1 
and a blockage site located at the site 1. The hopping probability in the bulk is 
given hy q = 0.75 and at the blockage site hy q = 0.25. 

compensate for the reduced local velocity at the blockage so that the flux 
around the blockage is identical to that far from it. However, if the global 
density exceeds ci, the particles pile up during the transient period leading to 
the phase-segregated steady state. Because of the particle-hole symmetry the 
phase-segregation does not take place if the particle density exceeds C2. 



We now develop a semi-phenomenological theoryf^ for the NaSch model with 
Vmax = 1, non-zero p and a single "impurity" site assuming the steady-state 
to be phase-segregated, as demonstrated by computer simulation (Fig. 49). 
Naturally, this theory cannot explain the underlying mechanism that gives 
rise to the phase-segregated structure of the steady-state. But, as we shall see 
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A microscopic approach for deterministic dynamics can be found in [176]. 
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soon, it provides a good estimate of the flux in the phase-segregated regime. 
Our calculations are based on arguments similar to those suggested originally 
by Janowsky and Lebowitz [214] in the context of TASEP with a single defect. 

Using equation (78), the flux in the high-density and low-density regions, 
far from their interface, are given by J/j = y^l — 4gc/i(l — Cft)j and 

Jt — I (l — y'l — 4gQ(l — q)^ and that across the defect bond is given by 

Jde/ ~ I (^1 — ^1 — 4gf(iC?i(l — Q)j. Since, in the steady state, the flux is same 
across the entire system, we must have qch{l — Ch) — ?q(1 — q) and, hence, 

Ch = Ci or Ch = l- Q. (114) 

The condition Ch = Q is satisfied by the uniform density profiles whereas the 
condition = 1 — q is satisfied by the phase-segregated density profile (see 
Fig. 49). Moreover, using the condition Jh = Jdef = Ji we get 

Ch{l - Ch) = q(1 - q) ~ rch{l - Ci) (115) 

where r = qa/q < 1 may be interpreted as the "transmission probability" or 
"permeability" of the blockage. Prom the (approximate) equations (115) we 
get 

Ch ^ = and Q ~ = (116) 

r + 1 p + Pd r + 1 p + Pd 



and, hence, 

iqr 



1- Wl- 



;i + r)2 



(117) 



The estimate (117) is in good agreement with the numerical data (Pig. 48) 
obtained from computer simulation [190]. Moreover, the estimates q and Ch 
are also in good agreement with q and Ch, respectively, in Pig. 49. 

Note that q and Ch depend only on r and are independent of c. Moreover, the 
estimates (116) of q and are in excellent agreement with ci and c^, respec- 
tively, in Pig. 48. At first sight, these two results may appear surprising and 
counter-intuitive. But, we'll now show that these are related to the mechanism 
of the phase-segregation. Conservation of the vehicles demand that 

cL = Chh + ce^ Chh + ce{L - h) (118) 

where h and i — L — h are the lengths of the high-density and low-density 
regions, respectively. Thus, 

h ^ c-Q ^ c{l + r)-r 
L Ch-Ci 1-r 
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The equation (119) shows that /i/L^Oasc— >q and /i/L — > 1 as c — > c/j. 
Therefore, keeping r fixed as the density is increased beyond Ci = q, the 
densities of the two regions remain fixed but the high-density region grows 
thicker at the cost of the length of the low-density region as more and more 
vehicles pile up and, eventually, ai c — C2 — Ch the low-density region occupies 
a vanishingly small fraction of the total length of the system signaling the 
disappearance of the phase segregation. Interestingly, recasting the expressions 
for Ch and q as q = 1/[1 + {At)imp\ and Ch = {At)imp/[1 + {At)imp], where 
{At)imp = ^/^-i we find close formal analogies with ci and C2, respectively, in 
equation (113) [215]. 

Schiitz [216] considered a TASEP with sublattice-parallel update (see App. 
A) where the motion of the particles is deterministic (i.e., q — 1) everywhere 
except at a defect site where they move with the probability 5^ < 1 (i.e., 
r = < 1). Exact solution is possible through a mapping on a 6- vertex 
model. Later, a solution using the matrix-product Ansatz was presented in 
[217]. Except for minor differences, the quahtative features of the results do not 
differ from the corresponding approximate results obtained for ^ 1 [214]. 
Qualitatively similar phase segregation phenomena have also been observed 
in a related model [202]. 

The qualitative features of the fundamental diagram do not change signifi- 
cantly if the "point-like" defect (or, impurity) is replaced by an "extended" 
defect [190], i.e. a few consecutive defect sites. However, with increasing length 
of the defect, the maximum value of the flux decreases monotonically and ap- 
proaches the maximum flow of the homogeneous system where the hopping 
probability associated with each of the bonds is identical to that associated 
with the defects in our model (Fig. 50). From Figs. 50 and 51 we conclude 
that the monotonic decrease of the flow with increasing length of the extended 
defects, leads to a larger difference Ch — ct between the densities of the high- 
density and the low-density regions of the phase-segregated steady-state. 

Next, instead of a single point-hke or extended "defect", let us consider the 
more general case of quenched "disorder" in the NaSch model with v^ax — 1 
where the quenched random hopping probabilities Q'jj+i = 1— Pjj+i are chosen 
independently from some probability distribution P(g), for the hopping from 
the cell j to the cell j -l- 1 (j = 1,2,...,L). For a given realization of the 
disordered system, every vehicle hopping from a given cell i to the next cell 
i -|- 1 must hop with the same probability and a given vehicle hops across 
different bonds, in general, with different probabilities assigned to these bonds 
as it moves forward with time. A similar generalization of the TASEP has also 
been studied [137]. We shall refer to this model as disordered TASEP (i.e., 
DTASEP). 
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Fig. 50. Fundamental diagram for the NaSch model with Vmax = 1 and defects of 
different length L^. Again p = 0.75 and pd = 0.25 is chosen. For comparison the 
fundamental diagram of the "fast" (p=0.75) and "slow" homogeneous systems are 
shown. 
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Fig. 51. Density profiles for different lengths of the defect in the NaSch model with 
Vmax = 1- Again p = 0.75 and p = 0.25 are used and the average density is given 
by c = 0.50. The defect extended over the first sites of the system. 

Suppose, q are chosen from the binary distribution 




P{qjj+i = Qd) = /, P{qj,j+i = g) = 1 - / (120) 

i.e., a fraction / of the bonds have a permeability r < 1 while the remaining 
fraction 1 — / have unit permeability. A mean-field theory has been developed 
[137] (see Appendix G for details) for computing the fundamental diagram 
of the DTASEP. The flux in this model has interesting symmetry properties 
under the operations of "charge conjugation" (which interchanges particles 
and holes), "parity" (which interchanges forward and backward hopping rates 
on each bond and "time reversal" (which reverses the direction of the current) 
[137,218]. 
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The quenched disorder in these "disordered" models can be viewed as "point- 
hke impurities" distributed randomly over the lattice. But, the qualitative fea- 
tures of the fundamental diagram of DTASEP are similar to those observed 
for a single point-hke defect and those for a single extended defect. Although 
the random distribution of the point-like impurities leads to a "rough" density 
profile for all densities, in an intermediate regime of density, phase-segregated 
steady-sates with macroscopic high- and low-density regions have been iden- 
tified. 

What is the underlying mechanism for the "macroscopic" phase segregation 
in all the models DTASEP [137] ? Let us denote the stretches of bonds with 
permeability 1 by X and the stretches of bonds with permeability r by Y . The 
two parabolas in Fig. 52 are the two steady-state fundamental diagrams for 
the two pure reference systems consisting of all X and all Y, respectively. Since 
the flux must be spatially constant in the steady-state, the possible densities 
are given by the four intersections of the line J = Jq with the two parabolas. 
If the average density is less (greater) than 1/2 then the two possible densities 
are Ci and C2 (ca and C4). The variation of density between ci and C2 (or C3 
and C4) in the X and Y stretches is merely micro phase-segregation while, on 
a macroscopic scale, the density remains uniform. For simplicity, we assume 
that the density in each stretch of like bonds is uniform. The global density 
of the system is approximately c ~ (1 — /)ci,4( Jq) + /c2,3( Jo) where / is given 
by equation (120). However, as the density increases the flux also increases 
till it attains the maximum allowed flux of the pure system consisting of all 
Y (this happens at a global density smaller than 1/2). What happens when 
the density increases further? According to the "maximum current principle" 
[162], no further increase of the flux is possible and the excess density is taken 
care of by increasing the density in some of the X stretches from Ci to C4 
(or, vice versa if c > 1/2). This conversion takes place adjacent to the largest 
stretch of Y bonds where the density also changes from C2 to C3 (or, vice versa 
if c > 1/2) to accommodate the additional particles added to the system. 
This leads to the macroscopic phase segregation as the system consists of two 
macroscopic regions of two different mean densities- one with lower densities 
Ci, C2 in the X and Y stretches and the other with the higher densities C3, C4 
in the X and Y stretches. 

It is not difficult to generalize the DTASEP to disordered NaSch (DNaSch) 
model by replacing the random sequential updating by parallel updating. How- 
ever, we face subtle conceptual difficulties in extending DNaSch model further 
to arbitrary Vmax] if the position of a vehicle at time i -|- 1 is decided by q at 
its current position at time t, it may be forced to move v sites downstream by 
hopping over sites of even smaller g if v > 1. 

Some alternative parameterization of the defect (or, disordered) sites in the 
NaSch model for arbitrary v^ax have also been suggested. In [219,190] local- 
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Fig. 52. Origin of phase separation in DTASEP. Shown are the fundamental dia- 
grams of two reference non-disordered TASEPs. 

ized defects have been investigated where the randomization parameter is 
larger than in the rest of the system. Csahok and Vicsek [220] have consid- 
ered the blockages as sites with a "permeabihty" smaller than unity whereas 
the permeability of all the other sites is unity. This effectively reduces v^ax 
while the vehicle is at a blockage. On the other hand, Emmerich and Rank 
[221,154] considered a model of where the velocity of every vehicle in the re- 
gion occupied by the blockage (or, more appropriately, hindrance) at the time 
step t -|- 1 is half of that at time t, i.e., f „(t + 1) = w„(t)/2 is the n-th vehicle 
is located within the hindrance region. Some effects of static hindrances on 
vehicular traffic have also been investigated following alternative approaches, 
e.g., car- following theory [222]. From a practical point of view, ramps have 
effects very similar to those of a static defect. For the NaSch model this has 
been investigated in [223]. 



11 Other CA models of highway traffic 

All the CA models of highway traffic described so far are basically generahza- 
tions of the minimal CA model proposed originally by Nagel and Schrecken- 
berg [20]. We now describe a few other alternative minimal CA models and 
the interesting features of the corresponding results. 

11.1 Fukui-Ishibashi model 

The update rules of the Fukui and Ishibashi (FI) [120] model of single-lane 
highway traffic are as follows: 

If Vmax or more sites in front of the n-th vehicle is empty at the time step t, 
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then it has a probabihty 1—p to move forward by Vmax sites and a probabihty 
p to move forward by v^ax ^ 1 sites in the time step t + 1. However, if only d 
sites {d < Vmax) in front of the n-th vehicle are empty at time t then it moves 
by d sites in the next time step. Since no site can be occupied simultaneously 
by more than one vehicle, a vehicle must not move forward in the time step 
i + 1 if the site immediately in front of it is occupied by a vehicle at the time 
step t. The model becomes deterministic in both the limits p — and p — 1. 

The FI model differs from the NaSch model in two respects: (a) the increase 
of speed of the vehicles is not necessarily gradual and (b) the stochastic delay 
applies only to high-speed vehicles. The FI model, obviously, reduces to the 
NaSch model if Vmax = 1. A site-oriented mean-field theory [224] and a car- 
oriented mean-field theory [225,226] for the FI model have been developed for 
arbitrary Vmax and p. Note that the FI model is equivalent to the deterministic 
CA rule 184 (in the notation of Wolfram [18]) in the hmit v^ax — l,p — 0. 
Generalizations of the deterministic limit p = of the FI model have also been 
proposed [227,228]. 

According to the classification of Sec. 7.3, the FI model belongs to class I, i.e. 
the high-acceleration limit where no spontaneous jamming exists. In [156] an 
alternative high-acceleration variant has been proposed. Here only the accel- 
eration step {Ul) of the NaSch model is changed to f ^ Vmax, i-e. all vehicles 
accelerate immediately to the maximal possible velocity. The other update 
steps of the NaSch model are left unchanged. In contrast to the FI model, all 
vehicles are subject to the randomization step. The behaviour of this variant 
is therefore similar to that of the NaSch model, e.g. one finds spontaneous jam 
formation. 



11.2 Galilei-invariant vehicle-vehicle interaction and metastahility 

In the NaSch model it is postulated that the gap between a pair of successive 
vehicles is adjusted according to the velocity of the leading vehicle alone. In 
contrast, often in real traffic, drivers tend to adjust the gap in front taking 
into account the difference between the velocity of their own vehicle and that 
of the leading vehicle. The latter aspect of real traffic is captured by a recent 
model developed by Werth, Froese and Wolf (WFW) [229,230]. 

If both the following vehicle and the leading vehicle move with constant ac- 
celeration 6, then a collision between the two can be avoided provided 

gap + £{vlv) > vpvTr + K'"fv) (121) 
where v^y a-nd vpy a-re the velocities of the leading vehicle and the following 
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vehicle, respectively, is the reaction time of the following vehicle and 

is the distance covered by a vehicle with initial velocity v before it comes to 
a stop by moving with a constant deceleration b. Using equation (122) the 
condition (121) for avoiding collision can be written as 

V 

gap > VFvTr + -^{vfv - Vlv) (123) 

where v — {vfv + vlv)/'^ if the average velocity of the pair of vehicles under 
consideration. Therefore, a sufficient condition for avoiding a collision is 



gap > 



\vFVTr for VFV < VLV, ^^24) 

[vFvTr + ^^{vfv - vlv) for VfV > VlV- 

In the limiting case b ^ oo the sufficient condition (124) reduces to gap > 
VFvTr which is identical to the form of vehicle- vehicle interaction in the NaSch 
model if one chooses r,. as the unit of time. In the opposite limit r = 0, the 
sufficient condition (124) reduces to 

. fo for vfv < Vlv, 

gap > s , , (125) 

[^(vfv - Vlv) for vfv > vlv- 

Since this type of vehicle-vehicle interaction involves the difference of the ve- 
locities Vfv ~ Vlv, it is clearly invariant under a Galilean transformation 
and, hence, the name. The vehicle- vehicle interactions in real traffic may be 
somewhere in between the two limiting cases of NaSch model and the Galilei- 
invariant model. 

Suppose the indices n — 1 and n label the leading vehicle and the following 
vehicle, respectively, of a pair. The update rules suggested by WFW [230] for 
implementing the Galilei-invariant vehicle- vehicle interaction are as follows: 



Step 1: Acceleration. 

t;^^) = min('y„ + 1, v^ax) 
Step 2: Deceleration (due to other vehicles). 

v^^^ = min{v^^\dn - 1 + Vn-l) 

Step 3: Randomization. 

v^^^ ma,-x{v^^ — 1,0) with probabihty p 
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Step 4' Deceleration (due to other vehicles). 

^i^) = min(t;(3),4-l + ^i'2i) 



Step 5: Vehicle movement. 



Thus, the rule for deceleration (due to other vehicles) is applied twice. Step 
4 makes sure collisions are avoided. Since also the new velocity v^^\ of the 
preceding car enters, this step can not be performed in parallel for all cars. 
Instead it is performed sequentially, but the final configuration is independent 
of the starting point of this sequential updating. Step 4 has then to be applied 
twice in order to determine all velocities v^^ consistently. 

The rules as given above define the retarded version of the Galilei-invariant 
model. In the non-retarded version, in step 2 Vn-i is replaced by the new 
velocity v^li- To determine v^^ consistently for all cars, step 2 has then to be 
iterated Vmax ~ 1 times. 

The most interesting feature of the Galilei-invariant model is that its funda- 
mental diagram has a metastable branch although its update scheme involve 
neither cruise-control nor slow-to-start rules. The mechanism leading to the 
existence of metastable states is different from the models with slow-to-start 
rules (see Sec. 9.1.2). The outflow from jams is the same as in the NaSch model 
since it is independent of the interaction between vehicles. However, due to 
the inclusion of anticipation effects (i.e. the driver knows the velocity of the 
preceding vehicle) the free-flow is less sensitive to fluctuations. 



11.3 CA versions of the optimal-velocity model 



The traffic jams appear spontaneously in both the OV models and the GA 
models. However, in the OV models spontaneous formation of the jams arc 
caused by the non-linearity of the dynamical equations whereas in the CA 
models it is triggered primarily by the stochasticity oi the update "rules". The 
mechanism for the spontaneous formation of jams in real traffic may be a 
combination of these two. 

In the OV model, the control of velocity is given by the control of acceleration 
through the OV function which gives the optimal velocity for the current 
distance-headway. Thus, unhke the GA models like the NaSch model, the 
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vehicles in the OV models get an opportunity to avoid crash without any need 
to exert unphysically large deceleration. In fact, collision of vehicles may take 
place in naive discretization of the dynamical equations for the OV models 
unless special care is taken in the discretization process (see below). 

In the following we present several CA model analogues of the OV model that 
have been proposed by different authors. In principle, the NaSch model is also 
an OV model, but with a linear OV function, v{d) = mm[d — l,Vmax]- The 
first attempt to generahze this relation is due to Emmerich and Rank [231]. 
The update rules of their model are as follows: 

Step 1: Find largest gap. 

Find the vehicle with the largest gap to the next vehicle ahead. 



Step 2: Acceleration. 



m.in{Vmax,Vn + 1) 



Step 3: Deceleration due to other vehicles. 

Vn Md„_i,^„ iidn-l < Vmax, 

i.e. a vehicle with velocity j and i empty cells in of it (i.e. a gap d„ 
reduces its velocity to Mj j (0 <i,j< v^ax)- 



= i + l) 



Step 4-' Randomization. 



max;(t'„ — 1,0) with probability p 



Step 5: Vehicle movement. 



Step 6: Next vehicle. 

Repeat steps 2-5 for the next vehicle behind, i.e. proceed in the direction 
opposite to the motion of the vehicles. 



For the NaSch model with i>„ 



(NaSch) 



5 the matrix Mj ^ is given by 

/O o\ 
11111 
12 2 2 2 

12 3 3 3 
12 3 4 4 
yO 1 2 3 4 5/ 



(126) 



114 



A general matrix Mij has to satisfy certain conditions (e.g. Mjj < min(i,_7) 
and Mij < Mi^k tor j < k) to guarantee e.g. the absence of colhsions in the 
model. In order to model the fact that faster vehicles keep a relatively larger 
headway to the preceding vehicle, Emmerich and Rank suggested the following 
matrix: 





0^ 





11111 





1 2 2 2 2 





1 2 2 3 3 





1 2 3 3 4 




1 2 3 4 4y 



(127) 



Using a parallel update scheme, the model shows unrealistic behaviour in the 
free- flow regime, especially in the deterministic limit p <S 1. Here the funda- 
mental diagram is non-monotonic [232], as can be seen from a simple example 

for p = 0. At density c = 1 /7, the stationary state is of the form 5 5 5 

where numbers denote the velocity of vehicles and '.'an empty cell. At density 

c — 1/6, on the other hand, the stationary state is 4 4 4 Comparing 

the corresponding flows, one finds J(c = 1/7) = 5/7 > J(c = 1/6) = 2/3. At 
density c = 1/5, the stationary state is 4. ...4. ...4.... with flow J(c = 1/5) = 4/5 
which is again larger than the flow at c = 1/6 and corresponds to the maximal 
possible flow. This kind of behaviour persists even in the presence of random- 
ness (p > 0) [232]. In order to circumvent this problem, Emmerich and Rank 
had to introduce a special kind of ordered-sequential update, where first the 
vehicle with the largest gap ahead is updated. Then, the position of the next 
vehicle upstream is updated, and so on, using periodic boundary conditions. 

Emmerich and Rank also investigated more general rules where even for gaps 
larger than Vmax the velocity of the vehicles is reduced to a value v < Vmax- 

Later, a similar model has been proposed by Hclbing and Schrcckcnbcrg [233]. 
It is closer to the spirit of the original optimal- velocity model (see Sec. 6.2). 



Step 1: Vehicle movement. 



Step 2: Acceleration. 



V'^t + 1) = Vr,{t) + \[V^t{dn{t)) - Vr,{t)\ 
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Step 3: Randomization. 



Vn{t + 1) = V'^{t + 1) 



1 with probability p 01 v'^{t + 1) > 0) 



otherwise. 



Here \_y\ denotes the floor function, i.e. the largest integer i < y. In [233] 
various optimal velocity functions Vopt{d) have been used in order to fit exper- 
imental data. The simplest, but unrealistic, choice was Vopt{d) = mm{d,Vmax) 
where d is the distance-headway. The parameter A corresponds formally to 
the sensitivity parameter in the OV model where it determines the timescale 
of relaxation towards the stationary fundamental diagram. Such an interpre- 
tation is not possible for discrete time models. Here the main effect of the 
parameter A is a rescaling of the OV function. 

The naive discretization of the OV function produces some undesirable fea- 
tures of the model, e.g. the flow corresponding to the OV function is non- 
monotonic in the free-flow region. Furthermore one finds - for certain initial 
conditions - a breakdown of the flow at a finite density c* < 1. This indicates 
a lack of robustness of model against small modifications of the rules. 

For A < 1 the model is not intrinsically collision-free [234], in contrast to most 
other models discussed in this review. Problems occur e.g. when fast vehicles 
approach the end of a jam. For the Vopt{d) = inm{d,Vmax) and the 'realistic' 
choice A = 0.77 [233] collisions can occur for Vmax > 5 [234]. For A > 1, on 
the other hand, a backward motion of vehicles is possible. For a given OV 
function it is possible to derive conditions for the parameter A which ensure 
the absence of collisions and backward motion [234] . 

Nagatani [235] has suggested a CA model which combines the OV idea with the 
TASEP. It is the discrete analogue of the simphfied OV model (51) presented 
in Sec. 6.2. Here the n-th vehicle moves ahead with probability Vn{t) where 
Vn(t) is interpreted as velocity. This velocity is obtained by integration of the 
OV function 



where a > 0, Axn{t) = Xn+i{t) — Xn{t) and Axc is a safety distance. Further- 
more the velocity is restricted to the interval < Vn{t) < Vmax < 1- 




for Axn{t) > Axc 
for Axn{t) < Axc 



(128) 
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11.4 CA from ultra-discretization 



Several CA which can be interpreted as traffic models have been derived using 
the so-called ultra-discretization method (UDM) [236]. This approach allows 
to establish a direct connection between certain differential equations and 
CA. The problem in the derivation of CA from differential equations lies in 
the discretization of the 'state' (or dependent) variable. E.g. in a numerical 
treatment, only space and time variables are discretized. 

The basic procedure of the UDM is as follows: 1) Start from a nonlinear wave 
equation, e.g. the KdV equation or Burgers' equation; 2) Discretize space- 
and time-variables in a standard way to obtain the discrete analogue of the 
wave equation which is still continuous in the state variable Uj{t); 3) The 
discrete analogue is now ultra-discretized. Defining Uj{t) — e\Q.{uj{t)) (where 
e depends on the discretization Aa; and At of space and time) one can use 
the identity lim^^o ^ log(e"^/^ + e^^"^ + ■ ■ ■ ) = max[74, . . . ] to derive the CA 
analogue fo the nonlinear wave equation. 

By applying the UDM to Burgers' equation vt — 2vVx + Vxx one obtains the 
(L + l)-state, deterministic CA [237,238]. 

nj{t -|- 1) = nj{t) + min[M, nj-i{t), L — nj{t)] — min[M, nj{t), L — nj^iit)]. 

(129) 

nj{t) is the occupation number of cell j at time t. In contrast to most other 
CA discussed in this review multiple occupations of cells are allowed. The 
maximum number of particles which can occupy the same cell is given by L, 
i.e. < nj{t) < L. The model defined by (129) might therefore by interpreted 
as a simple model for a highway with L lanes where the effects of lane changes 
are completely neglected. In [238] it has also been suggested to interprete 
nj{t)/L for large L as a coarse-grained density. The parameter M denotes the 
maximum number of vehicles that can move from cell j to cell j + 1- 

For M > L — 1 the model reduces to the rule-184 CA, i.e. the NaSch model 
with Vjnax — 1 and p = 0. In the case L < 2M the fimdamental diagram 
looks similar to that of rule-184 CA, but it is 'degenerate' in the sense that 
qualitative different stationary states with the same flow exist. For L > 2M 
the fundamental diagram the fundamental diagram resembles that of rule- 
184 CA with a blockage site (see Sec. 10.3, especially Fig. 48). In the region 
M/L < c < {L — M)/L the flow takes the constant value M/L, i.e. the 
parameter M can be interpreted as a ffow limiter. 

In [238] also a generalization of (129) to higher velocities v^ax > 1 has been 
suggested. For v^ax — 2 the generalized update rule is given by 
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nj{t + l)^ nj{t) + %_2(t) - aj{t) 

+ mm[bj-i{t) - aj-i(t), L — nj{t) - aj„2(i)] 

— mm[bj{t) — aj{t), L — nj+i{t) — aj-i{t)]. (130) 

Here bj{t) = mm[nj{t), L — nj+i(t)] is the maximum number of vehicles at site 
j that can move and %(t) = mm[nj{t), L — nj+i{t), L — nj+2it)] is the number 
of vehicles that move two cells forward. The idea behind this dynamics is that 
first the vehicles try to move two cells forward. This is only possible if the two 
cells ahead are not fully occupied, i.e. L — nj+i{t) > and L — nj+2{t) > 0. 
Then m.m[bj(t) — djit), L — rij+iit) — aj-i(t)] vehicles move forward one cell. 

The model defined by (130) can be considered as a generalization of the Fukui- 
Ishibashi model (see Sec. 11.1) to which it reduces for L = 1. Note that the 
flow limiter M has been dropped in the extended model. 

The fundamental diagram of the model (130) has a structure similar to that 
shown in Fig. 34, i.e. states of high flow exist. Due to the higher velocity the 
degeneracy found in the model (129) is lifted. In the case L = 2, the c = 1/2 
states ... 1111 .. . and . . . 2020 . . . are degenerate in the simple model. In the 
Vmax — 2 case, however, in the state ...1111... the vehicles can increase 
their velocity, while in the state . . . 2020 . . . they can not due to the presence 
of the fully-occupied cells '2'. The high-flow states are unstable against local 
perturbations. 

Another CA model obtained by ultra-discretization of the modified KdV- 
equation has been suggested in [239]. It is a second-order CA since the con- 
figuration at time t + 1 does not only depend on the configuration at time t, 
but also on the previous one at time t — 1. The update rule for the position 
Xj (t) of vehicle j is given by 

Xj{t + 2) = Xj{t) Axj(i) max[0, /^Xj{t -M) - M] - max[0, /^Xj{t -M) - L] 

(131) 

where Axj(i) = Xjj^i{t) — Xj{t) is the gap and M and L are constants. 

Although this model is deterministic it exhibits start-stop waves similar to 
those found in the NaSch model [239]. It seems that a second-order deter- 
ministic CA might produce effects similar to those of noise in a first-order 
stochastic CA. One should note, however, that the rules allow vehicles to 
move backwards. Other models obtained using the UDM have been discussed 
in [240]. 



12 Cellular automata models of city traffic and road networks 
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12.1 Biham-Middleton-Levine model and its generalizations 

In the BML model [21], each of the sites of a square lattice represent the 
crossing of a east- west street and a north-south street. All the streets parallel 
to the x-direction of a Cartesian coordinate system are assumed to allow only 
single-lane cast-bound traffic while all those parallel to the y-direction allow 
only single- lane north-bound traffic. Let us represent the east-bound (north- 
bound) vehicles by an arrow pointing towards east (north) . In the initial state 
of the system, vehicles are randomly distributed among the streets. The states 
of east-bound vehicles are updated in parallel at every odd discrete time step 
whereas those of the north-bound vehicles are updated in parallel at every 
even discrete time step following a rule which is a simple extension of the 
TASEP: a vehicle moves forward by one lattice spacing if and only if the site 
in front is empty, otherwise the vehicle does not move at that time step. 

Thus, the BML model is also a driven lattice gas model where each of the 
sites can be in one of the three possible states: either empty or occupied by 
an arrow t or Note that the parallel update rules of the BML model is 
fully deterministic and, therefore, it may also be regarded as a deterministic 
CA. The randomness arises in this model only from the random initial con- 
ditions [241]. Suppose, A^_^ and N-^ are the numbers of the east-bound and 
north-bound vehicles, respectively, in the initial state of the system. If peri- 
odic boundary conditions are imposed in all directions, the number of vehicles 
in every street is conserved since no turning of the vehicles are allowed by the 
updating rules. In a finite L x L system the densities of the east-bound and 
north-bound vehicles are given by = N^/L"^ and Cf = N-^/L^, respectively, 
while the global density of the vehicles is c = -|- C|. 

Computer simulations of the BML model with periodic boundary conditions 
demonstrate that a first order phase transition takes place at a finite non- 
vanishing density c*, where the average velocity of the vehicles vanishes dis- 
continuously signaling complete jamming; this jamming arises from the mutual 
blocking ("grid-locking") of the flows of east-bound and north-bound traffic at 
various different crossings (see [242] for the corresponding results of the BML 
model with open boundary conditions) . 

At concentrations just above c*, in the jammed phase, all the vehicles to- 
gether form a single cluster which is stretched along the diagonal connecting 
the south-west to the north-east of the system. In other words, the lowest- 
density jammed configurations consist of a single diagonal band where the — > 
and I occupy nearest-neighbour sites on the band in a zigzag manner. With 
further increase of density more and more vehicles get attached to the band in 
the form of off-diagonal branches and the infinite cluster of the jammed vehi- 
cles looks more and more random. Thus, in general, a typical infinite cluster of 
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the jammed vehicles consists of a "backbone" and "danghng vehicles" which 
are the analogs of the "backbone" and the "dangling ends" of the infinite per- 
colation clusters in the usual site/bond percolation [243]. However, in contrast 
to the infinite percolation cluster in the usual random site/bond percolation, 
the infinite spanning cluster of vehicles in the BML model emerges from the 
self-organization of the system. Nevertheless, concepts borrowed from percola- 
tion theory have been used to characterize the structure of the infinite cluster 
of jammed vehicles in the BML model at c > c^, [244,245]. The distribution of 
the waiting times of the vehicles at the signals (i.e., at the lattice sites) has 
also been investigated through computer simulations [246,247]. 



12.1.1 Poor man's mean- field estimates for the BML model 

If one is not interested in detailed information on the " structure" of the dy- 
namical phases, one can get a mean-field estimate of c^, by carrying out a 
back-of-the-envelope calculation. Suppose, and denote the average den- 
sities while and denote the average speeds of the east-bound and north- 
bound vehicles, respectively. In order to estimate one has to take into ac- 
count interaction of the east-bound (north-bound) vehicles not only with the 
north-bound (east-bound) vehicles [248] but also, in a self-consistent manner, 
with other east-bound (north-bound) vehicles [249]. Following the arguments 
of Appendix I, one can show that in the symmetric case = C| = c, a 
self- consistency equation for the speed is 




(132) 



for c < c^.. The critical density c* is determined by the condition that at c > c=k 
the equation (132) does not give a real solution. Hence we get c = c* = 6 — 
-\/32 ~ 0.343 which, in spite of the approximations made, is surprisingly close 
to the corresponding numerical estimate obtained from computer simulation 
[21]. Moreover, the mean-field estimate c* ~ 0.343 is also consistent with the 
more rigorous resultp^ [251,253] that is strictly less than 1/2. The BML 
model in three dimension, although not relevant for vehicular traffic, has also 
been studied [254]. 



Ishibashi and Pukui [250] claimed that complete jamming can occur in the BML 
model only for c = 1. However, a plausible flaw in their arguments was pointed out 
by Chau et al. [251]. In [252] it has been argued that oc i.e. c* = in the 

thermodynamic limit. 
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12.1.2 Mean- field theory of the BML model 

Recall that the occupation variables n{i;t) in the NaSch model describe the 
state of occupation of the sites i {i — 1,2,...) by the vehicles on the one- 
dimensional highway. A scheme for a truly microscopic analysis of the BML 
model begins [255] by introducing the corresponding generalized occupation 
variables n^{x,y;t) and n^{x,y;t), which describe the state of occupation of 
the sites {x, y) by the north-bound and east-bound vehicles, respectively, on 
the two-dimensional street-network. The analogs of the equations (64) and 
(65) (see Appendix I for details) have analogous physical interpretations. As 
usual, in the naive SOMF approximation one neglects the correlations between 
the occupations of different sites [255]. Since none of the sites is allowed to be 
occupied by more than one particle at a time, Pauli operators may be used to 
develop an analytical theory of vehicular traffic [256] but one should keep in 
mind that the particles representing the vehicles are purely classical and the 
system does not have any quantum mechanical characteristics. 

12.1.3 Generalizations and extensions of the BML model 

The BML model has been generalized and extended to take into account 
several realistic features of traffic in cities. 

• Asymmetric distribution of the vehicles: 

Suppose the vehicles are distributed asymmetrically among the east-bound 
and north-bound streets [248], i.e., c_> ^ C|. For convenience, let us write Cf = 
cfa and c^ — c{l — fa) where /„ is the fraction of the vehicles moving towards 
north. Clearly fa = 1/2 corresponds to the symmetric case = C| = c/2. On 
the other hand, /a = {fa = 1) correspond to the extreme asymmetric case 
where all the vehicles are east-bound (north-bound). Obviously, the absence of 
grid-locking in the extreme limits fa — and fa — i rules out the possibility 
of BML-like complete jamming transition, i.e., — 1 for both fa — and 
/a = 1. Moreover, decreases with decreasing asymmetry in the distribution 
of the vehicles; is the smallest for fa = 1/2, i.e., symmetric distribution of 
the vehicles. These results can be presented graphically by plotting the curve 
c*(/a) in the phase diagram on the c — fa plane [248]. 

• Unequal maximum velocities: 

In the BML model both east-bound and north-bound vehicles can move by a 
maximum of one lattice spacing at a time and, therefore, the average speeds of 
both types of vehicles can never exceed unity. On the other hand, recall that 
in the Fukui-lshibashi model [120] of highway traffic vehicles can move up to 
a maximum of M lattice sites at a time and, hence, have average velocities 
larger than unity. Incorporating similar high-speed vehicles Fukui et al. [257] 
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generalized the BML model where the east-bound vehicles are allowed torn 
ove by M sites at a single time step while the north-bound vehicles can move 
by only one lattice spacing. 

• Overpasses or two-level crossings: 

The BML model has been extended to take into account the effects of over- 
passes (or two- level crossings) [258] . A fraction fo of the lattice sites in the 
BML model are randomly identified as overpasses each of which can accommo- 
date up to a maximum of two vehicles simultaneously. The overpasses weaken 
the grid-locking in the BML model. Therefore, is expected to increase with 
increasing fo- Besides, is expected to be unity if fo = 1. Naturally, we ad- 
dress the question: does jamming disappear (i.e., c* is unity) only at fo = 1 
or for even smaller values of /„? In order to answer this question we extend 
the self-consistent mean-field arguments, which led to the equation (132), in- 
corporating the effects of the overpasses thereby getting the generalized self- 
consistency equation [259] 



1 

— Vi — V = - 

in the symmetric case = C| = c. The equation (133) reduces to the equation 
(132) in the limit = 0. It predicts that, if /o 7^ 0, a moving phase exists in the 
BML model with overpasses for vehicle densities c < c* = (6 — -\/32)/(l — fo) 
and that the jammed phase disappears altogether (i.e., c* becomes unity) for 
fo > 1 — 0.343 = 0.657, an underestimate when compared with the corre- 
sponding computer simulations. These results can be presented graphically by 
drawing the curve (/o) in the mean-field phase diagram on the c — fo plane 
not only for the symmetric distribution (i.e., for / = 1/2) but also for asym- 
metric distributions of the vehicles among the east-bound and north-bound 
streets [259]. 

• Faulty traffic lights: 

The effects of faulty traffic lights have been modeled by generalizing the BML 
as follows [260]: a fraction fu of the sites are identified randomly as the lo- 
cations of the faulty traffic lights. Both the north-bound vehicles currently 
south of the faulty traffic lights and the east-bound vehicles currently west of 
a faulty traffic light are allowed to hop onto the empty crossings where the 
faulty traffic lights are located, irrespective of whether the corresponding time 
step of updating is odd or even. If an east-bound vehicle and a north-bound 
vehicle simultaneously attempt to enter the same crossing, where the faulty 
traffic light is located, then only one of then is allowed to enter that crossing 
by selecting randomly. 



1-/0 



c + 



1 + 



fo 



- 4(1 - fo)c 



(133) 
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Since a north-bound (east-bound) vehicle will be able to move forward al- 
though only cast-bound (north-bound) vehicles would have moved forward 
had fti been zero, the average velocity of the vehicles is expected to increase 
with increasing fn. However, with the increase of fu there is an increasing 
likelihood that an east-bound (north-bound) vehicle would be blocked by a 
north-bound (east-bound) vehicle at a faulty traffic light. Thus, the increas- 
ing density of faulty traffic lights increases the effect of grid-locking thereby 
decreasing c*. 

• Static hindrances: 

The BML model has been extended to incorporate the effects of static hin- 
drances or road blocks (e.g., vehicles crashed in traffic accident), i.e., stagnant 
points [261,262]. A vehicle, which stays at a normal site for only one time step 
before attempting to move out of it, stays at a point-like blockage for Tp time 
steps before attempting to move out of it. Obviously, the longer is Tp, the 
lower is the corresponding c*. The time-dependent phenomenon of spreading 
of the jam from the blockage site during the approach of the system to its 
jammed steady-state configurations has also been investigated [261,263]. 

• Stagnant street: 

Let us consider the effects of a stagnant street, where the local density Cg of the 
vehicles is initially higher than that in the other streets [264], on the traffic flow 
in the BML model. The stagnant street, effectively, acts like a "line-defect", 
rather than a "point-defect". However, in contrast to the static roadblocks, a 
stagnant street offers a time-dependent hindrance to the vehicles moving in 
the perpendicular streets. As intuitively expected, the jamming transition has 
been found to occur at a lower global density when the local density in the 
stagnant street is higher. 

• Independent turning of the vehicles: 

Let us now generalize the BML model by assigning a trend or preferred di- 
rection of motion, Wn{x, y), to each vehicle n (n = 1, 2, ...) located at the site 
X, y. According to this definition, the vehicle n, located at x, y jumps to the 
next site towards east with probability Wn{x, y) while 1 — Wn{x, y) is the cor- 
responding probability that it hops to the next site towards north [265]. In this 
generalized model vehicles can take a turn but the processes of turning from 
east-bound (north-bound) to north-bound (east-bound) streets is stochastic. 
For simphcity, suppose, A^/2 vehicles are assigned Wn{x,y) = 7 while the re- 
maining A^/2 vehicles are assigned Wn{x,y) = 1 — 7 where < 7 < 1/2; this 
implies that N/2 vehicles move preferentially east-ward whereas the remain- 
ing N/2 vehicles move preferentially towards north. In the limit 7 = no 
vehicle can turn and we recover the original BML model with deterministic 
update rules. The most dramatic effect of the stochastic turning is that the 
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discontinuous jump A{v) of the average velocity (v) decreases with increasing 
7 and, eventually, the first order jamming transition ends at a critical point 
where A{v) just vanishes. 

• Jam-avoiding turn and drive: 

In the model of turning considered in references [265] a vehicle turns stochas- 
tically independently of the other vehicles. In real traffic, however, a vehicle 
is hkely to turn if its forward movement is blocked by other vehicles ahead of 
it in the same street. Therefore, let us now consider a model [266] where an 
east-bound (north-bound) vehicle turns north (cast) with probability ptum if 
blocked by another vehicle in front of it. Computer simulations of this model 
shows that c^iptum) increases with increasing pj„j.n- In a slightly different model 
[267], on being blocked by a vehicle in front, an east-bound (north-bound) ve- 
hicle hops with probability pja to the next east-bound (north-bound) street 
towards north (east). 

• A single north-bound street cutting across east-bound streets 

Let us now consider a special situation where only one north-bound street 
exists which cuts all of the L equispaced mutually parallel east-bound streets 
of length L [268]. This situation can be modeled as a L x L square lattice and 
each cell on the north-bound street can be in one of the three allowed states, 
namely, either empty or occupied by a — > or a j. But, in contrast to the 
BML model, there are two allowed states for each cell (outside the crossings) 
on the cast-bound street; these can be either empty or occupied by only — >. 
Since no grid-locking is possible with only one north-bound street, complete 
jamming occurs trivially in this case only if each of the cells cither on the cast- 
bound streets or on the north-bound street or on all the streets are occupied 
simultaneously by vehicles. Nevertheless, at any finite nonvanishing density, 
the crossings of the east-bound and north-bound streets act like hindrances 
with finite non-vanishing permeability for the fiow of the east-bound traffic. 
Obviously, the higher is the density of the north-bound vehicles, the lower is 
the permeability and the stronger is the rate-limiting factor of the bottleneck. 
A comparison of this problem with that of highway traffic in the presence of 
static hindrances [214,137] explains not only why the flux along the east-bound 
streets exhibits a flat plateau over an intermediate density regime but also why 
the plateau appears at lower values of the flux with increasing density of the 
north-bound vehicles [268]. A mean-fleld theory for the macroscopic phase 
segregation in this model has been developed by appropriately modifying that 
for the similar phenomenon in the one-dimensional models of highway traffic 
in the presence of static hindrances. It is worth emphasizing here that in [268] 
each site on the east-bound streets has been interpreted as a cell, which can 
accommodate one vehicle at a time, rather than as a crossing of the east- 
bound street with a north-bound street. Conceptually, this is an extension of 
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the BML model. Finally, we briefly mention that in [269] the phase diagram 
of a system consisting of one east-bound and one north-bound street with one 
crossing has been investigated. 

• Green-wave synchronization 
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Fig. 53. Comparison of the update procedure in the BML model (left) and the GW 
model (right) (from [270]). 

Often the traffic lights along the main streets in cities are synchronized to allow 
continuous flow; this is usually referred to as "green- wave" synchronization. 
A green-wave (GW) model has been developed [270] by replacing the parallel 
updating scheme of the BML model by an updating scheme which is partly 
backward-ordered sequential (see Appendix A for a general explanation of 
this update scheme). At odd time steps, an east-bound vehicle moves by one 
lattice spacing if the target site was empty at the end of the previous time 
step or has become empty in the current time step (this is possible because 
of the backward-ordered sequential updating at every time-step). Similarly 
the positions of the north-bound vehicles are updated at every even time step. 
The main difference between the BML model and the GW model (sec Fig. 53), 
arising from the different updating schemes, is that in the GW model vehicles 
move as "convoys" (a cluster of vehicles with no empty cell between them) 
thereby mimicking the effects of green- wave synchronization of the traffic lights 
in real traffic. The jamming transition in the GW model has been investigated 
by a combination of a mean-fleld argument and numerical input from computer 
simulations [270]. 

• More realistic description of streets and junctions 
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At first sight the BML model may appear very unrcaUstic because the vehicles 
seem to hop from one crossing to the next. However, it may not appear so 
unrealistic if each unit of discrete time interval in the BML model is interpreted 
as the time for which the traffic lights remain green (or red) before switching 
red (or green) simultaneously in a synchronized manner, and over that time 
scale each vehicle, which faces a green signal, gets an opportunity to move 
from j-th crossing to the {j + l)-th (or, more generally [257], to the (j -|-r)-th 
where r > 1). 

In the original version of the BML model the vehicles are located on the lat- 
tice sites which are identified as the crossings. Brunnet and Goncalves [271] 
considered a modified version where, instead, the vehicles are located on the 
bonds and, therefore, never block the flow of vehicles in the transverse di- 
rection. Consequently, in this version of the CA model of city traffic jams of 
only finite sizes can form and these jams have finite lifetime after which they 
disappear while new jams may appear elsewhere in the system; an infinitely 
long-lived jam spanning the entire system is possible only in the trivial limit 
c = 1. In contrast, Horiguchi and Sakakibara [272] generalized the BML by 
replacing each of the bonds connecting the nearest-neighbour lattice sites by a 
bond decorated with an extra lattice site in between. In [273] a generalization 
to s extra lattice sites between crossings has been presented. However, the ve- 
hicles are still allowed to hop forward by only one lattice spacing in the model 
of Horiguchi and Sakakibara. Moreover, generalizing the rules for turning of 
the vehicles in Ref. [265], Horiguchi and Sakakibara also allowed probabilistic 
turning of the vehicles in their model. The model exhibits a transition from 
the flowing phase to a completely jammed phase. 

The streets in the original BML model were assumed to allow only one-way 
traffic. This restriction was relaxed in a more realistic model proposed by 
Preund and Poschel [274] which allows both- way traffic on all the streets. 
Thus, each east-west (north-south) street is implicitly assumed to consist of 
two lanes one of which allows east-bound (north-bound) traffic while the other 
allows west-bound (south-bound) traffic. Moreover, each site is assumed to 
represent a crossing of a east-west street and a north-south street where four 
numbers associated with the site denote the number of vehicles coming from 
the four nearest- neighbour crossings (i.e., from north, south, east and west) 
and queued up at the crossing under consideration. So, in this extended version 
of the BML model, each site can accommodate at most AQ particles if each of 
the four queues of vehicles associated with it is allowed to grow to a maximum 
length Q. 

In the model proposed by Freund and Poschel [274], initially, each of the 
vehicles is assigned a site selected at random, the queue to which it belongs 
(i.e., whether it is approaching the crossing from north,south, east or west) 
and the desired direction (i.e., left, right or straight) for its intended motion 
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at the next time step. At each discrete time step a vehicle is allowed to move 
forward in its desired direction of motion by one lattice spacing provided (a) 
it is at the front of the queue in its present location and (b) there are fewer 
than Q vehicles queued up at the next crossing in the same desired direction 
of motion. Once a vehicle moves to the next crossing it finds a place at the 
end of the corresponding queue at the new crossing while the vehicles in the 
queue it left behind are moved "closer to the crossing" by one position (by 
mere relabeling as no physical movement of the vehicles in the queue takes 
place exphcitly) . Various reasonable choices for the rule which determines the 
desired direction of each vehicle at every time step have also been considered. 

The finite space of the streets in between successive crossings do not appear 
explicitly in the extension of the BML model suggested by Freund and Poschel 
[274] although it is more reahstic than the BML model because it implicitly 
takes into account the possibility of formation of queues by vehicles approach- 
ing one crossing from another. Chopard et al. [275] have developed a more 
realistic CA model of city traffic where the stretches of the streets in between 
successive crossings appear explicitly. In this model also each of the streets 
consist of two lanes which allow oppositely directed traffic. The rule for imple- 
menting the motion of the vehicles at the crossing is formulated assuming a 
rotary to be located at each crossing. Depending on the details of the rules to 
be followed by the vehicles at the rotary, the system can exhibit a variety of 
phenomena. For example, the flow can be metastable at all densities if each of 
the vehicles on the rotary is required to stop till the destination cell becomes 
available for occupation [275]. Moreover, the bottleneck created by the vehi- 
cles on the rotaries at the junctions can lead to a plateau in the fundamental 
diagram which is analogous to that caused by a static hindrances on a highway 
[214,137]. 

12.2 Marriage of the NaSch model and the BML model; a "unified" CA model 
of city traffic 

If one wants to develop a more detailed "fine-grained" description of city traffic 
than that provided by the BML model then one must first decorate each bond 
[268,272] with D — 1{D>1) sites to represent D — 1 cells in between each pair 
of successive crossings [268,275] thereby modeling each segment of the streets 
in between successive crossings in the same manner in which the entire highway 
is modelled in the NaSch model. Then, one can follow the prescriptions of the 
NaSch model for describing the positions, speeds and accelerations of the 
vehicles [275,276] as well as for taking into account the interactions among 
the vehicles moving along the same street. Moreover, one should fiip the color 
of the signal periodically at regular interval of T (T 3> 1) time steps where, 
during each unit of the discrete time interval every vehicle facing green signal 
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should get an opportunity to move forward from one cell to the next. Such a 
CA model of traffic in cities has, indeed, been proposed very recently [277,278] 
where the rules of updating have been formulated in such a way that, (a) a 
vehicle approaching a crossing can keep moving, even when the signal is red, 
until it reaches a site immediately in front of which there is either a halting 
vehicle or a crossing; and (b) no grid-locking would occur in the absence of 
random braking. 

Let us model the network of the streets as a, N x N square lattice. For simplic- 
ity, let us assume that the streets parallel to x and y axes allow only single- lane 
east-bound and north-bound traffic, respectively, as in the original formula- 
tion of the BML model. Next, we install a signal at every site of this N x N 
square lattice where each of the sites represents a crossing of two mutually 
perpendicular streets. We assume that the separation between any two suc- 
cessive crossings on every street consists of D cells so that the total number 
of cells on every street is L = N x D. Each of these cells can be either empty 
or occupied by at most one single vehicle at a time. Because of these cells, the 
network of the streets can be viewed as a decorated lattice. However, unlike 
the BML model [21], which corresponds to D — 1, and the model of Horiguchi 
and Sakakibara [272], which corresponds to D — 2, D {< L) in this model is 
to be treated as a parameter. Note that D introduces a new length scale into 
the problem. 

The signals are synchronized in such a way that all the signals remain green for 
the east-bound vehicles (and simultaneously, red for the north-bound vehicles) 
for a time interval T and then, simultaneously, all the signals turn red for the 
east-bound vehicles (and, green for the north-bound vehicles) for the next 
T timestcps. Clearly, the parameter T introduces a new time scale into the 
problem. Thus, in contrast to the BML model, the forward movement of the 
individual vehicles over smaller distances during shorter time intervals are 
described exphcitly in this "unified" model. 

If no turning of the vehicles is allowed, as in the original BML model, the 
total number of vehicles on each street is determined by the initial condition, 
and does not change with time because of the periodic boundary conditions. 
Following the prescription of the NaSch model, we allow the speed v of each 
vehicle to take one of the v„iax + 1 integer values w = 0, 1, ...^Vmax- Suppose, 
Vn is the speed of the n-th vehicle at time t while moving either towards east 
or towards north. At each discrete time step t — > t + 1, the arrangement of N 
vehicles is updated in parallel according to the following "rules": 

Step 1: Acceleration. 
Vn mm{vn + l,v 

max)- 

Step 2: Deceleration (due to other vehicles or signal). 
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Fig. 54. A segment of a east-west street, where the n-th and the (n + l)-th vehicles 
are located, is shown schematically together with one of the nearest signals. 

Suppose, dn is the gap in between the n-th vehicle and the vehicle in front of 
it, and s„ is the distance between the same n-th vehicle and the closest signal 
in front of it (see Fig. 54). 

Case J: The signal is red for the n-th vehicle under consideration: 
Vn niin(t;„, dn - 1, s„ - 1). 

Case II: The signal is green for the n-th vehicle under consideration: 
Suppose, r is the number of the remaining time steps before the signal turns 
red. Now there are two possibilities in this case: 

{i) When dn < Sn, then Vn — > min(f„, — 1). The motivation for this choice 
comes from the fact that, when dn < the hindrance effect comes from 
the leading vehicle. 

(ii) When dn > Sn, then w„ — > min(f„, dn — 1) if min(f„, (i„ — 1) x r > s„; 
else Vn — > min(f„, — 1). The motivation for this choice comes from the 
fact that, when dn > Sn, the speed of the n-th vehicle at the next time step 
depends on whether or not the vehicle can cross the crossing in front before 
the signal for it turns red. 

Step 3: Randomization. 

Vn — > max(f„ — 1, 0) with probability p (0 < p < 1); p, the random deceleration 
probability, is identical for all the vehicles and does not change during the 
updating. 

Step 4-' Vehicle movement. 

For the east-bound vehicles, Xn —>■ Xn + f n while for the north-bound vehicles, 

Un^ yn + Vn- 

The rule in case II of step 2 can be simplified without changing the overall 
behaviour of the model [278]: 

Case IF: 

If the signal turns to red in the next timestep (r = 1): 

Vn min(f„, Sn - 1, rf„ - 1) 

else 

Vn mm{vn,dn - 1). 
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These rules are not merely a combination of the rules proposed by Biham et 
al. [21] and those introduced by Nagel and Schreckenberg [20] but also involve 
some modifications. For example, unlike all the earlier BML-type models, a 
vehicle approaching a crossing can keep moving, even when the signal is red, 
until it reaches a site immediately in front of which there is either a halting 
vehicle or a crossing. Moreover, if p = every east-bound (north-bound) 
vehicle can adjust speed in the deceleration stage so as not to block the north- 
bound (east-bound) traffic when the signal is red for the east-bound (north- 
bound) vehicles. 

Initially, we put A^_^ and vehicles at random positions on the east-bound 
and north-bound streets, respectively. We update the positions and velocities 
of the vehicles in parallel following the rules mentioned above. After the initial 
transients die down, at every time step t, we compute the average velocities 
{vx{t)) and {vy{t)) of the east-bound and north-bound vehicles, respectively. 
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Fig. 55. A typical jammed configuration of the vehicles (N = 5, D = 8). The 
east-bound and north-bound vehicles are represented by the symbols — > and |, 
respectively. 

A phase transition from the "free- flowing" dynamical phase to the completely 
"jammed" phase takes place in this model at a vehicle density c^{D). The 
dependence on the dynamical parameters p, v^ax and T is not clear yet [278]. 
The data obtained so far from computer simulations do not conclusively rule 
out the possibility that the transition density only depends on the structure 
of the underlying lattice, similar to the percolation transition [243], and is 
independent of the dynamical parameters. The intrinsic stochasticity of the 
dynamics, which triggers the onset of jamming, is similar to that in the NaSch 
model, while the phenomenon of complete jamming through self-organization 
as well as the final jammed configurations (see Fig. 55) are similar to those in 
the BML model. The variations of {v^) and {vy) with time as well as with c, 
D, T and p in the flowing phase are certainly more reahstic that in the BML 
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model [277]. 
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Fig. 56. Time-depcndcncc of average speeds of vehicles. The symbols +, x, * and 
□ correspond, respectively, to {vx), {vy), and the fractions of vehicles with in- 
stantaneous speed V = 0, fxo and fyo, respectively. The common parameters are 
Vmax = 5,p = 0.1, D = 100, T = 100 and c = 0.1. The continuous line has been 
obtained from heuristic arguments given in [277]. 

The "unified" model has been formulated intentionally to keep it as simple 
as possible and at the same time capture some of interesting features of the 
NaSch model as well as the BML model. We believe that this model can be 
generalized (i) to allow traffic flow both ways on each street which may consist 
of more than one lane, (ii) to make more realistic rules for the right-of-the-way 
at the crossings and turning of the vehicles, (iii) to implement different types 
of synchronization or staggering of traffic hghts [276] (including green- wave) , 
etc. 



12.3 Practical applications of the models of vehicular traffic; On-line simu- 
lation of traffic networks 

A large fraction of the available resources are spent by the governments, par- 
ticularly in the industrialized developed countries, to construct more highways 
and other infrastructural facilities related to transportation. The car-following 
models, the coupled-map lattice models as well as the CA models have been 
used for computer simulation with a hope to utilize the results for on-line 
traffic control. For planning and design of the transportation network [19], for 
example, in a metropolitan area [279-281], one needs much more than just 
micro-simulation of how vehicles move on an idealized linear or square lat- 
tice under a specified set of vehicle- vehicle and road-vehicle interactions. For 
such a simulation, to begin with, one needs to specify the roads (including 
the number of lanes, ramps, bottlenecks, etc.) and their intersections. Then, 
times and places of the activities, e.g., working, shopping, etc., of individual 
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drivers arc planned. Micro-simulations are carried out for all possible different 
routes to execute these plans; the results give informations on the efficiency 
of the different routes and these informations are utilized in the designing 
of the transportation network [19]. Some socio-economic questions as well as 
questions on the environmental impacts of the planned transportation infras- 
tructure also need to be addressed during such planning and design. For a 
thorough discussion of these aspects we refer to the recent review by Nagel et 
al. [19]. 



13 Some related systems, models and phenomena 

In this section we will briefly present some stochastic models and phenomena 
which are somehow related to the main topic of this review. For some models 
the relation to traffic flow is rather obvious, e.g. there are toy models which 
share some features with traffic models, but which can be solved exactly. 
We also point out similarities in the description of other phenomena, e.g. 
granular flow or surface growth. Similarities also exist to systems from solid 
state physics, namely ionic conductors. This point will not discussed here, 
instead we refer to refs. [282,283]. Finally we hke to mention that there are 
some resemblances between the CA models of vehicular traffic and the CA 
models of driven diffusive Frenkel-Kontorova-type systems [284] . 

13.1 Generalizations of the TASEP 

The TASEP is probably one of the best studied models in nonequilibrium 
physics. Using powerful methods hke MPA or Bethe Ansatz recently it be- 
came possible to calculate not only simple expectation values in the station- 
ary states, but also more complicated quantities like diffusion constants or the 
large-deviation function. For a more detailed discussion and a list of references 
we refer to the recent reviews on the TASEP [22,174] (see also [15]). Several 
variants of the TASEP have been proposed. Most of those preserve the exact 
solvability. In the following we will discuss briefly some of the variants and 
generahzations of the TASEP discussed so far in the literature. 

• Partially asymmetric exclusion process 

An obvious generahzation of the TASEP is to allow hopping processes in 
both directions [285]. Here only results for the random-sequential update ex- 
ist, since for parallel dynamics ambiguities in the updating appear when two 
particles attempt to hop onto the same site. One flnds that the phase dia- 
gram looks essentially like that of the TASEP with three different phases (see 
Sec. 8.5). Recently it has been shown [286], however, that the high- and low- 
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density phases can be divided into three subphases (AI, All, AIII and BI, 
BII, Bill) instead of two in the TASEP case. Again the phase boundaries of 
these subphases are determined by the behaviour of the density profiles and 
the corresponding correlation lengths. 

• Multispecies models I 

Karimipour and collaborators [287] have developed a multispecies generaliza- 
tion of the TASEP which retains the solvability by MPA. It is similar to the 
disordered model discussed in Sec. 10.1 where each particle is characterized 
by a hopping rate (also called "velocity" in this context) Vj {j = 1, . . . ,Ns), 
i.e. there are Ng different "species'" of particles. In contrast to the models 
with quenched disorder discussed in Sec. 10, however, overtaking of particles 
is possible, i.e. the ordering of the particles is not fixed. 
More specifically the dynamics of the model is defined as follows: A particle 
is chosen at random (random-sequential update). If the particle is of species 
j and the cell to its right is empty, it hops there with rate vj. If the cell is 
occupied by a particle of species I and vi < vj, then they interchange there 
positions with rate Vj — Vi. This means that a fast particle can overtake a 
slower one with a rate proportional to their velocity difference. 

• Multispecies models II 

In several papers multispecies generalizations of the ASEP have been sug- 
gested which exhibit phase separation and spontaneous symmetry breaking. 

Arndt et al. [288] considered a system of positive and negative charged par- 
ticles diffusing on a ring in opposite directions. Positive particles move to an 
empty right neighbour and negative particles move to an empty left neighbour 
with the same rate A. Furthermore positive and negative particles on neigh- 
bouring sites can exchange their positions. The process — I > H — occurs with 

rate 1, and the inverse process H > — h with a different rate q. For equal 

densities of positive and negative particles the system exhibits three phases. 
For g < 1 in the thermodynamic limit, the system organizes itself into configu- 
rations consisting of blocks of the type 000 ■■■ + + + ■■ ■■ The dynamics 

out of theses states is extremely slow. Translational invariance is broken and 
the current vanishes. This phase is called 'pure phase'. For 1 < q < qc the sys- 
tem is in the 'mixed phase' which consists of two coexisting phase, the dense 
phase and the fluid phase. The dense phase where no vacancies exist covers 
a macroscopic region which shrinks to for q q^. It is remarkable that the 
current for 1 < q < q^ takes the value J = (g — l)/4 independent of the total 
density and the rate A. For q > qc the fluid phase extends through the whole 
system. There is no charge separation and density profiles are uniform. 

Similar results have been found by Evans et al. [289] in a shghtly different 
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model. The dynamics of their 3-species modelp^ is given by 

AB ^ BA, BC ^ CB, CA ^ AC, (134) 
1 1 1 

i.e. the rates are cychc in A, B and C and the numbers Na, Nb and Nc of 
particles of each species is conserved. 

In the case = Nb = Nc the dynamics satisfies the detailed balance 
condition with respect to a Hamiltonian with long-range asymmetric inter- 
actions. Stationary states are of the form A - ■ ■ AB ■ ■ ■ BC ■ ■ - C and exhibit 
phase separation, i.e. for large separations r the two-point function satisfies 

\imL-.oc[{AiAr) - {Ai){Ar)] > 0. 



13.2 Surface growth, KPZ equation and Bethe Ansatz 



In Sec. 8.1.2 we have explained how the NaSch model with Vmax = 1 can be 
mapped onto a stochastic surface growth model. This connection can be em- 
ployed to calculate several properties of the noisy Burgers and KPZ equation 
exactly. 

Gwa and Spohn [290] used the Bethe Ansatz (see e.g. [291]) to determine the 
spectrum of the stochastic Hamiltionian (see App. F) 

^ = -7 E h ■ '^.+1 - 1 + (^l^rfi - ^I^J+i)] (135) 

corresponding to the ASEP with random-sequential updating and periodic 
boundary conditions, aj = {aj,aj,aj) are the standard Pauli matrices at 
site j and e is the asymmetry of the hopping rates, qright = ^{^ + ^) 
qieft = 1(1 -e)- 

The "ground state" of Ti. has eigenvalue and is L-fold degenerate. For a 
fixed number N of up-spinsp^ every configuration has equal weight in the 
ground state. In order to determine the dynamical scaling exponent of the 
noisy Burgers and KPZ equations, Gwa and Spohn investigated the finite-size 
behaviour of the energy gap of Ti. Since 7i is non-hermitian its spectrum is 
complex. The first excited state is then defined as the eigenvalue with the 
smallest (positive) real part Egap. In [290] it was shown that Egap oc L~^/^ 
for e = 1 and = L/2. This implies that the dynamical exponent z for 
the stationary correlations of the KPZ equation is given by z = 3/2. The 



One species may be interpreted as a vacant site, 
i.e. a fixed number N of particles 
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dynamical exponent relates temporal and spatial scaling behaviour on large 
scale. Generalizations and related results can be found in [292]. 



By an extension of the Bethe Ansatz method of Gwa and Spohn, Derrida 
and Lebowitz [293] calculated the large deviation function (LDF) of the time- 
averaged current of the TASEP. The LDF is related to the total displacement 
Y{t), i.e. the total number of hops to the right minus the total number of hops 
to the left between time and time t. In the corresponding growth model Y{t) 
is the total number of particles deposited until time t. The LDF is then defined 
as 



f{y) = lim - In 



Prob 



'Y{t) 



t 



(136) 



where v = \imt^oo{Y{t)/t = ec(l — c)L'^/ {L — 1) is the mean current for a ring 
of finite size L and density c = N/L. The results of [293] have been extended 
and generalized in [294]. 



Apart from the treatment of finite systems, the BA can also be used to solve 
the master equation for an infinite system with a finite number of particles 
[295]. This allows e.g. to study the collective diffusion of two single particles. 



13.3 Protein synthesis 



You must have noticed in the earher sections that some of the models of traffic 
are non-trivial generalizations or extensions of the TASEP, the simplest of 
the driven-dissipative systems which are of current interest in nonequilibrium 
statistical mechanics [14,15,17]. Some similarities between these systems and 
a dynamical model of protein synthesis have been pointed out [296,15]. 

In a simplified picture of the mechanism of biopolymerisation ribosomes read 
the genetic information encoded in triplets of base pairs. They attach to 
one end of a messenger-RNA and then move along the chain after adding 
a monomer to a biopolymer attached to the ribosome. The type of monomer 
added depends on the genetic information read by the ribosome. When the 
ribosome reaches the other end of the m-RNA the biopolymer is fully synthe- 
sized and the ribosome is released. 

MacDonald et al. [297] have described the kinetics of this process using an 
ASEP-type model. The m-RNA is represented by a chain of L sites where 

each site corresponds to one triplet of base pairs. The ribosome is given by 
a hard-core particle covering r neighbouring sites (r ^ 20 — 30) which moves 
forward by one site with rate q. At the beginning of the chain particles are 
added with rate aq and at the end they are released with rate (3q. In the 
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idealized case r = 1 this is exactly the TASEP of Sections 8.1.1 and 8.5. The 
relevant case for the experiments is a = /5 < 1/2. The exact sohition of the 
TASEP allowed for an explanation of many aspects of the experiments [296] . 

13.4 Granular flow 

Another quasi-one-dimensional driven-dissipative system, which is also receiv- 
ing wide attention of physicists in recent years, is the granular material flowing 
through a pipe [3,4]. Since the fascinating phenomena (e.g. size segregation, 
convection, standing waves, localized excitations) found in granular materials 
have been subject of several excellent reviews [298] we discuss only briefly 
the similarities between the clustering of vehicles on a highway and particle- 
particle (and particle-cluster) aggregation process [299,300,209]. 

Obviously both highway traffic and granular flow through a pipe may be de- 
scribed as quasi-onedimensional systems consisting of discrete elements (ve- 
hicles, grains). The dynamics of these elements is determined by an intricate 
interplay between a driving force (driver, gravitation) and dissipation (brak- 
ing, inelastic scattering processes). These similarities already show that both 
systems can be described by similar approaches. One important difference be- 
tween traffic flow and granular flow exist, however. In granular flow density 
waves can move both in and against the direction of the flow whereas in traffic 
flow they only move backwards. 

An important success of a description of granular flow using the optimal veloc- 
ity model (Sec. 6.2) is the explanation of the experimentally observed f~^^^- 
behaviour of the power spectrum [301]. 

13.5 The bus route model 

The bus route model (BRM) [302] has been formulated as a one-dimensional 
lattice with periodic boundary conditions. The sites represent cells, each of 
which may be thought of as a bus stop and are labeled by an index i (i — 
1, 2, L) [302]. Two binary variables CTj and are assigned to each cell i: (i) 
If the cell i is occupied by a bus then cij = 1; otherwise (Xj = 0. (ii) If cell i 
has passengers waiting for a bus then = 1; otherwise Tj = 0. Since a cell 
cannot have simultaneously a bus and waiting passengers, let us impose the 
condition that a cell cannot have both ai — 1 and Ti — 1 simultaneously. Each 
bus is assumed to hop from one stop to the next. 

Next, let us specify the update rules (see Fig. 57): a cell i is picked up at 
random. Then, (i) ii — and — (i.e, cell i contains neither a bus nor 
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Fig. 57. The bus-route model (from [302]). 



waiting passengers), then r ^ 1 with probabihty A, where A is the probabiUty 
of arrival of passenger (s) at the bus stop, (ii) If cTj = 1 (i.e., there is a bus 
at the cell i) and CTj+i = 0, then the hopping rate /j, of the bus is defined as 
follows: (a) if Tj+i = 0, then /j, — at, but (b) if Tj+i = 1, then n — Pb, where 
ab is the hopping rate of a bus onto a stop which has no waiting passengers 
and (3b is the hopping rate onto a stop with waiting passenger (s). Generally, 
f3b < Oib, which reflects the fact that a bus has to slow down when it has to 
pick up passengers. We can set ab = I without loss of generality. When a bus 
hops onto a stop i with waiting passengers Tj is reset to zero as the bus takes 
all the passengers. Note that the density of buses c = N/L in a conserved 
quantity whereas that of the passengers is not. 

An ideal situation in this bus-route model would be one where the buses are 
evenly distributed over the route so that each bus picks up roughly the same 
number of passengers. However, because of some fluctuation, a bus may be 
delayed and, consequently, the gap between it and its predecessor will be larger 
than the average gap. Therefore, this bus has to pick up more passengers than 
what a bus would do on the average, because during the period of delay more 
passengers would be waiting for it and, as a result, it would get further delayed. 
On the other hand, the following bus has to pick up fewer passengers than 
what a bus would do on the average and, therefore, it would catch up with the 
delayed bus from behind. The slowly moving delayed bus would slow down 
the buses behind it thereby, eventually, creating a jam. In other words, once 
a larger-than-average gap opens up between two successive buses, the gap is 
likely grow further and the steady state in a flnite system would consist of a 
single jam of buses and one large gap. This is very similar to the Bose-Einstein- 
condensation-like phenomenon we have observed earlier in particle-hopping 
models with slow impurities. On the basis of heuristic arguments and mean- 
field approximation it has been argued [302] that this model exhibits a true 
phase transition from an inhomogeneous low-density phase to a homogeneous 
(but congested) high-density phase only in the limit A — > 0. Finally we mention 
that the BRM with parallel dynamics has been recently been studied in [303] 
where also its connection with the NaSch model has been elucidated. 
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13.6 Mobile directional impurities 



We have considered the effects of random hopping probabihties, assigned either 
to the lattice sites [214,137,190] or to the particles [204-206] in the TASEP 
and in the NaSch-typc models of vehicular traffic, on the nature of the corre- 
sponding steady states as well as their approach to the steady-states starting 
from random initial conditions. Toroczkai and Zia [304] solved analytically a 
model with one "mobile directional impurity"; this model is also an extension 
of the ASEP. In this model, N particles, labeled by integers ItoN (from left 
to right), occupy the sites of a one-dimensional lattice of length N where 
periodic boundary conditions are applied. Thus, there is a single "hole" (i.e., 
empty site) in this model. The shifting of the hole from the site in between the 
particles n and n-|-l to the site in between particles n-\-l and n-|-2 is described 
by the statement "hole jump from position n to n -|- 1". In the absence of any 
impurity, the hole at position n can exchange position with either the particle 
on its left (with probability Wn-i,n) or the particle on its right (with proba- 
bility Wn+i^n)- These probabilities are arbitrary and direction-dependent (i.e., 
in general, VTn-i.n 7^ W^n+i,n) but time- independent. Note that the hopping 
probabilities of the hole is determined by the particles (rather than the lattice 
sites) in front and behind it. The general case, where Wn-i,n 7^ l^n+i,n and 
the probabilities Wn~i,n as well as Wn+i,n for different n are chosen randomly, 
is referred to as the random asymmetric case. As the hole wanders, the string 
of particles also shifts as a whole. However, this system is translationally in- 
variant in the sense that the jump rate of any particular particle-hole pair is 
independent of its location on the lattice. 

The "directional impurity" is introduced by identifying a specific bond (whose 
position is fixed with respect to the lattice) as a defect bond such that the 
time-independent rates of particle-hole exchanges across it are fixed at, say, q 
and q' irrespective of the particle-hole pair involved. In other words, when the 
hole is in between the particles n and n + 1, the hopping probability to the 
right (left) is always Wn+i,n (W^n-i,™), except when the defect bond is involved. 
This model can be represented as in Fig. 58 where the defect bond is shown 
as a kink; the motivation for such a kink came from an earlier model [305] of 
gel electrophoresis [306] . 

More specifically, suppose, the probabilities for the exchange of the hole with 
a particle across the defect bond are q if the hole moved upward and q' if 
the hole moved downward, independent of the particle involved. In the so- 
called pure limit, we have, for all n, Wn+i,n — and Wn-i^n — W^. In this 
limit, the hole can also be regarded as a particle undergoing biased diffusion (if 
7^ W^) everywhere except across across a specific defect bond. If we make 
the further assumption that no backward motion of any particle is allowed, 
then this model reduces to the model of TASEP with a single blockage [214]. 
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Fig. 58. (a) The translationally-invariant Toroczkai-Zia model without any defect 
bond, (b) The translational invariance is broken by the defect bond (represented by 
the kink denoted by the dashed hne). 

13.7 Computer networks 



Inspired by the recent success of the methods of statistical mechanics out- 
side the traditional domain of physics, tools of statistical mechanics have also 
been applied to analyze fundamental properties of information traffic on the 
international network of computers (Internet) [307-310]. Messages in the form 
of information packets are continuously being emitted from the hundreds of 
millions of host computers and transported to their destination computers 
through this network. Each of these packets is relayed through the so-called 
routers on its way. The routers can deal with the packets one by one. Each 
router has a finite buffer where the arriving packets get queued up and for- 
warded one by one from the head of the remaining queue to their respective 
next destinations. Since packets run with the velocity of light through the 
cables, information congestion does not take place inside the cables. It is the 
routers which give rise to the information congestion on the internet. Measur- 
ing the fiuctuations in the round-trip time taken by a message on the internet 
(using the ping command of the UNIX operating system), l//-like power spec- 
trum has been observed [307,308]. 

In the square lattice model of a computer network developed by Ohira and 
Sawatari [310] information packets are generated at the sites on the boundary 
at a rate A with the corresponding destination addresses chosen randomly from 
among the boundary sites. The packets can form queues of unlimited length 
at the the inner nodes, which act as routers of the network. At every time 
step, the packets from the heads of the queues at the routers are forwarded 
to the tail ends of the queues at the next router. Both deterministic and 
probabilistic strategies have been considered for selecting the next router to 
which the individual packets are to be forwarded. On reaching their individual 
destinations the packets die. The average number of time steps between the 
birth and death of a packet is referred to as the average lifetime of a packet. 
Computing the average lifetime as a function of the birth rate A of the packets, 
Ohira and Sawatari [310] observed a transition from a low-congestion phase 
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to a high-congestion phase at a non-zero finite value A*. 



14 Summciry and conclusion 

In this section we summarize our conclusions regarding the current status of 
understanding of the statistical physics of vehicular traffic. We also speculate 
on the future trends of research in this area. 

As we stated in the beginnning, one of the main aims of basic research on 
vehicular traffic, from the point of view of statistical physics, is to understand 
the nature of the steady-states of the system. We have summarized the empiri- 
cal evidences available at present in support of the occurrence of three distinct 
dynamical phases, namely, those corresponding to free-flow, synchronized flow 
and stop- and- go traffic. Our critical review of the theoretical works has made 
it clear that, at present, the physical mechanisms at the "microscopic" level, 
which give rise to the synchronized traffic, are not as well understood as those 
responsible for the free-flow and stop-and-go traffic. There arc, however, strong 
indications that for a complete theory, which would account for all these three 
phases, one must take into account not only the vehicles on a given stretch of 
the highway but also on the on- and off-ramps. 

We have explored the possibility of transitions from one dynamical phase to 
another in the NaSch model (and TASEP) with periodic boundary conditions. 
Moreover, we have also presented the generic phase diagram of the TASEP 
with open boundary conditions and explained the notion of boundary-induced 
phase transitions in such one-dimensional driven- diffusive lattice gases which 
are far from equilibrium. Furthermore, we have found that while some models 
exhibit first-order phase transitions, some others exhibit second order phase 
transitions and the signatures of "criticality" while in some rare situations, 
e.g., in the cruise-control limit of the NaSch model, the system is found to 
exhibit even "self-organized" criticality. 

In the beginning, we stated that one of the aims of basic research on vehicular 
traffic is to understand the nature of the dynamical fluctuations around the 
steady-states. The time-dependent correlations functions and the distributions 
of the relaxation times have been computed for the NaSch model and some 
other models but the general questions of the validity of dynamic scaling and 
dynamic universality classes have not been addressed. Another aim of the 
statistical mechanical approach to vehicular traffic, as we stated also in the 
beginning, is to investigate how the system evolves from initial states which 
are far from the corresponding steady-state. The phenomenon of "coarsening" 
of the platoons of vehicles during evolution from random initial states have 
been studied in some models. But, the questions of "universality", if any, of 
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the growth exponents have not been addressed so far. Mctastable states have 
been observed in several CA models. But, to our knowledge, the mechanisms 
of spontaneous decay of such states (analogue of homogeneous nucleation) 
has not been investigated so far. Besides, to our knowledge, so far it has not 
been possible to develope any powerfuU analytical technique for calculating 
the dynamical properties of the traffic models. 

While stating the aims of basic research on vehicular traffic, we also mentioned 
the need to understand the effects of quenched disorder on the steady-states 
as well as on the dynamical properties of the systems. We have seen that the 
randomization of the hopping probability of the vehicles can lead to some 
exotic platooning phenomena which are close analogs of the "Bose-Einstein- 
condensation" . An alternative prescription for introducing quenched disorder 
into the traffic models is to install random bottlenecks on the road and assign a 
time-independent hopping probability (or, equivalently, a "permeability") for 
hopping across bonds in such locations. It has been found that such localized 
bottlenecks can lead to global phase-segregation. 

The NaSch model is the most extensively studied minimal CA model of vehic- 
ular traffic on idealized single-lane highways. We have explained the concep- 
tual framework, and illustrated the use of the mathematical formalism, of the 
cluster-theoretic analytical calculations for the NaSch model. This formalism, 
which yields exact results for the NaSch model with Vmax = 1, gives quite 
accurrate estimates of various quantities of interest as long as Vmax is not too 
large. It would be desirable to develope a new formalism to carry out exact 
analytical calculations for higher velocities too. Here also the limit Vmax — 
is interesting since it shows a rather peculiar behaviour. 

In addition to the detailed discussions on the NaSch model and its various gen- 
eralizations, we have also mentioned briefly some other alternative CA models 
of single-lane highway traffic, e.g., the VDR model, the Fukui-Ishibashi model, 
the Werth-Froese-Wolf model, etc. so that one can appreciate the ongoing ef- 
forts to formulate the most satisfactory minimal model. 

Prom the point of view of practical applications, modelling vehicular traffic 
on multi-lane highways are more relevant than that on idealized single-lane 
highways which are, nevertheless, interesting from the point of fundamental 
understanding of truly non- equilibrium phenomena in drivcn-diffusive lattice 
gases. At present, there are several different alternative prescriptions for formu- 
lating the CA rules for lane-changing of the vehicles on multi-lane highways. 
But, in order to pick out the most appropriate one from among these CA 
models theorists would require input from careful further observations of the 
phenomenon of lane-changing on real multi-lane highways. Empirical observa- 
tions may also indicate modifications or extensions of the CA rules necessary 
for more realistic modelling of the multi-lane traffic. 
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The generalizations of the CA models of traffic on idealized single-lane high- 
ways to those on multi-lane highways may be regarded as extensions of one- 
dimensional model chains to one-dimensional strips. The BML model of ve- 
hicular traffic in cities may be regarded as a further generalizations of these 
models from one-dimensional chains to two-dimensional lattices, or further, 
to decorated lattices. A few different CA rules have been considered so far for 
taking into account the effects of the traffic lights at the crossings of streets in 
such idealized street networks. We have emphasized the intrinsic differences 
between percolation clusters and the cluster of jammed vehicles in the BML 
(and similar) models in spite of some apparent similarities between them. 

We have focused attention mainly on the progress made in the recent years 
using "particle-hopping" models, formulated in terms of cellular automata, 
and compared these with several other similar systems. Although this may be 
a slightly biased overview (as all reviews usually are) of the theory of vehic- 
ular traffic, we have also discussed the main ideas behind all the major ap- 
proaches including the fluid-dynamical, gas-kinetic and car-following theories 
of vehicular traffic. At present the relationships between different approaches 
of modelling have not explored in great detail. It would be very useful if the 
phenomenological parameters of the macroscopic theories can be estimated by 
utilizing the mathematical formulae relating these with those of the "micro- 
scopic" models. 

It is now quite clear that, in order to make significant further progress, we 
not only need more realistic models and better techniques of calculation but 
we also need more detailed and accurate empirical data from real traffic on 
highways as well as more careful re-analysis of the existing data in the light 
of recently developed concepts. So far as the observation of the real traffic is 
concerned, a lot can be learnt from a systematic analysis of aereal pictures 
or video photographs. Alternatively, as the second best choice, a series of 
counting loops along the highway can give more insight, by providing detailed 
information on, for example, time-headways (flux), velocity and local density. 

If you are a critical thinker (or a pragmatist) you may ask: "armed with the 
theoretical tools at our disposal now, can we predict the occurrence of a traffic 
jam at a specific place on a given highway (or street) at a particular instant of 
time"? This question sounds similar to questions often asked in the context of 
some other interdisciplinary topics of current research in the area of complex 
systems, e.g., "can we predict an earthquake", or, "can we predict a stock 
market crash"? Of course, we know that, at present, the best we can hope 
for is to predict (if at all possible) probabilities of occurrences of all these 
phenomena. But, we must admit that, we have a long way to go before we 
come even close to this goal. Nevertheless, we hope you have enjoyed the 
fascinating twists and turns of the way we have covered so far. Our endeveour 
will be more successful if your interest has been stimulated by the intellectual 
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challenges posed by the open problems and if you are willing to uncover the 
current mysteries as well as anticipating new surprises that may lie ahead. We 
are just at the beginning of a long road! 



"The volume of vehicular traffic in the past several years has rapidly out- 
stripped the capacities of the nation's highways. It has become increasingly 
necessary to understand the dynamics of traffic flow and obtain a mathemati- 
cal description of the process" - H. Greenberg (1959) 
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A Definition of update orders 



A dynamical model is not fully defined just by its local transition rules. In 
addition one has to specify the order in which the rules are applied to the 
different particles, i.e. the update ordering (sometimes also called 'dynamics'). 
This is an essential part of the definition of the model since the transient and 
even the stationary state may differ dramatically [311,312]. 

For the NaSch model one uses a parallel update scheme [140] where the rules are 
applied to all particles (i.e., vehicles) at the same time. This kind of ordering 
is sometimes also called synchronous updating. 

Among the various types of asynchronous update schemes most frequently the 
the so-called random- sequential update is used. Here one picks one particle at 
random and applies the transition rules to it. Then one makes another random 
choice (which can also be the same particle again) and so on. This update is 
sometimes called continuous since it can be described by a master equation in 
continuous time. 

Apart from the parallel update there are other updates which are discrete 
in time. We just mention the ordered- sequential updates. Here one starts by 
applying the transition rules to one particle. After that the rules are applied 
to the other particles in a fixed order, e.g. one might continue with the next 
particle ahead of the first one (forward-ordered) or the next particle behind it 
[backward- ordered). We would like to point out that, in principle, one has to 
distinguish two different types of ordered-sequential updates which one could 
name site- ordered- sequential and particle- ordered- sequential^ respectively. In 
contrast to the particle-ordered-sequential update described above, in the site- 
ordered-sequential update the rules are applied to all sites consecutively. This 
might have a strong effect, since a particle might move to next cell ahead 
which then is updated next (for the forward-particle-ordered-sequential up- 
date). Then this particle might move again and so on. This is different from 
the particle-ordered-sequential case where a particle at most moves once dur- 
ing a sweep through the lattice. As an example consider the extreme case 
of the NaSch mode with v^ax — 1 where only one particle is present which 
moves with probability g = 1 to an empty cell in front. This particle will 
move through the whole lattice during one sweep! By looking at a lattice with 
two particles, one can already see that the two different updates might in- 
troduce rather different correlations. Starting with particles separated by d 
empty sites, in the site-ordered-sequential update the left particle will move 
to the right until it reaches the right particle, which then starts to move. On 
the other hand, in the case of particle-ordered-sequential update the particles 
will stay always d oi d — 1 sites apart. For general values of q the situation is 
similar. 
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There are several other updates which can be defined. We refer to the htera- 
ture (see e.g. [166,311-315]) for a comparison of different update procedures. 
The parallel update usually produces the strongest correlations and is used for 
traffic simulations [140] . Note that the forward-particle-ordered-sequential up- 
date is almost identical to the parallel update. In the case of periodic boundary 
conditions a difference only occurs during the update of the last particle. In 
the forward- ordered case the particle in front of it (i.e. the first particle) might 
already have moved since it has been updated earlier. Although this difference 
appears to be minor it can have a large effect. The difference between parallel 
and forward-particle-ordered update can be viewed as a dynamical defect. 



B TASEP 



This simple model of driven systems of interacting particles is one of the most 
exhaustively studied prototype models in nonequilibrium statistical mechan- 
ics [14-16,316,317]. This model can be divided into four classes on the basis 
of the boundary conditions and the update scheme for the implementation of 
the dynamics. In this appendix we consider the TASEP with only random- 
sequential dynamics. 

Let us consider the TASEP with periodic boundary conditions and random- 
sequential dynamics. Since only two states, namely empty and occupied, are 
allowed for each site we can use a two-state variable rii to denote the state 
of the i-th site where = if the i-th site is empty and rii — 1 ii the i- 
th site is occupied. For any given initial configuration {ni(0)}, we can write 
the equations governing the time evolution of (jii{t)) (and all the correlation 
functions) by taking into account all the processes during the elementary time 
interval dt. It is not difficult to establish that 



ni{t -\- dt) 



ni{t) with probability 1 — 2dt 

rii-iit) + ni{t) — ni-i{t)ni{t) with probability dt 
ni{t)ni+i{t) with probability dt 



(B.l) 



-r^ ^ Ui-iil - rii) - ni{l - rii+i) (B.2) 
dt 

which, upon averaging over the history between times and t, leads to the 
equation 

= (rii-i) - {ui) - {ni-iUi) {niUi+i) (B.3) 
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for (nj), the average occupation of the i-th site. Note that the equation for 
{rii) involves two-site correlations. Similarly, it is straightforward to see that 
the equations for the two-site correlations involve three-site correlations and 
so on. Thus, the problem is an intrinsically A^-body problem! The probability 
distribution for this system in the steady-state is given by [22] 

N\{L-N)\ 

Psteady-stateyV^is) — 

where L and refer to the total number of sites and the total number of 
particles, respectively. Prom this distribution it follows that {ni) — N/L and 
('^) = iJzY which lead to {rii) — c and (v) — {1 — c) in the thermodynamic 
limit. 

The stationary state of the TASEP with open boundary conditions and random- 
sequential dynamics has been determined exactly using the so-called Matrix 
Product Ansatz (MPA) (see Appendix P for a more technical introduction) in 
[167] and in [168] using recursion relations. This solution has been generalized 
to different types of discrete dynamics in [172,173,166]. The solution for par- 
allel dynamics was obtained recently in [164] and [165] using generalizations 
of the MPA technique. 



C Naive site-oriented mean-field treatment of the NaSch model 

Suppose, Cy{i,t) = Probability that there is a vehicle with speed v {v — 
0,1,2, ...,v,nax) at the site i at the time step t. Then, obviously, c{i,t) — 
X^J^o"" Cjih t) = Probability that the site i is occupied by a vehicle at the time 
step t and d{i, t) = 1 — c{i, t) is the corresponding probability that the site i 
is empty at the time step t. Using the definition 

v=l 

for the flux J{c,p) one can get the mean-field fundamental diagram for the 
given p provided one can get Cy in the mean-field approximation. 

Step I: Acceleration stage (t — ti) 

(C.l) 

(0 < ^; < Vmax) (C.2) 
iii,t) (C.3) 

Step II: Deceleration stage {ti ^2) 



Co(i,ti) = 0, 
Cy{i,tx) = Cy_i{i,t), 
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Co{i, t2) = Co{i, h) + c{i + 1, ti) Cy{i, ti) 



v=l 



Cy{i,t2) = c{i + v + l,ti)'[[d{i+j,ti) 

j=l v'=v+l 



'^max 



(0 < f < V„ 



(C.4) 



(C.5) 
(C.6) 



Step III: Randomization stage (^2 — ^s) 



Cy{i,t3)^qcy{i,t2) +pcv+i{i,t2), 



{0 <V < Vmax) 



(C.7) 
(C.8) 
(C.9) 



Step IV: Movement stage (^3 — > t + 1) 

Cy{i,t + 1) = c^{i - v,t3), 



{0<V < Vmax) 



(C.IO) 



In the special case Vmax = 1 the equations (C.l-C.lO) get simphfied and, hence, 
we get 



co(i, t + l) = c{i, t)c{i + 1, t) + pc{i, t)d{i + 1, t), 
Ci(i, t + 1) = qc{i — \,t)d{i, t). 

Similarly, in the case of Vmax — 2, one gets 



(C.ll) 
(C.12) 



Co(i,i + 1) = [c{i,t) +pd(i,i)]co(i,i) + [1 ^ pd{i,t)\c{i,t)[ci{i,t) + C2(i,t)], 

(C.13) 

ci(i,i + l) = d{i,t)qcQ{i,t) + d{i,t)[qc{i,t) pd{i,t)][ci{i,t) + C2(i,t)], 

(C.14) 

C2{ht + l) = qd'^{i,t)[ci{i,t)+C2{i,t)]. (C.15) 



D Pciradisical mean-field theory 



For Vmax = 1 the question, whether a state is a GoE state or not, can be decided 
locally by investigating just nearest-neighbour configurations. By analysing 
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the update rules one finds that all states containing the local configurations 
(0, 1) or (1, 1), i.e. configurations where a moving vehicle is directly followed 
by another car, are GoE states. This is not possible as can be seen by looking 
at the previous configurations. The momentary velocity gives the number of 
cells that the car moved in the previous timestep. In both configurations the 
first car moved one cell. Therefore it is immediately clear that (0, 1) is a GoE 
state since otherwise there would have been a doubly-occupied cell before 
the last timestep. The configuration (1, 1) is also not possible since both cars 
must have occupied neighbouring cells before the last timestep too. Therefore, 
according to step 2, the second car could not move. 

The MET equations (64) and (65) have to be modified to take into account the 
existence of GoE states. In general, one has to follow the procedure outlined in 
Appendix C. A quicker way to derive the paradisical mean-field (pME) equa- 
tions is to analyse the MF equations (64) and (65). In (64) the contribution 
c{i;t)c{i + l]t) appears. Since we know that site i + 1 can never be occupied 
by a car with velocity 1 if site i is not empty, this contribution has to be mod- 
ified to c{i;t)co{i + l;t) in pMET. All other contributions are left unchanged 
compared to MET. 

Due to this modification and the corresponding reduction of the configuration 
space the normalization Cq + Ci = c is no longer satisfied automatically. There- 
fore a normalization constant J\f has to be introduced. The final equations for 
a homogeneous stationary state are than given by 



co = Mico+pd)c, (D.l) 
ci=j\fqcd, (D.2) 



with the normalization 



^f^-—,. (D.3) 

Co + a 



Since Cq + Ci = c we have only one independent variable for fixed density c, 
e.g. Ci. Solving (D.l), (D.2) for ci we obtain 

ci = i (l - ^l-4g(l-c)c) . (D.4) 

The fiow is given by /(c) = ci and we recover the exact solution for the case 



In the case Vmax = 2 more GoE states exist. In order to identify these it is 
helpful to note that the rules step 1 - step 4 imply dj{t) = dj{t — l)+Vj^i{t) — 
Vj{t) and therefore 
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d,{t)>v,+,{t)-vjit), (D.5) 
vj{t)<dj{t-l) (D.6) 

The second inequality (D.6) is a consequence of step 2. 

In the following we list the elementary GoE states, i.e. the local configurations 
which are dynamically forbidden (cars move from left to right): 



(0,1), (0,2), (1,2), (0,.,2), (D.7) 
(1,1), (2,1), (2,2), (1,.,2), (2, .,2), (D.8) 

(0,.,.,2). (D.9) 

The numbers give the velocity of a vehicle in an occupied cell and • denotes 
an empty cell. 

The elementary GoE states in (D.7) violate the inequality (D.5), and the 
configurations in (D.8) violate (D.6). The state in (D.9) is a second order 
GoE state. Going one step back in time leads to a first order GoE state since 
(0, •, •, 2) must have evolved from (0, v) (with v — 1 or v — 2). 

Again we can derive the pMF equations by modifying the method for the 
derivation of the MET. Taking into account only the first order GoE states 
(D.7) and (D.8) one obtains the following pMF equations: 

co=M[coc + pd{co + cic)], (D.IO) 
Cl ^Af [pd'^ici + C2) + qd{co + Cic)] , (D.ll) 
C2=Mqd\ci + C2). (D.12) 

The normalization Af ensures co + ci + C2 — c and is given explicitly by 

Co + dci + d'^C2 Co + d{l-C2)' (D.13) 



These equations have been analysed in [142]. After expressing C2 through cq 
and c by 

+ d-^{co + dy-Aqd^{c-Co)y (D.14) 



''^2d^'' 

Inserting this result into (D.IO) we obtain a cubic equation which determines 
Cq in terms of the parameters c and p [142]. Results for different values of p 
are shown in Fig. 19. These results are only slightly modified when also the 
second order GoE state is taken into account [142]. 
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E Equations of car-oriented theory of NaSch Model and COMF 
approximation 



In terms of p, g, g and g{t) = 1 — g{t), the equations describing the time 
evolution of the probabihties Pn{t) for the NaSch model with Vmax — 1 are 
given by 



Po{t+l)=g{t)[Poit) + qPiit)], (E.l) 
Pi(t + l)=git)Poit) + [qg{t)+pg{t)]Piit) + qgimit), (E.2) 
Pn{t + 1) =pg{t)Pn-i{t) + [qg{t) + pg{t)]Pn{t) + qg{t)Pn+i{t) {u > 2). 

(E.3) 



It is worth mentioning here that, for the NaSch model with v^ax = Ij the 
2-cluster probabilities P2{ai,aj) of the the SOMF theory are related to the 
probabilities P„ of the COMF theory through 



P2(l,l) = cPo, 
P2(1,0) = C(1-C)P1, 

P2(0,0) = (l-c)%i 



[n > 1) 



(E.4) 
(E.5) 

(E.6) 



F The matrix-product Ansatz for stochastic systems 



For the stochastic systems considered here the time-evolution of the prob- 
ability P{T,t) to find the system in the configuration r = (ti,... ,ti) is 
determined by the master equation. For random-sequential dynamics it has 
the form 

^^t^ = Y^wiT - T)P(f,t) - Y^wiT - f)P(r,t) (F.l) 

T T" 

with transition rates 'w{f r) from state r to state r. Eq. (F.l) can be 
rewritten in the form of a Schrodinger equation in imaginary time [318], 

^\P{t))^-n\p{t)), (F.2) 

with the state vector \P{t)) = Yl^^ P{T,t)\T). The vectors |t) = \ti,...,tl) 
corresponding to the configurations r form an orthonormal basis of the con- 
figuration space. The stochastic Hamiltonian Ti is defined through its matrix 
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elements 



{r\n\T) = -w{f ^ t), {T\n\T) =Y,w{t ^ t) {t^t). (F.3) 

The stationary state of the stochastic process corresponds to the eigenvector 
|Po) olH with eigenvalue 0. 

For discrete-time dynamics the master equation takes the form 

P(t, t+l)=Y.W{T^ T)P{f, t) (F.4) 

where W{f — > r) = w{f — > r) • are transition probabilities. This can be 
rewritten as 

\P{t + l))=T\P{t)). (F.5) 

Here the stationary state corresponds to the eigenvector \Pq) of the transfer 
matrix T with eigenvalue 1. 

A very powerful method for the determination of stationary solutions of the 
master equation is the so-called matrix-product Ansatz (MPA). For a system 
with open boundaries the weights P{t) in the stationary state can be written 
in the form 

Pin, ...,rL) = n {t,D + (1 - T,)E\ \V). (F.6) 

For periodic boundary condition the MPA takes the form 

P(ri, . . . , T^) = Tr + (1 - ^.)^] j • 

For simplicity we have assumed a two-state system where e.g. = corre- 
sponds to an empty cell j and = 1 to an occupied cell. Zi^ is a normalization 
constant that can be calculated as Zl = (iy|C^|y). In (F.6), (F.7) E and D 
are matrices and and \y) are vectors characterizing the boundary condi- 
tions. The explicit form of these quantities has to be determined from the con- 
dition that (F.6) or (F.7) solves the master equation. This leads in general to 
a algebraic relations between the matrices E and D and the boundary vectors 
(W^l and \y). Once one these have been determined one has a simple recipe for 
determining P(ti, . . . , r^): First, translate the configuration ri, . . . , into a 
product of matrices by identifying each empty cell [jj = 0) with a factor E and 
each occupied cell (tj = 1) with D. In that way the configuration 011001 ■ ■ • 
corresponds to the product EDDEED ■■■ = ED'^E'^D The weight of the 
configuration is then just the matrix element with the vectors {w\ and \v). 
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A simple example is the ASEP discussed in Sect. 8.1.1. Here these quantities 
have to satisfy 



pDE = D + E, 
a{W\E={W\, 
(3D\V) = \V). 



(F.8) 
(F.9) 
(F.IO) 



If one is able to find explicit representations for this algebra one can determine 
in principle all expectation values in the stationary state exactly. For (F.8)- 
(F.IO) one can show that all representations are infinite-dimensional [167]. 
Only on the line a + f3 = p one-dimensional representations (with E, D and 
{W\, \V) being real numbers) are possible^H- 

At this point it might appear that the MPA works only in very special cases. 
However, it can be shown that the stationary state of one-dimensional stochas- 
tic processes is generically of matrix-product form [319-321]P^. Even if it is 
not straightforward to find general representations of the resulting algebras, 
one can at least search systematically for finite-dimensional representations 
on special lines in the parameter space of the model. Furthermore, since the 
mathematical structure of the stationary state is known it is sometimes possi- 
ble to derive rather general results. As an example in [320] interesting relations 
between expectation values for ordered- sequential and sublattice-parallel dy- 
namics have been derived. 

For a more detailed description of the MPA for different types of dynamics 
and its relation with the MPA technique for quantum-mechanical spin systems 
[322] we refer to [166] . A review of the treatment of the ASEP using the MPA 
is given in [22]. The MPA has also been extended to treat the full dynamics, 
not only the stationary state [323]. Using time-dependent matrices D{t) and 
E{t) one obtains the Bethe Ansatz equations for the corresponding stochastic 
Hamiltonian [323-325]. 



G Two schemes for solving the mean-field approximation of the 
DTASEP 

A mean-field approximation scheme for this model has also been developed 
[137]. The time-averaged steady-state current Jjj+i in the bond (j, j + 1) 
is given by Jjj+i = qjj+i{nj{l — rij+i)). In the mean-field approximation. 

This can be seen easily from (F.8)--(F.10). 

We have to mention here that these results up to now do not include the case of 
parallel dynamics. 
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{nj{l — rij+i)) — {nj){{l — and, hence, 

J = = - Cj+i) (G.l) 

where Cj — (rij). In order to calculate the steady-state flux J as a function 
of the mean density c of the particles, Tripathy and Barma [137] used two 
different iteration schemes based on the equation (G.l). 

(i) Constant- current iteration scheme: In this scheme, for a given system length 
L and a fixed flux J = Jq, one starts with some value of Ci and, computes all 
the other cj (j > 1) using the equation (G.l), i.e., 

c,+i = l-— j = l,2,...,L (G.2) 

together with the periodic boundary condition Cj+L = Cj. If the iteration 
converges, i.e., one gets all the site densities in the physically acceptable range 
[0, 1], one accepts the average of these final site densities to be the global mean 
density of the particles corresponding to the flux Jq. (ii) Constant- density 
iteration scheme: In this scheme, for a given system length L and fixed global 
average density c, one begins by assigning the site densities < Cj(0) < 1 to 
the lattice sites subject to the global constraint j J2j Cj(0) = c. Then, the site 
densities are updated in parallel according to 

cjit + 1) = cj{t) + Jj^,j{t) - Jj,j+i{t), J = 1, 2, L (G.3) 

which follows from the equation (G.l). It is straightforward to verify that the 
this iteration scheme keeps the average global density c unchanged at every 
step of updating and hence the name. After sufficient number of iterations the 
set of densities converge to a set {cj} and the flux on each bond converge to 
the steady-state flux Jq. 



H Self-consistent equations for Vx and Vy in the mean-field approx- 
imation of the BML model 

Suppose, Vx and Vy denote the average speeds of east-bound and north-bound 
vehicles, respectively. Then, on the average, an east-bound vehicle spends a 
time 1/vx at a site whereas a north-bound vehicle spends a time 1/vy at a 
site. The north-bound vehicles lead to a reduction of the speed of the east- 
bound vehicles by Uy/vy. Moreover, because of the hindrance of the east-bound 
vehicles by other east-bound vehicles ahead of it there will be further reduction 
of the speed of the east-bound vehicles by n^; [— — 1] . Furthermore, if the density 
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Vy 





of the overpasses is /„, then 
Similarly, the corresponding equation for Vy is given by 

= 1 - (1 - fo] 



(H.l) 



(H.2) 



In the special case Hx = Uy = n/2 both the equations (H.l) and (H.2) reduce 
to the form 



1 1 
v~2 



n 



where = Vy = v. The solution of the quadratic equation (H.3) for v is 
1 



(H.3) 



1 H — n + 



\ 



l + i-^n] -4(l-/> 



H.4) 



I Derivation of the equations in the microscopic theory of the BML 
model 



By definition, ni{x,y;t) {n^{x,y : t)) is unity if the site {x,y) is occupied 
at time i by a north-bound (east-bound) vehicle and zero if the site {x, y) 
is not occupied by a north-bound (cast-bound) vehicle. Normalization re- 
quires n^{x, y; t) + n^{x, y;t) = 1 — Tiempt:y('^- IJ] t) where the two-state variable 
^empty{x,y]t) is Unity if, at time t, the site {x,y) is empty and zero if the 
site (x, y) is not empty. In order to describe the state of the signals at time 
t, one also defines a two-state variable S{t): S{t) — 1(0) if the signal is green 
(red) for the vehicles under consideration. The space-average of n-f (x, y; t) and 
n^{x,y;t) are c_>(i) and C|(i), respectively. Besides, the time-average of S{t) 
is 1/2. 

The updating rules of the BML model lead to the equations 



n-^{x,y;t + l) = n^{x,y;t) [n^ix + l,y;t) + n^{x + l,y;t)] S{t) 
+n-^{x - 1, y; t) [1 - n^{x, y; t) - n_(x, y; t)] S{t) + n^{x, y; t)[l - S{t)] 

(I.l) 

The first term on the right hand side of equation (I.l) describes the situation 
when the east-bound vehicle, which was at the site {x, y) at time finds a 
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green signal but cannot move because the next site towards cast is occupied by 
another vehicle. The second term on the right hand side of (1.1) corresponds 
to the situation where the east-bound vehicle, which was at the site {x — l,y) 
at time t, finds a green signal and moves to the next site towards east, which 
was empty. The last term on the right hand side of equation (I.l) arises from 
the possibility that the cast-bound vehicle, which was at (x, y) at time t, could 
not move because of a red signal, irrespective of the state of occupation of the 
next site towards cast. Following the similar arguments, one can also write 
down the corresponding equation for the north-bound vehicles. 
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